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Abstract. Strongly correlated quantum fluids are phases of matter that are 
intrinsically quantum mechanical, and that do not have a simple description in terms 
of weakly interacting quasi-particles. Two systems that have recently attracted a great 
deal of interest are the quark-gluon plasma, a plasma of strongly interacting quarks and 
gluons produced in relativistic heavy ion collisions, and ultracold atomic Fermi gases, 
very dilute clouds of atomic gases confined in optical or magnetic traps. These systems 
differ by more than 20 orders of magnitude in temperature, but they were shown to 
exhibit very similar hydrodynamic flow. In particular, both fluids exhibit a robustly 
low shear viscosity to entropy density ratio which is characteristic of quantum fluids 
described by holographic duality, a mapping from strongly correlated quantum field 
theories to weakly curved higher dimensional classical gravity. This review explores the 
connection between these fields, and it also serves as an introduction to the Focus Issue 
of New Journal of Physics on Strongly Correlated Quantum Fluids: from Ultracold 
Quantum Gases to QCD Plasmas. The presentation is made accessible to the general 
physics reader and includes discussions of the latest research developments in all three 
areas. 
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1. Introduction 

This review covers the convergence between three at first sight disparate fields: ultracold 
quantum gases, quantum chromodynamic (QCD) plasmas, and holographic duality. 
Ultracold quantum gases have opened up new vistas in many-body physics, from novel 
quantum states of matter to quantum computing applications. There are over one 
hundred experiments on ultracold quantum gases around the world on every continent 
but Antarctica. The QCD plasma, also called the quark-gluon plasma (QGP), has 
been the subject of intense experimental investigation over more than two decades, 
continuing now at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National 
Laboratory and the Large Hadron Collider (LHC) at the European Organization for 
Nuclear Research (CERN). A QGP is predicted to have occurred in the first microsecond 
after the Big Bang, and re-creation of the QGP on earth at present times allows us to 
probe the physics of the early universe. Holographic duality is a powerful mapping from 
strongly interacting quantum field theories, where the very concept of a quasiparticle can 
lose meaning, to weakly curved higher dimensional classical gravitational theories. This 
duality provides a new approach for modeling strongly interacting quantum systems, 
yielding fresh insights into previously intractable quantum many-body problems key to 
understanding experiments such as ultracold quantum gases and the QGP. 

What do these three fields have in common? All treat strongly correlated quantum 
fluids. Generically, strong interactions give rise to strong correlations. By strongly 
correlated we mean that we cannot describe a system by working perturbatively from 
non-interacting particles or quasiparticles. In the case of electrons in condensed matter 
systems this means that theories constructed from single-particle properties, like the 
Hartree-Fock approximation, cannot describe a material. In the case of fluids, we 
mean that kinetic theories based on quasiparticle degrees of freedom, in particular 
the Boltzmann equation, fail.* The natural candidates for building quasiparticles are 

* Fluids are materials that obey the equations of hydrodynamics. The word liquid refers to a phase 
of matter that cannot be distinguished from a gas in terms of symmetry, but exhibits short range 
correlations similar to those in a solid, and is separated from the gas phase by a line of first order 
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Figure 1. Temperature and pressure scales of extreme quantum matter. Ultracold 
quantum gases are the coldest matter produced to date, while the quark-gluon plasma 
is the hottest, together spanning about 18 orders of magnitude in temperature and 
more than 40 orders of magnitude in pressure. Yet these systems exhibit very similar 
hydrodynamic behavior, as characterized by the shear viscosity to entropy density ratio 
shown in Fig. 2). We also include two other well known quantum fluids, liquid helium 
and hot proto-neutron star matter, as well as a classical fluid, water, and a classical 
plasma, the Coulomb plasma in the sun. 



quark and gluons in the case of the QGP, neutrons and protons in the case of nuclear 
matter, and atoms in the case of ultracold atomic gases. In strongly interacting systems 
the mean free path for these excitations is comparable to the interparticle spacing, 
and quasiparticles lose their identity. Even though kinetic theory fails, nearly ideal, 
low viscosity hydrodynamics is a very good description of these systems. This is a 
central prediction of holographic duality, and it has been verified experimentally for 
both ultracold quantum gases and the QGP, as we will explore in this Review. 

As shown in Fig. 1, strongly correlated quantum fluids cover an enormous range 
of scales in temperature and pressure. s We remind the reader that temperature T and 

transitions that terminates at a critical endpoint. A plasma is a gas of charged particles. Gases, 
liquids, and plasmas behave as fluids if probed on very long length scales. Weakly coupled systems 
exhibit single particle behavior if probed on microscopic scales, but strongly coupled systems behave 
as fluids also on short scales. Liquids are typically more strongly correlated than gases, and more likely 
to behave as a fluid. 

E The points in Fig. 1 correspond to the range of temperatures for which the transport measurements 
shown in Fig. 2 have been performed. For the ultracold atomic Fermi gas experiments described in 
Sec. 2.1 the critical temperature is roughly 500 nK (the exact value depends on the trap geometry 
and the number of particles; Bose gases have been cooled to temperatures below 1 nK). The data 
points for helium and water are centered around the critical endpoint of the liquid gas transition. 
The point for the solar plasma corresponds to the geometric mean of the temperatures in the core 
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Figure 2. Transport properties of strongly correlated fluids. Ratio of shear viscosity g 
to entropy density s as a function of (T — T c ) /T c , where T c is the supcrfhiid transition 
temperature in the case of ultracold Fermi gases, the dcconfinement temperature in the 
case of QCD, and the critical temperature at the endpoint of the liquid gas transition 
in the case of water and helium. The data for water and helium are from [1], the 
ultracold Fermi gas data are from [2], the quark-gluon plasma point (square) is taken 
from the analysis of [3], the lattice QCD data (open squares) from [4], and the lattice 
data for the ultracold Fermi gas (open circles) are the 8 3 data from [5]. The dashed 
curves are theory curves from [6, 7, 8, 9]. The theories are scaled by overall factors 
to match the data near T c . The lines labeled "holographic bounds" correspond to the 
KSS bound h/(4irk B ) [11] and the Gauss-Bonnet bound (16/25)^/(4^) [10]. Similar 
compilations can be found in [11, 12, 13]. 



energy E are equivalent up to a factor of Boltzmann's constant, ks = 1.3806503 x 10 -23 
J/K, with E = ksT . We focus on fluids that can be studied in bulk, as opposed to 
quantum liquids that exist on lattices. We show ultracold Fermi gases, liquid helium, 
neutron matter in proto-neutron stars, and the QGP. For comparison we also show a 
classical fluid, water, and a classical plasma, the Coulomb plasma in the sun. 

Figure 2 shows that despite the large range in scale there is a remarkable universality 
in the transport behavior of strongly correlated quantum fluids. Transport properties of 
the fluid can be characterized in terms of its shear viscosity 77, which governs dissipation 
due to internal friction. A dimensionless measure of dissipative effects is the ratio 77/s 
of shear viscosity to entropy density in units of h/ks- Near the critical point, where 
the role of correlations is expected to be strongest, the ratio r]/s has a minimum. For 
classical fluids the minimum value is much bigger than h/ks, but for strongly correlated 

and the corona. The neutron matter point is at T = lMeV/fce = 1.2 ■ 10 10 K and at a density 
n = 0.1 no, where no = 0.14/fm 3 is nuclear matter saturation density. Neutron stars are born at 
T ~ 10 McV/fce, and they can cool to temperatures below 1 kcV/ks- The critical temperature of the 
QGP is T c ~ 150McV/fcB = 1-75 • 10 12 K. Experiments with heavy ions explore temperatures up to 
~3T C . 
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quantum fluids 77/s is of order H/ks, indicating that dissipation is governed by quantum 
effects. We observe, in particular, that 77/s for the QGP and ultracold Fermi gases 
is quite similar, even though the absolute values of 77 and s differ by many orders of 
magnitude." 

Remarkably, these values of r]/s lie near a lower bound, rj/s > h/(47ikB), which 
arises in the study of 4+1 dimensional black holes in classical Einstein gravity. These 
gravitational theories are conjectured to be dual to certain 3+1 dimensional quantum 
field theories; see Sec. 4. This lower bound is known to be non- universal; it can be 
violated in a more general class of theories dual to a gravitational theory known as Gauss- 
Bonnet gravity. Imposing basic physical requirements such as causality and positivity 
leads to a slightly smaller bound^ rj/s > (16/25)/i/(47r/c B ). The main feature of the 
results obtained using holographic dualities is that, at strong coupling, rj/s is both 
unusually small and relatively insensitive to the precise strength of the interactions, so 
long as they are strong. This is in sharp contrast to the predictions of kinetic theory for 
a weakly interacting gas. As a result, rj/s serves as a probe of the strength of correlations 
in a quantum fluid. 

We have chosen to focus on the fields of ultracold quantum gases and the QGP 
not only for their range of energy and density scales, but also because of their 
broad intrinsic interest. Ultracold fermions are connected to a wide variety of exotic, 
strongly interacting systems in nature, ranging from high-temperature superconductors 
to nuclear matter. They are incredibly flexible many-body systems that allow nearly 
arbitrary tuning of interactions, symmetries, spin structure, effective mass, and imposed 
lattice structures. The QGP on the other hand explores a very different regime from 
other particle physics experiments: thousands of particles are produced, and these 
particles form a hot and dense environment recreating the conditions of the early 
universe. 

Our review is outlined as follows. A key introductory or heuristic plot is presented 
in each of our three fields for the general physics reader. In Sec. 2, Fig. 3 shows the 
phase diagram for the Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein condensate 
(BEC) crossover, midway through which the unitary quantum gas, a strongly correlated 
quantum fluid, is obtained. In Sec. 3 the QCD phase diagram is shown in Fig. 11, 
including the strongly and weakly interacting QGP being explored presently at RHIC 
and the LHC. In Sec. 4, not only is an extensive heuristic description of holographic 
duality provided in Sec. 4.1, but a holographic dictionary is presented in Table 4.2.3, 
mapping quantities in strongly interacting quantum fields similar to the unitary 

" The theoretical curves, as well as the data for helium and water corresponds to systems in the 
thermodynamic limit. The lattice data for the QGP and the ultracold Fermi gas have finite volume 
corrections that have not been fully quantified. The experimental data point for the QGP is based on 
an analysis that assumes rj/s to be temperature independent. The data points for the ultracold Fermi 
gas show the ratio of trap averages of 77 and s. The local value of 77/s at the center of the trap is likely 
to be smaller than the ratio of the averages; see Sec. 5.1. 

^ Whether this value represents a true lower bound, or whether more general classes of fluids with even 
smaller values of 77/s can exist, is an active area of research; see Sec. 4.3. 
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quantum gas or the QGP to their gravitational duals. The three main sections also 
cover all 39 papers from three fields comprising this Focus Issue. Finally, in Sec. 5, 
rather than a summary, we instead conclude with open questions in each field. 



2. Ultracold Quantum Gases 

Ultracold quantum gases provide a unique table-top paradigm for exploring the 
properties of quantum many-body systems in nature [14, 15, 16, 17, 18, 19, 20], from the 
thermodynamics of high temperature superconductors to the hydrodynamics of QGPs. 
These gases are made of mainly alkali metal atoms but more recently other atoms as 
well as diatomic molecules. They can be fermionic or bosonic with a wide variety of 
internal hyperfine spin structures. They can be made strongly or weakly interacting and 
both attractive and repulsive. They are contained in a variety of magnetic and optical 
traps in one, two, and three dimensions, including optical lattices. The latter give rise 
to arbitrary lattice structures. Because these gases are dilute and very cold, they are 
described by first principles theories built up from low-energy binary scattering between 
atoms, and well-known interactions with magnetic and optical fields. The collection of 
atoms can be probed and manipulated by external laser beams and pulses, as well as 
external magnetic fields. 

In this review we focus on atomic Fermi gases [18, 17, 21, 22], in particular strongly 
interacting Fermi gases. Several have been cooled to degeneracy using evaporative 
cooling methods. The most widely studied atoms are 40 K [23] and 6 Li [24, 25, 26, 27]. + 
Experiments are carried out at temperatures in the nanokelvin to microkelvin range, 
with typical densities from 10 11 to 10 14 atoms/cc. For a 6 Li atom at nanokelvin 
temperatures, the thermal de Broglie wavelength 



V vhkbI 

is on the order of several /im. Quantum degeneracy occurs when the de Broglie 
wavelength is greater than or on the order of the particle spacing, AdB^ 3 ^ 1, where 
n is the density; this condition is equivalent to T '/Tp > 1, where Tp is the Fermi 
temperature.* Our interest is in ultracold Fermi gases that are quantum degenerate: 
T/Tp < 1. 

A cloud of trapped dilute fermions is typically about 100 /im in size, and is deformed 
by harmonic trapping fields into prolate or oblate forms, commonly called a cigar or a 
pancake. In the degenerate regime the cloud is stabilized against collapse by Pauli 
pressure [23, 26, 24]. The size of the cloud depends on the interplay between the 
trapping potential, the Pauli pressure, and interactions between the atoms. Because 

+ Degeneracy has also been achieved in metastable 3 He* [28], in 171 Yb and 173 Yb [29], and recently in 
87 Sr [30] and 161 Dy [31]. 

* In this review, Tp is always defined with respect to the non-interacting Fermi gas, Tp = h 2 k F /(2mkB) 
with kp = {'&iT 2 n) 1 l s ' for a two-component gas. 
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Figure 3. Ultracold Fermi gas phase diagram. Sketch of the Bardeen-Cooper- 
Schricffcr (BCS) to Bose-Einstein condensate (BEC) crossover for ultracold Fermi 
gases. When the scattering length a s passes through a pole, so that l/(kpa s ) —t 0, 
one obtains a strongly correlated fluid, the unitary gas. The critical temperature T c 
for the phase transition only approaches the pairing temperature T pa i r in the limit 
l/(fcpa) — > — oo. The crossover region is the strongly interacting regime, loosely 
defined by |l/(fepa s )| < 1. Note that we denote the scattering length by a in the text. 
Used with permission from Rcf. [33]. 



of the low density and ultracold temperatures these interactions are dominated by an 
effective s-wave contact interaction. The scattering amplitude is of the form 

f{k)= -l/a + roP/2-ifc' (2) 

where a is the s-wave scattering length and r is the effective range. Higher partial waves 
as well as short range corrections are suppressed by powers of ro/AdB and r^n 1 ^) The 
scattering length is widely tunable by a Feshbach resonance [32], an external magnetic 
field that brings a weakly bound excited molecular state into resonance with the unbound 
atomic scattering state. 

Each of the different trapped atomic elements used in ultracold quantum gas 
experiments has an internal spin structure due to hyperfine structure of the atom, that 

" The range of the atomic potential is on the order of the van der Waals length I = (mCs/ll 2 ) 1 / 4 , where 
Cq controls the van der Waals tail of the atomic potential, V ~ Ce/r e . We assume that the p-wave 
scattering length is natural, meaning a p ~ Tq. 
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is, the combination of the nuclear spin and, in the case of the alkali metals, the electron 
outside the closed shell. For instance, 6 Li has a nuclear spin of 1 and one unpaired 
electron. The two lowest hyperfme states have a total spin of 1/2, and the remaining 
four have a total spin of 3/2. By selecting out hyperfme states through appropriate 
laser-induced transitions and trapping and cooling methods, experiments can thus work 
with a variety of spin structures. The case in which two hyperfme states are trapped 
is effectively equivalent to a spin 1/2 atom. Tuning the scattering length by using a 
Feshbach resonance, one obtains three distinct regimes, shown in Fig. 3. The first is 
for weak attractive scattering, —kpa <C 1, where kp is the Fermi momentum. Then 
for temperatures well below the Fermi temperature Tp one obtains a BCS state [34], or 
s-wave superconductivity. We call this an atomic Fermi superfluid, since our systems 
are in fact neutral. In such a state fermions of opposite spin join to make Cooper pairs, 
but their average pair size £ c is greater than the interparticle spacing n -1 / 3 , so that they 
are overlapping: ^n 1 / 3 ^> 1. Tuning a as in the phase diagram, we observe that the 
scattering length passes through a pole; note that the figure shows temperature as a 
function of the inverse scattering length, l/{kpa). Thus as a — > ±oo, l/(kpa) — > 0, and 
one obtains a second regime, called the unitary gas. This middle regime is a strongly 
correlated fluid, and one finds ^n 1 / 3 ~ 1, i.e., the pair size is about equal to the 
interparticle spacing. Finally, for large positive scattering lengths, the paired fermions 
make much more tightly bound molecules, and one obtains a molecular BEC, similar 
to the well-known atomic BECs. This regime is depicted on the far right of Fig. 3. In 
practice these molecules are still quite large, thousands of Bohr or more, but still much 
smaller than the interparticle spacing, so that ^n 1 / 3 < 1. 

The upper curve on the figure depicts the pair formation temperature T pa ; r , which 
in general is distinct from the critical temperature for superfluidity, T c [35]. Note that 
superfluidity is associated with the spontaneous breakdown of a global symmetry, the 
U(l) phase symmetry of the wave function, and T c is therefore always well defined. 
T pair , on the other hand, is not associated with a symmetry or a local order parameter 
and may not be well defined. This remark is particularly relevant in the BCS regime, 
where the size of the pairs is large compared to the inter-particle spacing. Physically 
we expect that in the BCS regime there are no pre-formed pairs, and T c and T pa i r are 
very close together. 

Although we refer to these systems as ultracold, in terms of the dimensionless ratio 
T/Tp, and in comparison to solid state systems, they are quite hot. In the unitary 
regime the phase transition occurs at T C /T F = 0.167(23) [36], compared to typical solid 
state superconductors in which T c /Tp < 0.01. The unitary Fermi gas is a very high 
T c superfluid. As indicated in Fig. 3, in the BCS regime the temperature required to 
achieve a phase transition to an atomic Fermi superfluid is quite low. In this regime, 
the critical temperature is given by the weak coupling expression [37] 




(3) 
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where 7 ~ 0.577 is the Euler constant. The numerical value of the factor in front of 
the exponent is 0.277. We observe that even though Eq. (3) is formally valid only in 
the limit fcjr|a| <C 1 it also provides a reasonable estimate for T c at unitarity. This is 
somewhat of an accident, because higher order corrections in kp\a\ are divergent in the 
unitary limit. 

In the following, we explore the unitary regime of the BCS-BEC crossover for 
ultracold Fermi gases. In Sec. 2.1 we present an overview of experiments on these 
systems. In Sec. 2.2 we focus on universal aspects of unitary gases. In the strongly 
correlated regime the scattering length diverges, and the remaining length scales in the 
problem are the Fermi length \jkp and the de Broglie wave length AdB, given in Eq. (1). 
Thus many theoretical statements can be made despite the lack of a small parameter 
or perturbative calculations. In Sec. 2.3 and 2.4 we treat the thermodynamics and the 
structure of the phase diagram for unitary gases, and in Sec. 2.5 we focus on transport 
properties. Section 2.6 presents an overview of ultracold Fermi gases in optical lattices. 
Finally, in Sec. 2.7 we treat new directions in unitary gases as presented in this Focus 
Issue, including novel experimental probes, solitons, imbalanced systems and polarons, 
disorder, quantum phase transitions, Efimov physics, and the use of three hyperfme 
states to explore SU(3) physics and connections to the QGP. 

2.1. Ultracold Fermi Gas Experiments 

Historically, atomic Fermi gases were first brought to degeneracy at JILA in 1999 [23], 
using a mixture of two hyperfme states in 40 K to enable s-wave scattering in a magnetic 
trap. Dual species radio-frequency-induced evaporation produced a degenerate, weakly 
interacting sample, with T jTp ~ 0.5. Later, degeneracy was achieved in fermionic 6 Li by 
direct evaporation from a MOT-loaded optical trap [26] and by sympathetic cooling with 
another species [24, 25, 27], producing a lower T jTp. However, the minimum reduced 
temperature was initially limited to T jTp ~ 0.2, which may have been a consequence of 
trap-noise-induced heating [40] or, at the lowest temperatures, Fermi hole heating [38] 
in combination with evaporative cooling [39]. 

Optical traps enabled a dramatic improvement in the efficiency of evaporation and 
the creation of strongly interacting Fermi gases, through the use of magnetically-tunable 
collisional resonances, or Feshbach resonances [41]. Feshbach resonances in fermionic 
atoms were initially characterized in 2002 by several groups [42, 43, 44]. For a recent 
review of Feshbach resonances see Ref. [32]. In a Feshbach resonance, a bias magnetic 
field tunes the total energy of a pair of colliding atoms in the incoming open (triplet) 
channel into resonance with a molecular bound state in an energetically closed (singlet) 
channel. At resonance, the zero-energy s-wave scattering length a diverges and the 
collision cross section attains the unitary limit, proportional to the square of the de 
Broglie wavelength, i.e., a = Atr/k 2 , where hk is the relative momentum. The collision 
cross section therefore increases with decreasing temperature, enabling highly efficient 
evaporative cooling in optical traps and much lower reduced temperatures T jTp ~ 0.05. 
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Figure 4. Experimental images. Elliptic flow of a strongly-interacting Fermi gas as 
a function of time after release from a cigar-shaped optical trap: from top to bottom, 
100 (is to 2 ms after release. The pressure gradient is much larger in the initially 
narrow directions of the cloud than in the long direction, causing the gas to expand 
much more rapidly along the initially narrow directions, inverting the aspect ratio. 
Achieving nearly perfect elliptic flow requires extremely low shear viscosity, as is the 
case also for a quark-gluon plasma. The sequence of images is created by recreating 
similar initial conditions and destructively imaging the cloud at different times after 
release [46]. 

An optical trap is formed by a focused laser beam. Atoms are polarized by the field 
and attracted to the high intensity region, when the laser is detuned below resonance 
with the resonant optical transition, so that the induced dipole moment is in phase with 
the field. For large detunings, obtained using infrared lasers, the trapping potential is 
independent of the atomic hyperfine state, enabling several species to be trapped, which 
is ideal for Fermi gases [45]. Forced evaporation is accomplished by slowly lowering the 
intensity of the optical trap laser beam. Near a Feshbach resonance, a highly degenerate 
sample can be produced in a few seconds [46] . 

A milestone in the Fermi gas field was the observation in 2002 of a strongly 
interacting, degenerate Fermi gas [46], in the so-called BEC-BCS crossover regime, using 
this method. In contrast to Bose gases, which are unstable and undergo three-body loss 
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on millisecond time scales near Feshbach resonances, the two-component Fermi gas was 
found to be stable, as the Pauli principle suppresses three-body s-wave scattering [47]. 
Released from the cigar-shaped trapping potential of the focused beam, the cloud was 
observed to expand much more rapidly in the initially narrow direction, compared to 
the long direction, as a consequence of the much larger pressure gradient along the 
narrow axis. Consequently, the aspect ratio inverts from a cigar to an ellipse, as shown 
in Fig. 4. Remarkably, the same type of elliptic flow is also observed in the momentum 
distribution of an expanding quark-gluon-plasma produced in an off-center collision of 
two heavy ions; see Sec. 3.5. There, the temperature is nineteen orders of magnitude 
hotter and the particle density is twenty-five orders of magnitude greater than that of 
the Fermi gas. In both systems, however, this nearly perfect "elliptic" flow, Fig. 4, is 
a consequence of extremely low viscosity hydrodynamics, which persists in the normal, 
non-superfluid unitary gas and deeply connects these two apparently disparate fields. 

The creation of a degenerate Fermi gas near a Feshbach resonance was followed in 
2003 by the first measurements of the interaction energy [48, 49] and the creation of Bose- 
condensed dimer molecules [50, 51, 52, 53]. In 2004, condensed fermionic atom pairs 
were observed using a fast magnetic field sweep to project the pairs onto stable molecular 
dimers [54, 55]. Using this method, the first phase diagram in the crossover region was 
obtained (see Fig. 3 below) as a function of magnetic field and temperature, albeit using 
the temperature of the ideal gas obtained by an adiabatic sweep to a non-interacting 
regime [54]. Evidence of superfluidity in a Fermi gas was provided by measurements 
of collective mode frequencies and damping rates versus temperature [56] and magnetic 
field [57], and measurement of the pairing gap by radio-frequency spectroscopy [58]. The 
observation of a vortex lattice in 2005 provided direct verification of Fermi superfluid 
behavior [59]. 

Also in 2005, initial thermodynamic measurements were done by adding a controlled 
amount of energy to the cloud and measuring an empirical temperature from the 
corresponding cloud spatial profile [60]. However, the results were model-dependent, as 
calibration of the empirical temperature relied on comparing the measured cloud profiles 
with theoretical predictions. Model-independent measurements were soon to follow, 
based on universal behavior in the unitary regime, where the local thermodynamic 
quantities, such as the pressure, are functions only of the density n and temperature 
T [61]. 

Model-independent measurements of the total energy E of a resonantly interacting 
Fermi gas are based on the Virial theorem, which holds for a unitary gas as a consequence 
of universality and yields the energy directly from the cloud profile [62]. Using 
entropic cooling [63, 64], a model-independent measurement of the total entropy S 
was accomplished by an adiabatic sweep of the bias magnetic field from resonance 
to the weakly interacting regime, where the entropy can be calculated from the 
cloud profile [65]. As T = dE/dS, these measurements enabled the first model- 
independent temperature calibration and estimates of the critical parameters in the 
strongly interacting regime [66]. A refined temperature calibration is used in the 
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Figure 5. Ultracold quantum gas experimental apparatus. Left: sketch of the 
experimental apparatus for ultracold fermions. Right: Apparatus for the Duke 
experiments (currently at North Carolina State University). Compare to sketch of 
QGP experiment at the LHC in Fig. 12: the quantum gas experiment is about ten 
times smaller (2.5 meters vs. 26 meters), but the size of the trapped ultracold gas is 
11 orders of magnitude larger (a few hundred micrometers vs. a few femtometers). 
The ultracold quantum gas is at nanokelvin temperatures, or pico-eV, compared to 
the deconfinement temperature of ~ 2 x 10 12 K in the QGP, or 200 MeV, created by 
colliding gold nuclei at energies of 100 GcV/nucleon. 

measurement of universal quantum viscosity, as described in this Focus Issue [2]. 

Measurements of global thermodynamic quantities from the cloud profiles in the 
strongly-interacting regime are now superseded by model-independent measurements of 
local quantities [67, 68]. Using the Gibbs-Duhem relation 

dP = ndfi (4) 

at constant temperature, the local pressure is obtainable from the integrated column 
density, where the local chemical potential \x is determined by the known trap 
potential [69, 70]. Combined with a temperature measurement, the local equation of 
state P(fi, T) or P(n, T) is determined. The most precise local measurements avoid 
temperature measurement, which introduces the most uncertainty, by determining the 
pressure, density and compressibility from the cloud profiles. The resulting equation 
of state reveals clearly a lambda transition, and provides the best measurement of 
the critical parameters for a unitary Fermi gas [36], as described in detail in Sec. 2.3. 
Measurements of equilibrium thermodynamic quantities are now used as stringent tests 
of predictions, as described in this Focus Issue by Hu [71]. These thermodynamic 
measurements are connected to universal hydrodynamics and transport measurements, 
as described in [2]. 

We proceed to describe the all-optical methods developed at Duke in 2002 [26, 46], 
as one specific example of experimental techniques, which are closely tied to the theme of 
this Focus Issue, viscosity and transport measurements on Fermi gases in the universal 
regime [2]. A degenerate, strongly interacting Fermi gas is readily made by all-optical 
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methods [46]. As sketched in the left panel of Fig. 5, a magneto-optical trap (MOT) is 
used to pre-cool a 50-50 mixture of spin-up and spin-down 6 Li atoms, which is loaded 
into a CO2 laser optical trap and magnetically tuned to an s-wave Feshbach resonance. 
Atoms from the source, (lower right, green cylinder), form an atomic beam (blue arrow) 
that is slowed by radiation pressure forces from a resonant laser beam (top left, opposing 
red arrow). For 6 Li, the deceleration is 2 x 10 6 m/s 2 , slowing the atoms from oven 
thermal speeds of a km/s to tens of m/s sec over a distance of a fraction of a meter. 
Six laser beams (three thick red lines) then propagate toward the center of the MOT 
(point of intersection of three thick red lines), creating inward damping forces that cool 
the atoms. Opposing magnetic fields, created by two coils (stacked black circles, top 
and bottom,) spatially tune the atomic resonance frequency, causing the six beams to 
produce a harmonic restoring force at the MOT center. Typical MOT clouds are a few 
millimeters across and contain several hundred million atoms. A trapping laser beam 
(shown in yellow) is focused (lenses indicated by two light blue ovals) at the MOT 
center and loaded. After turning off the MOT beams and the MOT magnetic field, 
an additional bias magnetic field tunes the atoms to a collisional (Feshbach) resonance. 
Forced evaporation near resonance, by lowering the trap depth, rapidly cools the cloud 
to quantum degeneracy, i.e., T/T F <C 1, producing a cigar-shaped cloud that is typically 
a few microns in diameter and several hundred microns long, containing several hundred 
thousand atoms. In the right panel of Fig. 5 is shown the experimental apparatus from 
the Duke laboratory, currently located at North Carolina State University. From right 
to left in the photo: the oven assembly where hot fermions are produced (aluminum 
housing); camera to produce density images (blue device in foreground); Zeeman slower 
to bring atoms into MOT (middle, behind camera); bias field magnets containing MOT 
in ultrahigh vacuum (white plastic housings); ZnSe lens for the CO2 laser trapping beam 
and optical table (left). 

In the simplest case, the optical trap consists of a single laser beam, focused into the 
center of the MOT and detuned well below the atomic resonance frequency to suppress 
spontaneous scattering, which would otherwise heat the atoms. For an optical trapping 
laser detuned below the atomic resonance, the induced atomic dipole moment is in-phase 
with the trapping laser field, so that the atoms are attracted to the high field region 
at the trap focus, i.e., the effective trapping potential is U = —a{E 2 )/2, where the 
polarizability a > and (E 2 ) is proportional to the trap laser intensity, time-averaged 
over a few optical cycles. The effective potential then has the same spatial profile as 
the intensity of the focused trap laser. For ultracold atoms, the energy per particle is 
typically quite small compared to the depth of the optical trap, so that the atoms vibrate 
in a nearly harmonic confining potential. The vibration frequencies of the atoms in each 
direction Ui, % = x, y, z of the trap are readily determined by parametric resonance: the 
trap laser intensity is modulated and the size of the cloud is measured as a function of 
modulation frequency. When the modulation frequency is twice the harmonic oscillation 
frequency, the energy of the atoms increases, causing the density profile to increase in 
size. This method permits precise characterization of the trap parameters. 
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All information about the cloud is obtained by absorption imaging: a spatially 
uniform short (several /is) low intensity laser pulse is transmitted through the atom 
cloud, which partially absorbs the light. The shadow cast by the absorption is imaged 
onto a CCD (charge coupled device) array to record the image, which is analyzed 
to extract the column density, integrated along the line of sight. This "laser flash 
photography" method provides real space images with a resolution of a few microns, on 
a time scale short compared to the time scale over which the atoms move significantly 
compared to the spatial resolution. Both non-destructive and destructive imaging 
techniques are possible. In the latter case the entire cloud is destroyed by the laser pulse 
in order to make the most complete possible image. One then runs the same experiment 
multiple times, with nearly identical initial conditions, to obtain an average picture of 
temporal evolution. While the CCD measures just the density or correlations, in 
fact it is possible to extract density- density or correlations from the noise on an 
image [72, 73]. 

What do experimental measurements actually look like? In Fig. 4 are shown 
absorption images from the 2002 Duke experiment on elliptic flow [46]. In the 
experiments, N = 7.5 x 10 4 atoms in each of the two lowest hyperfine states were 
cooled to degeneracy in a C0 2 laser trap, with a reduced temperature T/Tpj between 
0.08 and 0.2. Here, Tpj = hio^GN) 1 ^ 3 jks is the Fermi temperature for an ideal gas at 
the center of a harmonic trap,^ where Co = 2n x 2160(65) Hz is the geometric mean of 
the trap oscillation frequencies, measured by parametric resonance as described above. 
For these parameters, the Fermi temperature is Tpi = 7.9 /iK. For an ideal gas in the 
trap, the Fermi radii are o~ x = 3.6 /xm in the narrow directions and o~ z = 103 fim in the 
long direction. 

2.2. Scale Invariance and Universality 

Studies of trapped ultracold Fermi gases have provided important information about the 
phase diagram, the equation of state, transport properties, and quasiparticle properties 
of strongly correlated Fermi gases. This is possible because under the conditions 
typically encountered in the experiments local properties of the trapped gas directly 
correspond to equilibrium properties of the homogeneous Fermi gas. Consider the 
ground state of N harmonically trapped fermions. Hohenberg and Kohn showed that 
the solution of the iV-body Schrodinger equation is equivalent to the minimum of the 
energy functional [74] 

E[n(x)} = J dx{£{n{x))+n{x)U(x)) , (5) 

ft Tpj is defined by Tpj = Tp(no(0)), where Tp(n) = h 2 kp(n) 2 /{2mkB) is the local Fermi temperature 
of a non-interacting gas evaluated at the center of the trap. An equivalent definition is that 
^bTfj = Epj, where Epi = hio^N) 1 ^ 3 is the Fermi energy of N non-interacting fermions in a harmonic 
trap. The advantage of Tpj is that it only depends on N and u>. 



Strongly Correlated Quantum Fluids 



16 



where n(x) is the density, subject to the condition J dxn(x) = N, S(n) is the energy 
density functional, and U (x) is the trap potential. If the density is sufficiently smooth 
we can write S(n) as a function of the local density and its gradients. On dimensional 
grounds we have 

£(n(x)) = ^ n(x)^ + ^! (Z^Wl! + o (V 4 n(x)) , (6) 

where Co,c\,... are dimensionless constants. At unitarity the coefficients c, are pure 
numbers, but for a finite scattering length they become functions of na 3 . To first 
approximation we can neglect the gradient terms. Then the density is given by 
n(x) = n eq (fi — U(x)), where n eq (/x) is the equilibrium density at the chemical potential fi 
and zero temperature. This approximation is known as the local density approximation. 
Gradient terms are suppressed by (tu±/ fi) 2 ~ l/(A z iV) 2 / 3 , where \ z = uj z /lo_l is the trap 
deformation. Typical experiments involve \ z ~ 0.025 to 0.1 and N > 10 5 , so corrections 
beyond the local density approximation are quite small. These arguments generalize to 
systems at non-zero temperature. In this case the density of the trapped system is 
n(x) = n eq O - t/(x),T). 

The equilibrium density can be determined from the equation of state, P = P(n, T), 
through the thermodynamic relation''' n = dP/dfi. In the following we also frequently 
refer to the relation P = P(n, T) as the equation of state. At unitarity the interaction 
is scale invariant and the only scales in the many-body system are the inter-particle 
distance n _1//3 and the de Broglie wave length, given in Eq. (1). Dimensional analysis 
implies that the equation of state must be of the form 

fc2 5/3 

P(n,T) = ^—f(n\l B ), (7) 
m 

where f(x) is a universal function. At zero temperature the pressure is proportional to 
n 5 / 3 jm. This implies, in particular, that the pressure is given by a numerical constant 
times the pressure of a free Fermi gas. The same is true for the energy per particle and 
the chemical potential. It has become standard to denote the ratio of the energy per 
particle of the unitary gas and the free Fermi gas as the Bertsch parameter £, 



o 



Bertsch posed the calculation of the parameter £ as a challenge problem to the many- 
body physics community in 1999 [75]. At the time, the problem was stated in the 
context of a toy model for dilute neutron matter; see Sec. 3.1. 

Using thermodynamic identities we can show that Eq. (7) implies P = |e, where e is 
the energy density. This relation is analogous to the equation of state of a scale-invariant 
relativistic gas, P = |e, as discussed in Sec. 3.2. For a trapped gas the relation between 

' Here and in the remainder of this review we have dropped the subscript "eq." 
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pressure and energy density implies a Virial theorem: in a harmonic trap, the internal 
energy of the system is equal to the potential energy of the trapping potential [76, 77], 

These universal relations have been extended in many ways; see [78] for a review. 
An important class of relations, discovered by Tan, connects the derivative of 
thermodynamic quantities with respect to 1/a to short range correlation in the gas. 
Tan defined the contact density C via [79, 80, 81] 

de h 2 „ . . 

■C, (10) 



d(a x ) Airm 

where the derivative is taken at constant entropy density. The contact density appears 
in a number of thermodynamic relations. The universal equation of state, for example, 
is given by 

2 h 2 

P=-e+—^C. (11) 

3 127rma v ; 

More remarkable is the fact that C controls short distance correlations in the dilute 
Fermi gas. The tail of the momentum distribution is given by 

n„{k)^tf (|a|- 1 <A;<r - 1 ) , (12) 

where C = J d 3 xC(x) is the integrated contact, n a (k) is the momentum distribution* 
in the spin state a, and r is the range of the interaction. There are similar expression 
for the asymptotic behavior of other correlation observables like the static and dynamic 
structure factors, and the dynamic shear viscosity; see [82] for a review. In this Focus 
Issue, Kuhnle et al. present a comprehensive set of measurements of the contact as 
a function of interaction strength and reduced temperature [83]. These results can be 
compared to new theoretical predictions discussed by Hu et al. [71]. 

Below the critical temperature for superfluidity the superfluid flow velocity v s can 
be viewed as an additional thermodynamic variable. The response of the pressure to 
the superfluid velocity defines the superfluid mass density 



o 2 p 

Ps = mn s = --— 



(13) 



where the derivative is taken in the rest frame of the normal fluid, meaning v n = 0. The 
superfluid mass density can be measured using rotating clouds [84]. The second moment 
of the trap integrated value of the superfluid mass density determines the quenching of 
the moment of inertia. New measurements of the moment of inertia can be found in [85] . 

For small values of n|a| 3 the equation of state P(n, T) can be computed in 
perturbation theory. This program was initiated by Lee, Yang, and Huang [86, 87]. 

* The momentum distribution is normalized as J d)i/(2TT) 3 n cr (k) = N a , where N a is the total number 
of atoms in state a. 
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At unitarity weak coupling methods can be used in the limit of high temperature. This 
is based on the observation that the binary cross section at unitarity is a = 4n/k 2 . At 
high temperature the mean momentum is large and the average thermal cross section 
is small. The equation of state can be written as an expansion in n\ dB , which is the 
well-known Virial expansion. We have 

P = nk B T{l + b 2 (n\ 3 dB ) + 0((nX 3 dB ) 2 )} , (14) 

with b 2 = — 1/(2 v^2) at unitarity [88, 89]. Analytic approaches in the non-perturbative 
regime n\ dB ~ 1 are based on extrapolating to the unitary limit from different regimes in 
the phase diagram. For this purpose the phase diagram has been studied as a function 
of the strength of the interaction, the number of species, and the number of spatial 
dimensions. The oldest theory of this type is Nozieres-Schmitt-Rink (NSR) theory 
[90, 91, 92], which is based on a set of many-body diagrams that correctly describe 
both the BCS and BEC limits. NSR theory works surprisingly well, despite the formal 
lack of a small parameter at unitarity. For example, the basic form of the critical 
temperature sketched in Fig. 3 is correctly reproduced. Modern theories of this type 
are typically based on self-consistent T-matrix approximations; see [93, 94]. Another 
idea is to generalize the unitary Fermi gas to 2N spin states [95]. Mean field theory 
is reliable in the limit iV — > oo, and the interesting case N = 1 can be studied by 
expanding in 1/N. This method is of interest in connection with holographic dualities, 
because the gravitational dual is expected to be classical in the limit that the number 
of degrees of freedom is large. Finally, it was proposed to use the number of dimensions 
as a control parameter. The unitary limit is perturbative in both D = 2 and D = 4 
spatial dimensions [96]. The interesting case D = 3 can be studied as an expansion 
around D = 2 + eorD = 4 — e dimensions [97]. 

These methods are promising, but currently the only techniques that provide 
reliable and systematically improvable results in the regime n\ dB ~ 1 are quantum 
Monte Carlo calculations. At zero temperature the standard technique is Green function 
Monte Carlo (GFMC) [98, 99]. This method relies on a variational initial wave function, 
which is used as the initial condition for an imaginary time diffusion process. The 
Monte Carlo method suffers from a fermion sign problem which is addressed using the 
fixed node approximation. At finite temperature a number of groups have performed 
imaginary time path integral Monte Carlo (PIMC) calculations [100, 101, 102, 103]. 
These calculations do not rely on variational input, and they do not suffer from a sign 
problem, but they are formulated on a space-time lattice and require an extrapolation 
to zero lattice spacing. A new technique is bold diagrammatic Monte Carlo (BDMC), 
which is based on sampling the sum of all Feynman diagrams. The method suffers 
from a sign problem, but convergence in the regime above the critical temperature for 
superfluidity was found to be very good [104]. 
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2. 3. Experimental Determination of the Equation of State 

The equation of state describes a functional relation between key thermodynamic 
variables, such as the pressure T) as a function of chemical potential \i and 
temperature T. In Fermi g stated in Sec. 2.1, what we can actually measure 

is the density profile of a trapped cloud. There are various techniques for translating 
density measurements into thermodynamic quantities, all relying on the local density 
approximation. 

The first thermodynamic measurements took place at Duke in 2005, where global 
thermodynamic variables were measured. In the initial experiments [60], a controlled 
amount of energy was added to the trapped cloud by abruptly releasing it from 
the optical trap, allowing it to expand hydrodynamically by a known factor, and 
then recapturing the cloud after a selected expansion time, thereby increasing the 
potential energy. After allowing the gas to equilibrate, the cloud profile was measured 
to determine an empirical temperature, which was later calibrated by comparing to 
theoretical cloud profiles, predicted as a function of reduced temperature. The resulting 
energy versus temperature curve was compared to that measured for an ideal Fermi 
gas in the same trap, and showed a departure from ideal gas behavior at a certain 
temperature, which yielded an estimate of the superfluid to normal fluid transition 
temperature. However, the results suffered from being model-dependent, as calibration 
of the empirical temperature relied on comparison of the measured cloud profiles with 
theoretical predictions. To avoid this model dependence, in 2006 the JILA group 
measured the potential energy of the strongly interacting cloud of 40 K as a function 
of the ideal Fermi gas temperature that was obtained after an adiabatic sweep of the 
bias field to the noninteracting regime above resonance [105]. 

Model-independent determination of thermodynamic quantities was done by the 
Duke group in 2007 [65], where the total energy E and total entropy S of a trapped 
cloud were measured from cloud profiles, by exploiting the universal behavior of a unitary 
Fermi gas at a Feshbach resonance. From Eq. (9) we know that for a harmonic trapping 
potential the total energy is twice the average potential energy, E = 2(U) = 3mu^(z 2 ). 
Hence, by measuring the harmonic oscillator frequency and mean square cloud size, the 
total energy is readily determined from cloud images. This method was demonstrated 
experimentally in [62]. With the energy in the strongly interacting regime measured 
from cloud images, the entropy is determined by adiabatically sweeping the magnetic 
field to a weakly interacting regime. There, the entropy and mean square cloud size are 
readily calculated as a function of the temperature, yielding the entropy as a function of 
the mean square cloud size of the weakly interacting gas. Adiabatic behavior is verified 
by a round-trip sweep. Fig. 6 shows the energy per particle as a function of the entropy 
per particle, in the universal regime. The temperature is determined by fitting a smooth 
curve [66, 2]. Hu, Liu, and Drummond combined the demonstrated universal behavior 
of the global thermodynamic quantities by reanalyzing the measurements in 40 K and 
showing that the 40 K and 6 Li data fit on a single thermodynamic curve [106]. 
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Figure 6. Total energy per particle of a strongly interacting Fermi gas in the universal 
regime versus the entropy per particle. The blue dots show the entropy obtained 
by adiabatically sweeping the magnetic field from 840G to 1200 G, where the gas 
is weakly interacting. The red dots are the theoretical calculations using a second 
Virial coefficient approximation. The green curve is a fit using two power laws, which 
determines the temperature T = dE/dS. From [2]. 



As already mentioned briefly in Sec. 2.1, model- independent measurement of global 
thermodynamic variables was superseded in 2010 by model-independent measurement of 
local thermodynamic quantities, which can be directly compared to predictions, within 
the local density approximation. Equation 4, determines the local pressure P from 
the local density n and local chemical potential, \x = fi g — U, where \x g is the global 
chemical potential and U is the known trapping potential. Absorption imaging directly 
yields the column density n(x,z) = J dyn(x,y,z), for an imaging beam propagating 
along y. In a cylindrically symmetric harmonic trap U = muj_(x 2 + y 2 )/2 + moo 2 z 2 /2, 
where u± and u z are the radial and axial trapping frequencies, respectively, we can 
write —dfi = dll = muj^pdp = mcu 2 ^ dx dy / (2tt) from which we see that the pressure 
is determined from the doubly integrated density, i.e., the integrated column density 
n(z) = j dx dyn(x, y, z). For a 50-50 mixture of spin-up and spin-down fermions with n 
the total density, the pressure is determined as a function of z, P(z, T) = mui\ n(z) / (2tt). 
For x = y = 0, the corresponding chemical potential is /i(z) = \i g — muj 2 z 2 /2, so 
that P(fi, T) is determined if the temperature and global chemical potential can be 
determined. 

To determine the temperature, the Tokyo group [67] used the temperature 
calibration by the Duke group to obtain T from the total energy and hence from the 
mean square cloud size [66]. An improved version of this calibration is described in 
Ref. [2] of this Focus Issue. The ENS group [68] directly determined the temperature 
by using a weakly interacting 7 Li impurity to measure the temperature of a strongly 
interacting 6 Li gas. The global chemical potential was determined from the wings of the 
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Figure 7. Measured energy versus the temperature obtained from the calibration of 
Ref. [2] (red dots); For comparison, we show the data obtained by the ENS group [68] 
(black dots) and the theory of Hu et al. [71] (green curve). From [2]. 



density profile using a fit based on a Virial expansion, yielding P(/a, T), from which all 
other local thermodynamic quantities, such as the energy density and entropy density, 
were determined. These results enable a direct comparison with predictions. Further, 
the total entropy S and total energy E of the trapped cloud are readily determined by 
integration, and can be compared with the corresponding global quantities measured 
by the Duke group. The total energy versus temperature is shown in Fig 7, where the 
improved temperature calibration of Ref. [2] from this Focus Issue is displayed for the 
Duke data. As quite different methods were employed to make the measurements, their 
close agreement indicates the correctness of the data. 

The latest studies of local thermodynamics by the MIT group [36] use refined 
methods, where measurements of the isothermal compressibility k directly from the 
density profiles replaces temperature measurements, yielding an equation of state 
n(P,n). This eliminates the determination of the temperature and local chemical 
potential from fits at the edges of the cloud, which produced the most uncertainty in 
previous work. For this method, the three dimensional density n(x, y, z) is determined 
by tomographic imaging, using an inverse Abel transform [107] to determine n from the 
measured column density. The trap potential is carefully characterized by determining 
the surfaces of constant density, hence constant chemical potential, in a very shallow 
trap where the axial z trap potential is almost perfectly harmonic, and therefore known. 
Equation 4 yields the pressure 

/fi POO 
dfin(fi',T)= dUn(U,T), (15) 
-oo Ju 

where the unknown global potential \x g is not needed, since the integral is over the 
known trap potential [36]. The compressibility is the change in density with respect to 
a change in the local trapping potential U, and is determined from the density profile 
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Figure 8. Experimentally measured thermodynamics for a unitary Fermi gas at MIT. 
Left panel: energy, free energy, chemical potential, as a function of temperature. Right 
panel: entropy, as a function of temperature. See the text for an explanation of the 
curves. Used with permission from Ref. [36]. 
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since dP = ndfi = —n dU at constant T. Thus, the actual observed equation of state 
is the functional relation n(ft, P), measured directly from the density distribution, the 
clear and direct experimental observable for ultracold quantum g discussed in 

Sec. 2.1. From n(P, k), all other local thermodynamic quantities are determined [36], 
as shown in Fig. 8. 

These experiments provide the best current value for the Bertsch parameter of 
£ = 0.376(5), which is consistent with the value obtained in measurements of global 
quantities, £ = 1 + (3 = 0.39(2) [66]. These measurements can also be compared 
to theory predictions. The two most recent quantum Monte Carlo calculations give 
f < 0.383(1) [108] and f = 0.39681°;°°™ [109]. See [109] for an extensive compilation of 
analytic results and earlier Monte Carlo calculations. 

Figure 8 shows several representations of the equation of state, and provides a 
glimpse of the level of precision that can be achieved in present experiments. In the left 
panel are shown the chemical potential /i, energy E, and free energy F. The right panel 
shows the entropy per particle S/^Nks) versus T/T F . The chemical potential (red solid 
circles) is normalized by the Fermi energy; energy (black solid circles) and free energy 
(green solid circles) are normalized by E Q = ^NEp, which is the energy per particle in 
a uniform Fermi gas. At high temperatures all quantities approximately track those for 
a non-interacting Fermi gas, shifted by £ n — 1 with £ n ~ 0.45 (dashed curves). The peak 
in the chemical potential roughly coincides with the onset of superfluidity. In the very 
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Figure 9. Experimental Fermi gas phase diagram. Experimental measurements of 
the Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein condensate (BEC) crossover 
for ultracold Fermi gases (compare to sketch in Fig. 3). Since a magnetic Feshbach 
resonance is used to tune the interaction strength, one axis is magnetic field. In this 
measurement, the BCS side is on the B > side, so Fig. 3 should be reversed for 
comparison; it is not possible to go deep into the BEC side because the lifetimes 
become too short in this experiment. AB < 0.6 contains the strongly interacting 
region kp\a\ > 1, and the dashed lines indicate uncertainty in the precise position of 
the Feshbach resonance. Used with permission from Rcf. [54]. 

low temperature regime, I^/Ep, E/E and F/F all approach £ (blue dashed line). At 
high temperatures, the entropy closely tracks that of a non-interacting Fermi gas (black 
solid curve). The open squares are from the self-consistent T-matrix calculation [94]. A 
few representative error bars are shown, representing mean ± standard deviation. 

2.4- Experimental Studies of the Phase Transition 

The first measurements of the phase transition at unitarity were made at JILA, shown 
in Fig. 9, by a pair projection technique. After the pairs were created in the BCS and 
unitary regimes, a rapid magnetic field sweep was used to pairwise project the fermions 
onto molecules, to protect them during subsequent expansion measurements, where 
the molecular momentum distribution was measured to determine the condensed pair 
fraction. The fraction of near-zero momentum molecules is interpreted as the percentage 
of condensed Fermi unitary or BCS pairs. In this early experiment 40 K was used, and 
the measured temperature was that of an ideal Fermi gas, rather than the temperature 
of the interacting gas. This ideal Fermi gas temperature was obtained by ballistic 
expansion after an adiabatic sweep to the weakly interacting regime above resonance, 
yielding the condensed pair fraction versus ideal gas temperature and magnetic field. 
Despite the uncertainties regarding the temperature and the pair conversion efficiency we 
observe that the shape of the transition line is qualitatively consistent with theoretical 
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Figure 10. Superfluid phase transition of a unitary Fermi gas. Specific heat per 
particle, Cy / {Nkg), as a function of quantum degeneracy, T/Tp. The phase transition 
is clearly evident at T c /Tf = 0.167(13). Adapted with permission from Rcf. [36]. 



expectations, as summarized in Fig. 3. 

Initial estimates of the critical parameters of the trapped gas were done by 
the Duke group, first based on model-dependent measurements of the energy versus 
temperature [60] and later based on model-independent measurements of the total 
energy E and entropy S, obtained as described above. Assuming that a phase-transition 
would be manifested in a change in the scaling of E with S, two power-laws were used to 
fit the E(S) data, one to fit the high temperature data and one to fit the low temperature 
data, which joined at a point S c . The continuity of the energy and temperature were 
used as a constraint, and the critical entropy S c was estimated from the joining point 
that minimized the x 2 for the fit [66]. While the results obtained by this method were 
consistent with predictions, they suffered from a dependence on the form of the fit 
function, making the uncertainty difficult to quantify. Further, the global observables 
suffer from the averaging that masks an abrupt local phase transition, initially near the 
trap center. 

Using the refined local measurements described in Sec. 2.3, the MIT group has 
traced out the phase transition in the unitary regime, without any need for rapid 
magnetic sweeps, fitting parameters, thermometry, or indeed any kind of theory besides 
elementary thermodynamic considerations. The main goal of these experiments was 
to obtain a cusp-like signal of the phase transition to superfluidity in a unitary Fermi 
gas, by focusing in on second order derivatives of the pressure, where a cusp appears. 
Although superfluidity was established by creation of vortex lattices [59], the actual 
phase transition itself had only been indirectly observed. 

In Fig. 10 is shown the specific heat per particle, clearly displaying the superfluid 
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phase transition. Experimentally, the specific heat is derived from the compressibility 
k and the pressure P [36], 

C v 5 Tp { _ 1\ 

^iV = 2^ v ^J' (17) 
where R = k/kq and p = P/Pq are normalized to the non- interacting Fermi gas 
compressibility Ko = f^^: an d pressure P = ^nE F , respectively. Using n = (dP/dfi)T, 
the compressibility can be written k = {l/n 2 ){dn/dfi)T = {l/n 2 ){d 2 P/d 2 ^)T- As k is a 
second derivative of the pressure, the specific heat shows a clear cusp-like signature, 
Fig. 10. Qualitatively, the behavior of Cy can be understood as follows: As one 
approaches the phase transition from above, T/T c > 1, the compressibility increases 
due to the attraction between fermions; below the phase transition, T/T c < 1, the 
compressibility decreases because fermions are bound into pairs, and it becomes more 
difficult to squeeze the gas, i.e., to change the single particle density. 



2.5. Universal Hydrodynamics and Transport 

Transport properties of the unitary Fermi gas are of interest for several reasons. The first 
reason is related to the main theme of this review: holographic dualities suggest a new 
kind of universality in the transport properties of strongly interacting quantum fluids. 
We expect, in particular, that the shear viscosity to entropy density ratio is close to the 
value rj/s ~ h/ {kirks) originally discovered in the QGP, and first obtained theoretically 
using the AdS/CFT correspondence [11], where AdS is a special maximally-symmetric 
spacetime described in detail in Sec. 4.2, and CFT stands for conformal field theory. 
The second reason is that transport properties are very sensitive to the strength of 
the interaction, and the types of quasiparticles present in the system. The Bertsch 
parameter, which characterizes the effect of interactions on the energy per particle, 
varies by about a factor of two between the weak coupling (BCS) and strong coupling 
(unitarity) limits. The shear viscosity, on the other hand, changes by many orders of 
magnitude. Finally, quantum limited transport has also been observed in systems that 
are of great practical significance, in particular in the strange metal phase of high T c 
compounds; see the contribution to this Focus Issue by Guo et al. [110]. 

Transport properties have been studied experimentally by exciting hydrodynamic 
modes, such as collective oscillations [56, 111, 112, 113, 114], collective flow [46, 115], 
sound [116], and rotational modes [117]. In a system that can be described in 
terms of quasiparticles the hydrodynamic description is valid if the Knudsen number 
Kn = l m tp/L, the ratio of the mean free path l mfp to the system size L, is small. E In the 
unitary gas the mean free path is Z m f p = l/{na), where n is the density and a = 4ir/k 2 
is the universal cross section. In the high temperature limit the thermal average cross 

E Criteria for the validity of hydrodynamics can also be formulated if there is no underlying quasiparticle 
description, a situation that is of great interest in connection with holographic dualities. In this case, 
hydrodynamics is based on a gradient expansion of the conserved currents. The ratio of the 0(v) to 
0(dv) terms in the stress tensor is known as the Reynolds number, Re = vLmn/rj. Validity of the 
gradient expansion requires that the Reynolds number is large. 



Strongly Correlated Quantum Fluids 



26 



section is a = 4A^ B . The Knudsen number of a unitary gas confined in a cigar-shaped 
harmonic trap is 



where we have taken L as the radius in the narrow, or z direction. Here, N is the 
number of particles, \ z was defined previously as the aspect ratio of the trap, and Tpi is 
the global Fermi temperature for a harmonically trapped ideal gas; see Sec. 2.1. Using 
N = 2 • 10 5 and A 2 = 0.045 as in [113] we conclude that hydrodynamics is expected 
to be valid for T < 5Tpj. For the Fermi gas viscosity measurements described in this 
Focus Issue [115, 2] the maximum temperature is T ~ 1.5 Tpj and Kn < 0.09. 

Nearly ideal hydrodynamic behavior was first observed in the expansion of a unitary 
Fermi gas after release from a deformed trap [46]; see Fig. 4. For a ballistic gas 
the expansion reflects the isotropic local momentum distribution in the trap. As a 
result the gas expands in all directions and the cloud slowly becomes spherical. For a 
hydrodynamic system the expansion is driven by gradients in the pressure. In the case 
of a deformed cloud the gradients are largest in the short direction of the trap, and the 
expansion takes place mostly in the transverse direction. As a result the cloud eventually 
becomes elongated along what was originally the short direction. This phenomenon is 
analogous to the elliptic flow observed in heavy ion collisions, as described in Sec. 3.5. 
What is also remarkable is the fact that even though the gas becomes more dilute as 
it expands this effect is compensated by the growth in the mean cross section. As 
a consequence, the gas remains hydrodynamic throughout the expansion. Ballistic 
behavior sets in eventually only because of imperfections, such as the fact that the 
scattering length is not truly infinite. 

The role of dissipative effects, in particular shear viscosity, was first studied in 
collective modes. The radial breathing mode can be excited by removing the confining 
potential, letting the gas expand for a short period of time, and then restoring the 
potential. Hydrodynamic behavior can be established by measuring the frequency of 
the breathing mode. For an ideal fluid lu = a/10/3 ujq whereas in a ballistic system 
oj = 2ujq [118, 119]. The transition from ballistic behavior in the BCS limit to 
hydrodynamics in the unitary limit was observed experimentally in [111, 112]. In the 
hydrodynamic regime damping is expected to be dominated by dissipative terms in the 
equations of fluid dynamics. The energy dissipation is given by 



where V{ is the fluid velocity, rj is the shear viscosity, ( is the bulk viscosity, k is 
the thermal conductivity, and all derivatives, divergences, and gradients are spatial. 
At unitarity the system is scale invariant and the bulk viscosity is expected to 
vanish [120, 78]. This prediction was experimentally checked by Cao et al. [2, 121]. 






(19) 
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Thermal conductivity is not important because the system remains isothermal in ideal 
hydrodynamics. Temperature gradients only appear due to shear viscosity, and their 
contribution to dissipation is higher order in the gradient expansion. This means that 
damping is dominated by shear viscosity. 

Collective modes in ideal fluid dynamics are described by scaling solutions of Euler's 
equation. This means that shape of the density profile does not change during the 
evolution, and that the velocity is linear in the coordinates. The solution is analogous 
to Hubble flows in cosmology, and the Bjorken expansion of a QGP, as discussed in 
Sec. 3.5. In the case of a scaling solution the shear stresses d{Vj are spatially constant, 
and the energy dissipated only depends on the spatial integral of r\. On dimensional 
grounds we can write 77 = hna n . For a scale invariant system a n is only a function of 
the dimensionless variable n 2//3 ^, 2 /(mfcsT), i.e., a function of the reduced temperature 
T/Tp. The reduced temperature varies across the trap, but for a given fluid element 
it remains approximately constant during the hydrodynamic evolution of the system. 
This implies that the damping constant of a collective mode is related to the spatial 
average of the shear viscosity in the initial equilibrium state, 



The extracted values of (a n ) can be converted to the trap averaged shear viscosity to 
entropy density ratio by using the measured entropy per particle. This type of analysis 
was originally carried out in [122, 123], where it was observed that rj/s < 0.5h/kB in 
the vicinity of T c . More recently, Cao et al. showed that in the high temperature 
regime a n exhibits the scaling behavior expected from the solution of the linearized 
Boltzmann equation [115]. The shear viscosity due to elastic 2-body scattering has the 
form 77 ~ npl m { p , where p ~ Ajg is the mean quasiparticle momentum. As we saw 
above, Z m f p ~ l/(na) ~ l/(n\^ B ). This implies that at high temperature 77 ~ A^. The 
coefficient of proportionality was determined by Bruun et al. [8, 124]. They find 



There is an important problem related to Eq. (21) that affects the extraction of the 
shear viscosity from experiments with scaling flows. Equation (21), which is reliable in 
high temperature or low density part of the cloud, is independent of the density". As 
a consequence the integral in Eq. (20) diverges in the low density region. A solution 
to this problem was proposed in [126, 127]: In the low density regime the viscous 
relaxation time ~ 77/(n/ceT), which is the time is takes for the dissipative stresses to 
relax to the Navier Stokes form r](ViVj + V jVi — |5ijVfcWfc), becomes very large. Since 
the dissipative stresses are initially zero, taking relaxation into account suppresses the 

II This is a general property of the viscosity of dilute gases, and was first noticed by Maxwell. The 
result was experimentally confirmed by Maxwell himself, who measured the damping of oscillating discs 
in a partially evacuated container [125] 
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contribution from the dilute corona. A simplified version of this approach was used in 
Cao et al. [115, 2]. The relaxation time and its relation to the spectral function of the 
shear tensor is discussed in the contribution by Braby et al. [128]. 

The shear viscosity drops with temperature and is expected to reach a minimum 
near T c . In this regime quantum effects are important. T-matrix calculations can 
be found in [129] and in the contributions by LeClair and Guo et al. in this Focus 
Issue [130, 110]. It is also possible that in this regime quasiparticle descriptions break 
down completely, and the most efficient description of the unitary gas near T c is in 
terms of a suitable weakly coupled holographic dual. Progress towards constructing 
holographic duals of non-relativistic quantum fluids is summarized in Sec. 4.3.4. 

In the low temperature superffuid regime the appropriate description is superfluid 
(two-fluid) hydrodynamics [131]. Superfluid hydrodynamics predicts the existence of 
additional hydrodynamic modes, in particular second sound, and contains additional 
transport coefficients that come into play if there is relative motion between the 
superfluid and normal components of the fluid. There are proposals for exciting second 
sound modes in the literature [132], but these ideas have not been confirmed yet. At very 
low temperature the shear viscosity is expected to be dominated by phonons, similar 
to liquid helium or dilute Bose gases. Elastic phonon scattering gives r] ~ 1/T 5 [133], 
whereas inelastic processes can give a slower increase at low temperature, rj ~ 1/T [9]. 
These predictions are difficult to verify experimentally because the phonon free path is 
quite large. 

The short mean free path in the strongly interacting normal fluid suggests that not 
only the viscosity, but also other transport coefficients may exhibit universal behavior. 
The spin diffusion constant was recently studied by Sommer et al. [134, 135]. In the 
first paper Sommer et al. observed collisions between polarized Fermi gas clouds. The 
colliding clouds are initially very far from equilibrium, but at late times the system 
relaxes diffusively. The corresponding relaxation time can be used to measure the spin 
drag and the spin diffusion constant D s . The spin diffusion constant in the homogeneous 
system is defined by Fick's law, 



where j s is the spin current, and M = — nj, is the polarization. In this Focus Issue 
Sommer et al. follow up on these studies by measuring the damping of the spin dipole 
mode in strongly polarized gases [135]. A theoretical study of the spin drag relaxation 
rate for a repulsive gas is presented by Duine et al. in this Focus Issue [136]. They show 
that spin fluctuations enhance the spin drag in the vicinity of the Stoner ferromagnetic 
transition. Similarities in the spin transport in the unitary gas and graphene are studied 
in this Focus Issue by Miiller and Nguyen [137]. 



A calculation of the diffusion constant based on the two-body Boltzmann equation 
can be found in [138, 134]. They find 



Js = -D S VM, 



(22) 




(23) 
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Similar to the shear viscosity, the spin diffusion constant drops with decreasing 
temperature. Near the critical temperature D s is expected to approach the universal 
value D s ~ H/m. This behavior is indeed observed in the experiment [134]; see also 
the recent analysis of Bruun and Pethick [139]. It is interesting to compare this result 
to the observed minimum of the shear viscosity. Shear viscosity governs the rate of 
momentum diffusion. The associated diffusion constant is D v = r)/(mn). Near T c we 
have t]/s ~ 0.5h/kB and s/n ~ ks- This implies D rj ~ 0.5h/m, comparable to what is 
seen in spin diffusion. A similar correlation between shear viscosity and diffusion was 
observed in the QGP, as discussed in Sec. 3.7. 

2.6. The BCS-BEC Crossover in Lattices 

A question quite distinct from that of the BCS-BEC crossover in ultracold Fermi gases 
in the continuum is how fermions pair, and reach the unitary regime, in a discrete 
context. Optical lattices, or crystals of light, are created from interfering laser beams. 
They make a sinusoidal standing wave that traps fermions via the AC Stark effect. The 
strongly discretized regime, i.e., in the lowest band(s) and tight binding approximations, 
is quite distinct from that explored by both the QGP and continuum ultracold quantum 
gas experiments. Essential features of weakly interacting fermions in this context have 
already been explored and understood [140]; however, the strongly interacting regime 
remains in question. Moreover, the lattice regime presents an additional challenge for 
holographic duality, as it too has a unitary point where l/(kpa) — > 0. The question 
of how to correctly model this problem in the context of a lattice is subtle, as simply 
attaching a band index to fermionic fields leads to hundreds of bands deep in the BEC 
regime, and is therefore numerically and practically intractable. 

In the continuum, a single channel model treating only fermions qualitatively 
reproduces the phase diagram of Figs. 3 and 9. This qualitative approach can be 
taken for instance by using the generalized BCS ansatz, technically valid only for 
weak interaction and high density, at the mean-field level for arbitrary interaction 
and finding the chemical potential self-consistently by fixing the average number of 
particles [141]. This method was first used by Eagles in the context of superconductivity 
in low carrier concentration systems [142] and later by Leggett [143] and Nozieres and 
Schmitt-Rink [90] explicitly for the BCS-BEC crossover at zero temperature and finite 
temperature, respectively, together called NSR theory. In NSR one begins with a lattice, 
but then quickly takes a continuum or low momentum limit. Interestingly, if this limit 
is not taken one obtains a completely incorrect prediction: instead of the BEC critical 
temperature tending to a constant non-zero value as l/(kpa) — > oo, as in Fig. 3, the 
critical temperature tends to zero algebraically. This unphysical tendency is because 
in order to hop or tunnel between lattice sites, a pair of fermions must be broken and 
tunnel one by one. Then for very strong pair binding energy the process is prohibitively 
expensive. In fact, even on the BCS side such models utilizing a single band have been 
shown to fail quantitatively for quite small values of kpa < [144, 145]. This model is 
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called the Hubbard Hamiltonian [146] and is also a proposed model for high-temperature 
superconductivity [147]. 

A solution to this quandry is to introduce a two-channel model incorporating 
both fermions and bosons [148, 149, 150, 151]. The fermions then represent the 
unbound atoms scattering at threshold (open channel) while the bosons represent 
the weakly bound molecule brought into resonance (closed channel). Such two- 
channel models are mentioned briefly in the condensed matter context as far back 
as 1985 [152] and explored seriously in partial form starting in 1995 [153]; they 
continue to be explored as the Cooperon model, in current research on high temperature 
superconductivity [154]. However, only in ultracold quantum gases have all hopping, 
interaction, and interconversion terms been included to make a Fermi-Bose Hubbard 
Hamiltonian in the lattice [155]. A series of papers attempted this approach with 
steady improvement over time, starting from simply attaching a band index to 
each channel [155, 156] up to performing a renormalization procedure to produce 
effective molecules in the lattice and minimize the number of bands that must be 
included [157, 158, 159, 160, 161]. The latter method has led to a Hamiltonian that is 
so far from the original Hubbard Hamiltonian that it has been given a new name, the 
Fermi Resonance Hamiltonian. 

The Fermi Resonance Hamiltonian is 

H*=-t f E E a L%-+^ E X>*- EE^4a»(24) 

+ E Va E^« + E E^-i.«-* [4<At«*U + h - c -] , 

aGM i a&M ijk 

where a\ a creates a particle with spin a in the lowest open channel band Wannier state 

centered at lattice site i; d\ a creates a particle in the a th dressed molecule Wannier 

"(f) 

state centered at site i; n\J is the number operator for fermions in the lowest Bloch 
band; and hf^ is the number operator for the a th dressed molecule state. The set of 
dressed molecules A4 which are included dynamically can be determined on energetic 
and symmetry grounds from the two-particle solution. In order, the terms in Eq. (24) 
represent tunneling of atoms in the lowest Bloch band between neighboring lattice sites 
i and j; the energy Eq = ^qEi^/N 3 of a fermion in the lowest band with respect to 
the zero of energy; tunneling of the dressed molecular center of mass between two lattice 
sites i and j, not necessarily nearest neighbors; detunings of the dressed molecules from 
the lowest band two-particle scattering continuum; and resonant coupling between the 
lowest band fermions at sites j and k in different internal states and a dressed molecule 
at site %. The Fermi Resonance Hamiltonian is a two-channel resonance model, between 
unpaired fermions in the lowest band, and dressed molecules nearby in energy [161]. 
Among other unusual predictions it makes, in contrast to usual Hubbard physics, is 
significant diagonal hopping, pairing between atoms which do not lie along a principal 
axis of the lattice to form a dressed molecule, and multiple molecular bound states 
induced by the lattice. Thus the lattice problem is indeed quite different both from the 
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continuum crossover problem and well-known Hubbard physics. 

In practice, the solution of the crossover problem in the lattice must occur in 
three steps. First, for a given lattice strength and interacting strength the two-body 
problem must be solved exactly, including a renormalization procedure to get effective 
Wannier states for the molecules. Then the coefficients in Eq. (24) are calculated. Then 
the Hamiltonian itself must be solved. Potential solution methods range from matrix 
product state methods in ID to mean field, quantum Monte Carlo, and dynamical 
mean field theory methods in higher dimensions. Holographic duality may offer new 
approaches to this Hamiltonian in the future. The crossover problem on the lattice 
remains very much an open problem and therefore we return to it in Sec. 5. 

2. 7. Recent and New Directions in Crossover Physics 

In this subsection, we cover some of the new directions in unitary Fermi gases and 
related systems not discussed in previous sections, as explored in this Focus Issue. 

2.7.1. New Experimental Probes We require new experimental probes into unitary 
gases in order to better understand their behavior. Thus nearly all papers on ultracold 
fermions in this Focus Issue suggest an experimentally measurable effect. However, 
three papers realize new probes directly in experiments. First, Riedl et al. [85] measure 
the quenching of the moment of inertia in the unitary regime. This quenching is a 
well-known signal of superfluidity, as starting with liquid 4 He it has been demonstrated 
that the moment of inertia drops below its classical value as one decreases a system 
through its critical temperature for the superfluid transition. Unlike other methods of 
measurement in ultracold Fermi gases such as frequency and damping rate of collective 
excitations in response to small perturbations, quenching of the moment of inertia 
zeroes in specifically on superfluid effects. Thus one can distinguish between nearly 
ideal fluid dynamics in the normal phase, which is chracterized by very low viscosity 
but allows rotational flow, and irrotational superfluid hydrodynamics. A second new 
experimental probe is spin diffusion measurement, covered by three sets of contributors, 
one experimental [135] and two theoretical [136, 138], as we already touched on in 
Sec. 2.5. The spin diffusion constant D s and the momentum diffusion constant D v 
exhibit a similar temperature dependence, and an analogous relation between heavy 
quark and momentum diffusion is expected in the QGP. Overall, diffusion constants 
offer another point of comparison for bounds predicted by holographic duality. Finally, 
a third new experimental probe is Tan's contact density, given in its original form in 
Eq. (10), and covered in this Focus Issue in both an experimental [83] and a theoretical 
contribution [71]. Tan's contact density connects microscopic scattering properties to 
macroscopic thermodynamic observables throughout the BCS-BEC crossover. Finally, 
as shown in Fig. 7, an improved temperature calibration is realized in this Focus Issue 
to better characterize the key ratio of the viscosity to the entropy [2]. 
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2.7.2. Solitons to Polarons We turn now to new directions in theory. The simplest 
and oldest approach to describing the BCS-BEC crossover is BCS mean field theory, 
in the form of the Boguliubov-de Gennes (BdG) equations describing quasiparticles in 
the BCS phase. Despite the fact that these equations are technically only valid for 
weak interactions —kpa <C 1, and are based on a quasiparticle picture, one can attempt 
an extrapolation through the crossover by self-consistently solving the gap equation, 
BdG equations, and density and normalization of quasiparticles. Both Spuntarelli et 
al. [162] and Baksmaty et al. [163] describe this approach in detail in this Focus Issue; 
the method is numerically challenging as it requires solution of simultaneous nonlinear 
partial differential equations. Spuntarelli predicts features of dark solitons through the 
crossover. Solitons appear as robust local minima in the gap with a strong indicator in 
the density as well, even at unitarity. Solitons are also described by holographic duality 
in this Focus Issue by Keranen et al. [164], providing a point of comparison to mean 
field theory. Further recent work on solitons through the BCS-BEC crossover can be 
found in Refs. [165, 166, 167]. Baksmaty et al. use the BdG approach to shed light on 
experimental data in imbalanced Fermi gases, i.e., those with more spin up than spin 
down fermions. Imbalanced Fermi gases in elongated traps are shown to display strong 
violation of the local density approximation in the form of phase separation, a quite 
different perspective from the thermodynamic one. 

Imbalanced Fermi gases, sometimes also called spin-polarized Fermi gases, open up 
new territory for observing quantum phase transitions and have been used recently to 
explore a Fermi liquid description of the strongly interacting normal phase [168] and 
to observe a Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) state in one-dimensional Fermi 
gases [169]. The FFLO state is a pairing between two different Fermi surfaces, resulting 
in a finite center-of-mass momentum of the pairs. Spin-dependent scattering is suggested 
as a novel experimental probe in a theoretical contribution by Sheehy [170], to measure 
the pairing spin gap in locally imbalanced Fermi gases, where the imbalance is induced by 
coupling to a spin-dependent potential. Imbalanced gases are also used in spin diffusion 
measurements, as described in Sec 2.5. In the limit in which a Fermi gas is strongly 
imbalanced, a polaron model can be used, where the minority species is dressed by the 
Fermi sea of the majority species and the resulting quasiparticle is called a polaron; in 
the extreme case a single atom of the minority species is considered, leading to a single 
polaron. Exactly how such concepts apply at unitarity where quasiparticle concepts 
tend to breakdown is an open question. In this Focus Issue Sadeghzadeh et al. [171] 
predict the existence of a new metastable state consisting of a Fermi sea of polarons. 

2.7.3. Disorder and Quantum Phase Transitions As discussed in Sec. 2.6, lattice 
physics is only barely beginning to be treated in Fermi gases through the BCS-BEC 
crossover. Two new directions are treated for ultracold Fermions in optical lattices in 
this Focus Issue. First, Han and Sa de Melo [172] treat the BCS-BEC crossover in the 
presence of disorder. They find the superfluid near unitary to be much more robust 
against disorder than in the BCS and BEC regimes, and give some useful observations 
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on how the different practical realizations of disorder in optical lattices can lead to 
different considerations for crossover physics. Second, lattices give one access to strongly 
correlated systems which have nothing to do with the BCS-BEC crossover. That is, near 
quantum phase transitions one generally obtains scale-invariant phenomena that can 
potentially be studied with holographic duality. Zhang et al. [173] cover experiments 
on one such system in this Focus Issue, ultracold bosons near the Mott insulator to 
superfluid transition in the Hubbard model in two dimensions. Then the same kinds 
of issues as we have discussed throughout Sec. 2 reappear, from a universal equation of 
state to transport. The first half of this paper forms a very readable overview for the 
general reader, while the second half gives a status report on experiments. 

The massive ground state degeneracy induced by lattices is the principle underlying 
the achievement of strongly correlated states in bosons. One can obtain similar 
degeneracy either by rotating a system to achieve such degeneracy in the lowest Landau 
level and get quantum Hall physics, or by taking advantage of the spin degree of freedom. 
The latter strategy is pursued by Shlyapnikov and Tsvelik [174]. Chromium, a Bose- 
condensed atom with spin 5 = 3, provides a practical working example, and they find a 
commensurate-incommensurate quantum phase transition, among other possibilities in 
a rich phase diagram. 

2.7.4- Trimers: From Efimov States to SU(3) Although we very briefly mentioned 
three-body losses in Sec. 2.1, in fact three-body losses are the dominant loss mechanism 
at unitarity, whether for fermions or for bosons. Three body physics shows some real 
surprises in ultracold fermions: even when there is no two body bound state at all, 
one finds a denumerably infinite and universal^ set of three-body bound states in the 
strongly interacting regime, called Efimov states. Efimov states can occur for various 
combinations of bosons and fermions of same or different masses. The evidence for 
such series of states is found in three-body losses. Braaten and Hammer wrote a review 
back in 2006 [175]. However, this area continues to develop rapidly, motivated strongly 
by ultracold experiments. In this Focus Issue, Wang and Esry cover this important 
area [176], extending our knowledge of Efimov physics beyond the low-energy threshold 
collisions typical of ultracold gases to the higher energy regime relevant to colliding 
clouds. They treat both positive and negative scattering length and both broad and 
narrow Feshbach resonances [32]: in narrow Feshbach resonances Efimov effects are 
less universal in nature, requiring a second scattering potential parameter beyond the 
scattering length, namely, the effective range. Losses in fermions also show a puzzling 
feature, a maximum at a magnetic field below the pole of the Feshbach resonance where 
the scattering length diverges. Zhang and Ho suggest an explanation for this anomaly 
in terms of an interplay between atom-atom, atom-dimer, and dimer-dimer interactions. 
Based on these considerations, and incorporating temperature and trapping parameters, 
they develop rate equations which match all experimental loss features, covering data 

^ By universal we mean that the details of the scattering potential don't matter for Efimov states. 
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in four different labs taken over seven years. 

Finally, so far we have discussed solely SU(2) physics. However, the hyperfine spin 
manifold of ultracold fermions allows for a variety of spin structures to be created. In 
particular, the preparation of degenerate three-state mixtures with three-fold (SU(3)) 
symmetric attractive interactions opens up new territory in the connection between 
cold atoms and nuclear matter. Such a three-state gas can be used to explore 
pairwise superfluidity and spontaneous symmetry breaking, color superconductivity, and 
a superfluid to trimer gas transition that mimics the deconfmement to hadronization 
transition in quark matter [177]. In symmetric three-component gases, it is possible 
that more than one pairing field can simultaneously become nonzero [178]. Magnetism 
and domain formation also can be studied [179]. A key issue is to perform experiments 
at low density in order to reduce the 3-body decay rates oc n 3 , which are not Pauli- 
suppressed in three-state systems. Three-body decay arises from recombination, which 
is enhanced and modified near a Feshbach resonance, where Efimov states exist with 
binding energies well beyond the threshold region [176]. By tuning away from resonance 
and using radio-frequency dressing to precisely symmetrize the pair-wise interactions, 
and confining the atoms in a lattice to increase the critical reduced temperature, the 
goal of creating cold-atom analogs of quark matter appears to be within reach [177]. 

3. Quantum Chromodynamics, the Quark-Gluon Plasma, and Heavy Ion 
Collisions 

In this section we will give an introduction to the physics of the quark-gluon plasma 
(QGP). This is a vast subject and we will not attempt to provide a comprehensive 
review here. Recent reviews of various issues related to quark matter and the QGP can 
be found in [180, 181, 182, 183, 20, 184] and the results from the experimental program 
at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Laboratory are 
summarized in [185]. Some results from the Large Hadron Collider (LHC) program 
at the European Organization for Nuclear Research (CERN) can be found in the 
contributions [186, 187, 189] and in the recent review [190]. 

We will focus on observables and experimental results that are important in 
establishing the presence of strong correlations in the QGP, and that provide direct 
connections to the physics of ultracold atomic gases and holographic dualities. In Sec. 3.1 
we introduce quantum chromodynamics (QCD), the theory of the strong interaction, 
and discuss the QCD phase diagram. The differences between weakly coupled QCD 
plasmas and strongly coupled QCD fluids are discussed in Sec. 3.2. The experimental 
study of high temperature QCD matter is based on heavy ion collisions. We provide 
an overview of the experimental program in Sec. 3.3 and discuss specific observables 
(multiplicities, flow, energy loss and heavy quarks) in Sec. 3.4-3.7. Finally, a discussion 
of open issues can be found in Sec. 5.2. 
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3.1. Quantum Chromodynamics and the Phase Diagram 

QCD describes the behavior of strongly interacting matter, such as protons, neutrons, 
and nuclei, as well as the hot and dense matter created in heavy ion collisions, in terms of 
quarks and gluons and their interactions. The complicated phenomenology of the strong 
interaction is encoded in a deceptively simple Lagrangian. The Lagrangian is formulated 
in terms of quark fields q c a j and gluon fields A a . Here, a = 1, . . . , 4 is a Dirac spinor index 
(corresponding, in the Dirac representation, to quarks and anti-quarks with spin up and 
down), c = 1, . . . , N c is a color index, and / = up, down, strange, charm, bottom, top is a 
flavor index. In our world the number of colors is N c = 3, but as a theoretical laboratory 
it is useful to consider theories with N c ^ 3. The dynamics of the theory is governed 
by the color degrees of freedom. The gluon field A ^ is a vector field (like the photon) 
labeled by an adjoint color index a = 1, . . . , 8. The octet of gluon fields can be used to 
construct a matrix valued field A^ = A a ^, where A a is a set of traceless, Hermitian, 
3x3 matrices. Repeated indices are assumed to be summed throughout our treatment. 
The QCD Lagrangian is 

N f 

C = --Gl v Gl v + Qfi^D, - m f )q f , (25) 
/ 

where G® u is the QCD field strength tensor 

G% = d,Al - d u Al + gr hc A\Al , (26) 

and f abc = 4i Tr([A a , A b ]A c ) is a set of numbers called the SU(3) structure constants. 
The covariant derivative acting on the quark fields is 

iD,q=(^d, + gAl^q, (27) 

and raj is the mass of the quarks. The different terms in Eq. (25) describe the interaction 
between quarks and gluons, as well as nonlinear three and four-gluon interactions. It is 
important to note that, except for the number of flavors and their masses, the structure 
of the QCD Lagrangian is completely fixed by a local SU(3) color symmetry. 

For many purposes we can consider the light flavors (up, down, and strange) to 
be approximately massless, and the heavy flavors (charm, bottom, top) to be infinitely 
massive. In this limit the QCD Lagrangian contains a single dimensionless parameter, 
the coupling constant g. If quantum effects are taken into account the coupling becomes 
scale dependent. At leading order the running coupling constant is 

9/ 9N 16tt 2 11 2 



where q is a characteristic momentum and Nf is the number of active flavors (Nf = 3 
if we consider charm, bottom and top to be infinitely heavy). The result in Eq. (28) 
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Figure 11. Schematic phase diagram of QCD as a function of temperature T and 
baryon chemical potential /i. QGP refers to the QGP, and sQGP/wQGP denote the 
regions of the phase diagram where the plasma is strongly /weakly coupled. The 
CFL (color-flavor locked) phase is the color superconducting phase that occurs at 
asymptotically large chemical potential. LHC/RHIC denote the regions of the phase 
diagram that are being explored by the experimental heavy ion programs at the LHC 
and RHIC. The red and black points denote the critical endpoints of the chiral and 
nuclear liquid-gas phase transitions, respectively. 



implies that, as a quantum theory, QCD is not characterized by a dimensionless coupling 
but by a dimensionful scale, the QCD scale parameter Aqcd- This effect is called 
dimensional transmutation [191] . We also observe that the coupling decreases with 
increasing momentum. This is the phenomenon of asymptotic freedom [192, 193]. The 
flip side of asymptotic freedom is anti-screening, or confinement: the effective interaction 
between quarks increases with distance. 

In massless QCD the scale parameter is an arbitrary parameter (a QCD "standard 
kilogram"), but once QCD is embedded into the standard model and quarks acquire 
masses by electroweak symmetry breaking the QCD scale is fixed by the choice of units 
in the standard model. The numerical value of Aqcd depends on the renormalization 
scheme used to derive Eq. (28). Physical masses, as well as the value of b , are 
independent of this choice. In the modified minimal subtraction (MS) one finds 
Aqcd ^ 200 MeV; see the QCD section in [194]. 

Asymptotic freedom and the symmetries of QCD determine the main phases of 
strongly interacting matter that appear in the QCD phase diagram shown in Fig. 11. In 
this figure we show the phases of QCD as a function of the temperature T and the baryon 
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chemical potential /x. The chemical potential [i controls the baryon density p, defined 
as 1/3 times the number density of quarks minus the number density of anti-quarks. 

At zero temperature and chemical potential the interaction between quarks is 
dominated by large distances and the effective coupling is large. As a consequence, 
quarks and gluons are permanently confined in color singlet hadrons, with masses of 
order Aq C d- For example, the proton has a mass of m p = 935 MeV. + . 

If we think of the proton as composed of three constituent quarks this implies that 
quarks have effective masses tuq ~ mp/3 ~ Aqcd- This should be compared to the 
bare up and down quark masses which are of the order 10 MeV. 

Strong interactions between quarks and anti-quarks lead to a vacuum condensate of 
qq pairs, (qq) ~ — Aq CD [195, 196, 197]. This vacuum condensate spontaneously breaks 
the approximate chiral SU(3)l x SU(3)r flavor symmetry of the QCD Lagrangian. 
Chiral symmetry breaking implies the existence of Goldstone bosons, particles with 
masses smaller than Aqcd- These particles are known as pions, kaons, and etas.* 

The quark-anti-quark condensate (qq) is analogous to the di-fermion condensate 
(qq) that characterizes BCS superconductors [34]. BCS pairing involves particles with 
opposite momenta near the Fermi surface. The chiral condensate is an order parameter 
for pairing of fermions and anti-fermions near the surface of the Dirac sea. The most 
important difference between BCS pairing and chiral condensation is that because of 
the finite density of states near the Fermi surface BCS pairing can take place at weak 
coupling. Chiral condensation, on the other hand, only happens at strong coupling." 
We will see below that in QCD at large baryon density there is a transition between 
qq and qq pairing. This transition is driven by the competition between the stronger 
coupling in the qq channel and the growing density of states for qq pairing. 

At very high temperature quarks and gluons have thermal momenta p ~ T ^> Aqcd 
(Note that we are writing momenta in units of energy, again setting c = 1). Asymptotic 
freedom implies that these particles are weakly interacting, and that they form a plasma 
of mobile color charges, the QGP [199, 200]. We note that the argument that the QGP 
at asymptotically high temperature is weakly coupled is somewhat more subtle than 
it might appear at first sight. If two particles in the plasma interact via large angle 
scattering then the momentum transfer is large, and the effective coupling is weak 
because of asymptotic freedom. However, the color Coulomb interaction is dominated 
by small angle scattering, and it is not immediately clear why the effective interaction 
that governs small angle scattering is weak. The important point is that in a high 

+ Note that we are quoting the mass in units of energy, setting the speed of light equal to unity 
* The SU(3)l x SU(3)r is explicitly broken by quark masses and the mass of the charged pion is 
m v = 139 MeV, which is not much smaller than Aqcd- The lightest non-Goldstone particle is the rho 
meson, which has a mass m p = 770 MeV. 

3 In ultracold atomic gases in the BEC regime pair condensation involves prc-formed pairs. Whether 
or not pre-formed qq pairs exist in in the QGP above T c is still being investigated; see for example 
[198]. We should note that there is no completely rigorous criterion for the existence of a bound state 
embedded in a plasma. In practice, researchers have looked for peaks in the spectral function associated 
with the correlator (qq(x)qq(O)) . 
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temperature plasma there is a large thermal population [n ~ T 3 ) of mobile charges that 
screen the interaction at distances beyond the Debye length ~ l/(gT). We also note 
that even in the limit T ^> Aqcd the QGP contains a non-perturbative sector of static 
magnetic color fields [201]. This sector is strongly coupled, but it does not contribute 
to thermodynamic or transport properties of the plasma in the limit T — > oo. 

The plasma phase exhibits neither color confinement nor chiral symmetry breaking. 
This means that the high temperature QGP phase must be separated from the low 
temperature hadronic phase by a phase transition. The nature of this transition is very 
sensitive to the values of the quark masses. In QCD with massless u, d and infinitely 
massive s,c,b,t quarks the transition is second order [202]. In the case of massless (or 
sufficiently light) u, d, s quarks the transition is first order. Lattice simulations show 
that for realistic quark masses, m u ~ ~ 10 MeV and m s ~ 120 MeV, the phase 
transition is a rapid crossover [203]. The transition temperature, defined in terms of the 
chiral susceptibility, is T c ~ 151 ± 3 ± 3 MeV [204, 205]. 

The transition is believed to strengthen as a function of chemical potential, so that 
there is a critical fi at which the crossover turns into a first order phase transition [206] . 
This point is the critical endpoint of the chiral phase transition. Due to the fermion sign 
problem it is very difficult to locate the critical endpoint using simulations on the lattice. 
A number of exploratory calculations have been performed [207, 208, 209, 210], but at 
this point it is not even clear whether the idea that the transition strengthens as the 
baryon chemical potential increases is correct [211]. The critical endpoint is interesting 
because it is the only point on the phase transition line at which the correlation length 
diverges (there is a similar endpoint on the nuclear liquid-gas transition line). This 
means that the critical point may manifest itself in heavy ion collisions in terms of 
enhanced fluctuations [212]. The idea is that one can tune the baryon chemical potential 
by changing the collision energy, since lower beam energy allows for more stopping of 
the baryons in the initial state and therefore leads to higher chemical potential. As 
the collision energy is varied one looks for non-monotonic behavior of fluctuation and 
correlation observables. A typical observable is the variance (divided by the mean) 
of the net number of protons (protons minus anti-protons) in a finite sub-volume, 
((ANp) 2 ) I (N p ) with AN P = N p — (N p ). A beam energy scan is now under way at 
RHIC, and similar scans have been performed at lower energy as part of the CERN 
fixed target program [213, 214]. 

At low temperature the first phase one encounters as the chemical potential is 
increased from zero is nuclear matter, a strongly correlated superfluid composed of 
approximately non-relativistic neutrons and protons. It is interesting to note that 
nuclear matter in weak n -H- p + e + v e equilibrium is neutron rich,^ and that the 

Isolated neutrons are unstable with regard to the decay into protons, electrons, and anti-neutrinos. If 
we ignore the role of electrons then dense matter is composed of equal numbers of protons and neutrons 
because this configuration optimizes the total Fermi energy. Taking electrons into account we observe 
that for a given density of protons and electrons the Fermi energy of the electrons is much bigger than 
that of protons, and the lowest energy state is neutron rich. 
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neutron-neutron scattering length* a nn ~ 20 fm is much larger than the average inter- 
particle spacing, r nn ~ 1.9 fm at nuclear matter saturation density. As a consequence, 
dilute nuclear matter is closely related to the unitary atomic Fermi gases discussed in 
Sec. 2. 

At very large chemical potential we can use arguments similar to those in the 
high temperature limit to establish that quarks and gluons are weakly coupled. The 
main difference between cold quark matter and the hot QGP is that because of the 
large density of states near the quark Fermi surface even weak interactions can cause 
qualitative changes in the ground state of dense matter. In particular, attractive 
interactions between quark pairs lead to color superconductivity and the formation of a 
(qq) condensate. Since quarks carry color, flavor, and spin labels, many superconducting 
phases are possible. The most symmetric of these, known as the color-flavor locked 
(CFL) phase, is predicted to exist at very high density [215, 216]. In the CFL phase the 
diquark order parameter is (qaflpg) ~ e a/3t ABC £f g c- This order parameter has a number 
of interesting properties. It breaks the U(l) symmetry associated with baryon number, 
leading to superfluidity, and it breaks the chiral SU(3)l x SU(3)r symmetry. Except 
for Goldstone modes the spectrum is fully gapped; fermions acquire a BCS-pairing gap, 
and gauge fields are screened by the Meissner effect. This implies that the CFL phase, 
even though it arises from a superdense liquid of quarks, shares many properties of 
superfluid nuclear matter. 

The CFL phase involves equal pair-condensates (ud) = (us) = (ds) of all three 
light quark flavors. As the density is lowered effects of the non-zero strange quark 
mass become more important, and less symmetric phases are likely to appear [180]. 
Calculations based on weak coupling suggest that the first non-CFL phase is a CFL-like 
phase with a Bose condensate of kaons, and the second phase involves a standing meson 
wave superimposed on the kaon condensate [217, 218]. Other possibilities which may 
appear in strong coupling include a phase with up-down pairing only (2SC), or separate 
spin-one condensates of up, down and strange quarks. 

Some guidance for analyzing the competition between these phases may come 
from studying analog models of the quark-hadron phase transition based on ultracold 
atomic gases. A number of authors have pointed out to the possibility of a BCS- 
BEC crossover in dense quark matter [219]. The basic idea is that as the density is 
lowered, quark Cooper pairs become more strongly bound and form a diquark Bose 
condensate. At lower densities one may also encounter a mixture of condensed diquarks 
and unpaired quarks. The physics of this system can be studied using ultracold boson- 
fermion mixtures [24, 25, 220, 221]. Finally, it may be possible to study the transition 
from nucleons, bound states of three quarks, to color superconducting diquarks using 
three-species systems of ultracold fermions [222, 177, 178, 179]. 

Color superconductivity affects both the equation of state and transport properties 
of dense quark matter. These effects may manifest themselves in the structure and 

t 1 fm (Fermi)=10 -15 m is the typical unit of length in nuclear and particle physics. 1 fm is 
approximately equal to the radius of the proton and neutron. 
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evolution of neutron stars. Possible signatures of pairing appear in the mass-radius 
relation, the cooling curve, and the relation between the spin-down rate and the magnetic 
field. There are two contributions in this Focus Issue that touch on these issues. 
Shovkovy and Wang study the bulk viscosity of normal quark matter in the large 
amplitude regime [223]. Anglani, Mannarelli and Ruggieri investigate the interaction of 
collective modes in the CFL phase [224]. This calculation is a first step towards more 
accurate calculations of transport properties of the CFL phase [225, 226]. 

3.2. Weakly vs. Strongly Coupled Plasmas 

At asymptotically high temperature the coupling is weak and properties of the QGP 
can be systematically computed. One of the most basic properties is the equation of 
state. In the following we will consider a related quantity, the entropy density. The 
perturbative expansion of the entropy density is [227, 228] 

s = T 3 {c + c 2 g 2 + c 3 g 3 + ...} , (29) 

where g is the QCD coupling constant from Eq. (28) and we have chosen units such 
that h = ks — 1. The expansion is performed with the coupling g = g{jx) evaluated at 
a fixed scale ft. Higher order terms contain logarithms of the form log(2irT / p) , which 
combine with lower order terms to eliminate the dependence on the arbitrary scale fi, 
and lead to an expansion in terms of the running coupling g(q) evaluated at a scale 
q ~ 2tiT. The first term in Eq. (29) corresponds to the Stefan-Boltzmann law, and Co 
is proportional to the number of degrees of freedom, 

Co = ^(2(iV c 2 -l)+4iV c iV / ^ , (30) 

where 2(N 2 — 1) is the number of degrees of freedom from the gluons and 4N c Nf from 
the quarks. The naive perturbative expansion is an expansion in powers of g 2 . Odd 
powers of g appear because of infrared divergences. As we approach the phase transition 
the coupling becomes large and higher order terms are no longer small. Indeed, the 
convergence properties of the expansion in Eq. (29) are extremely poor: the series 
shows no signs of converging unless the coupling is taken to be much smaller than one, 
5 < 1, corresponding to completely unrealistic temperatures on the order of 1 TeV. 
The convergence can be improved significantly by using a self-consistent quasiparticle 
expansion [229]. This means that the perturbative expansion is formulated not in terms 
of free quarks and gluons, but in terms of quasi-quarks and quasi-gluons which have 
effective masses and effective interactions. 

Quasi-particle expansions fit lattice data quantitatively down to temperatures 
T ~ 2T C . In this regime s/sq ~ 0.85, where Sq is the entropy density of a non-interacting 
plasma. In the past this was frequently taken as evidence that quasiparticles are not 
strongly coupled at temperatures relevant to the early stages of heavy ion collisions 
at RHIC or the LHC A new perspective on this questions is provided by holographic 
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duality; see Sec. 4. Holographic duality maps strongly interacting quantum field theories 
onto one higher dimensional classical gravity, in particular anti-de Sitter space. In the 
original version of the mapping the quantum field theory is a QCD-like theory known 
as N = 4 supersymmetric (SUSY) Yang-Mills theory [230]. This theory is conformal. 
That means its coupling constant does not run, g 2 (q 2 ) = const, and there is no analog of 
Aqcd- It also does not have a phase transition, and so the theory is in the plasma phase 
for all values of the coupling. The particle content of N = 4 SUSY Yang-Mills theory 
differs from that of QCD. The theory has no quarks, and in addition to gluons it contains 
supersymmetric fermionic partners of gluons called gluinos as well as additional colored 
scalar fields. Nevertheless, the plasma phase shares many features of the QGP, and 
by focusing on a ratio like s/sq we can remove the difference in the number of degrees 
of freedom. Later versions of holographic duality relax conformal and supersymmetric 
requirements, but still do not map precisely onto QCD. 

Using holographic duality it was found that in the limit of strong coupling 
and a large number of colors the entropy density of SUSY Yang-Mills theory is 
s/sq = 0.75 [231]. This value is remarkably close to the result observed in lattice 
QCD calculations near the phase transition, casting doubt on an interpretation of the 
data in terms of weakly coupled quasiparticles. More generally, holographic duality 
demonstrates that thermodynamic properties of the plasma need not be very sensitive 
to the strength of the interaction. 

Holographic duality also shows that, in contrast to equilibrium properties, transport 
properties of the plasma are very sensitive probes of the strength of the interaction. In 
perturbative QCD the shear viscosity of three flavor QCD is [232, 6, 233] 

kT 3 

v = —, — T~*l v ( 31 ) 

g 4 log{/j*/m D ) 

where k = 106.67, /i* = 2.96T and vtld ~ gT is the screening mass. Shear viscosity 
is related to the rate of momentum diffusion. A simple estimate of the shear viscosity 
is [234] 

1 i 1 P foo\ 

V-o n P l ^P-o > ( 32 ) 

O 6 0~t 

where n is the density, p is the mean momentum, Z m f p is the mean free path, and 
is the transport cross section. In a perturbative QCD plasma the typical momentum 
is p ~ T, and the transport cross section is o~t ~ g 4 log(g)T~ 2 . The time scale for 
momentum diffusion is r)/(sT) ~ l/(g 4 log(g)T). We note that in the weak coupling 
limit this number is parametrically large. In a QGP at T = 200 MeV, just above the 
phase transition, we have T _1 ~ 1 fm/c and rj/s ~ 9.2/ g 4 . A typical value of the 
coupling is g ~ 2 (corresponding to a s ~ 0.3), which implies t]/(sT) ~ 0.6 fm/c. 

The strong coupling limit of rj/s in the SUSY Yang-Mills plasma was studied by 
Policastro, Son, and Starinets [235]. They find 

^ = ^(1 + (A-/*)), (33) 
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where A = g 2 N c is the 't Hooft coupling. This result implies that at strong coupling 
momentum equilibration is almost an order of magnitude more rapid than it is in weak 
coupling. The ratio rj/s is also close to a bound that was argued to arise from the 
uncertainty relation [234]. The idea is that, because of quantum mechanics, the product 
plmtp in Eq. (32) cannot become smaller than Planck's constant H (note that we are 
using units in which K = 1). This implies rj/s < §(7) ~ 0-09, where we have used the 
entropy per particle of a free gas [234]. The uncertainty argument has never been made 
precise, since the mean free path estimate is not valid at strong coupling. However, 
the holographic duality calculation was discovered to be quite general. It was shown 
that the strong coupling limit of rj/s is universal in a large class of theories that have 
gravitational duals, and that the 0(A~ 3 / 2 ) corrections are positive [236, 237]. These 
observations led Kovtun, Son and Starinets to make the conjecture that 

n - > ^ (34) 

is a universal bound [11] that applies to all fluids. The status of this conjecture 
is discussed in more detail in Sec. 4, and sketched in Fig. 2. There are some 
known counterexamples involving theories with gravitational duals described by higher 
derivative gravity. However, it seems clear that Eq. (34) applies to a large class of 
theories that include generalizations of QCD. 

Experiments on collective flow in heavy ion collisions at RHIC and the LHC 
(discussed in Sec. 3.4) indicate that the viscosity of the QGP near the critical 
temperature is indeed close to the proposed bound, and that equilibration must be 
very rapid. Together with the observation of a large energy loss of highly energetic 
probes of the plasma these results led to the conclusion that the QGP produced at 
RHIC must be strongly coupled [238, 185]. 

In the case of the SUSY Yang-Mills plasma one can show that strong coupling 
implies the absence of well defined quasiparticles. This can be seen most clearly by 
studying the spectral function piuS) of the stress tensor correlation function. Kubo's 
formula relates the value of the spectral function at zero energy to the shear viscosity. 
The spectral function at finite energy carries information about the physical excitations 
that contribute to momentum relaxation. In weak coupling the spectral function has 
a peak at zero energy. The width of the peak, r ~ g 4 log(g)T, is related to the 
quasiparticle life time. The spectral function in the strong coupling limit was computed 
in [239, 240]. It was found that the spectral density is completely featureless: the 
intercept at zero energy smoothly connects to the continuum contribution p(u) ~ u 4 . 

In the case of the QGP the evidence regarding the existence of quasiparticles 
is ambiguous. The stress tensor correlation function in lattice QCD was studied 
by H. Meyer [4]. The reconstructed spectral function is smooth, but the resolution 
was not sufficient to exclude the presence of quasiparticles. There are a number of 
results that have been interpreted as favoring the existence of quasiparticles. One 
is the fact that lattice calculations of fluctuations of conserved charges, like baryon 
number and strangeness, are compatible with the behavior of a free gas, even near 
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Figure 12. Schematic layout of the ALICE detector at the Large Hadron Collider 
(LUC). The beam lines are the thin yellow tubes entering from both sides. The 
outermost (red) layer is a magnet. The magnetic field helps to discriminate between 
particles of different charge-to-mass ratios. The central gray barrel is a time projection 
chamber and the outer layers inside the magnet contain calorimeters, transition 
radiation and time-of-flight detectors. The innermost layers shown in the blow-up 
contain silicon detectors for precise determination of the initial production vertices. 

T c [241]. An experimental observation that has been cited as evidence for the presence 
of quasiparticles is the approximate quark-number scaling of the elliptic flow parameter 
v-i [242]. Another experimental observable that may shed some light on the quasiparticle 
structure is the heavy quark diffusion constant D, which we will discuss in Sec. 3.7. The 
main observation is that kinetic theory and holographic duality make very different 
predictions about the relation between the momentum relaxation time r)/(sT) and the 
heavy quark relaxation time rriQD/T [181]. 

3.3. Nuclear Collisions: Initial Conditions 

The experimental study of heavy ion collisions at relativistic energies, which were long 
seen as the best environments for the production of the QGP, began in the late 1980's 
with fixed target programs both at CERN near Geneva, and Brookhaven National 
Laboratory (BNL), outside of New York City; see [243, 244, 245] for historical overviews. 
These two programs were essential in building up a vital experimental community, 
formed of larger and larger collaborations as the experiments became more and more 
sophisticated. The programs at CERN and BNL demonstrated that heavy ion collisions 
produce strongly interacting, approximately equilibrated matter. The next step was 
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taken with the turn-on of the Relativistic Heavy Ion Collider (RHIC) at Brookhaven in 
2000 which brought heavy ion physics into the collider era, with accessibility to processes 
calculable with perturbative QCD. Colliders are accelerators with counter-circulating 
beams in which the entire collision energy is available for particle production. RHIC 
collides gold ions with a beam energy such that every nucleon in each of the beams 
has an energy of 100 GeV. Since the rest mass energy of a nucleon is about 1 GeV, 
this implies a relativistic 7 factor of about 100. When two bunches of counter rotating 
beams overlap in each of the RHIC experimental halls, at a rate of thousands of times per 
second, the one billion ions in each bunch typically induce less than one collision. In a 
fraction of these collisions, particularly the ones where the nuclei are head-on, thousands 
of particles are produced by the conversion of kinetic energy into mass energy. 

The particles are recorded in detectors which were originally located in four 
experimental areas located around the ring. The early RHIC program had two large 
multipurpose detectors, STAR and PHENIX, with roughly 500 collaborators each, and 
two smaller detectors, BRAHMS and PHOBOS, with about 50 collaborators each [185]. 
STAR and PHENIX are both large spectrometers. STAR focuses primarily on charged 
hadrons using a time projection chamber of 4m diameter that creates a 3D image of 
collision events at a relatively low rate. PHENIX is a combination of drift chambers, 
particle identification counters and calorimeters that make precise measurements of 
both charged particles and photons. Both of these detectors were inspired by the collider 
detectors of the previous generation and make measurements mainly at large angles with 
respect to the beam directions. To make sure physics was not missed in other regions 
of phase space, BRAHMS and PHOBOS were both designed to make more limited 
measurements near the beam axis, BRAHMS with a narrow-band spectrometer and 
PHOBOS with a large-coverage single-layer silicon detector, sensitive to roughly 75% of 
the total number of charged particles produced per event. A range of the data from these 
experiments will be shown in later sections, and in several of the contributions to this 
Focus Issue [186, 187, 189]. As of 2011, RHIC has completed its eleventh experimental 
run, after colliding ions with a wide range of energies (from 7.7 GeV to 200 GeV) and 
nuclear species (protons, deuterons, copper and gold). 

As RHIC was taking data, the Large Hadron Collider (LHC) at CERN was being 
built, along with its three large detector systems ALICE, ATLAS and CMS. The LHC 
provides heavy ion collisions of primarily lead, slightly larger than gold, at center of 
mass energies up to 28 times that reached by RHIC. This translates into higher particle 
densities and temperatures and, even more significantly, into much higher rates for large 
momentum-transfer processes, such as jets, photons, and heavy flavor, discussed below. 
The ALICE detector was designed around another large time projection chamber, similar 
to but larger than the one in STAR at RHIC, with a set of additional detectors to identify 
different hadron species as well as photons and electrons in limited angular regions. 
The layout of the ALICE detector is shown in Fig. 12. The two larger experiments, 
ATLAS and CMS, were general purpose detectors deploying charged particle tracking 
and hermetic calorimetry over a large solid angle. However, the stringent design 
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Figure 13. Geometry of a high energy heavy ion collision. Left: collision of two 
Lorentz contracted gold nuclei. The beam direction is the z-axis. Right: same collision 
in the transverse plane. The impact parameter is along the x-axis, and the remaining 
transverse direction is the y-axis. 



requirements to search for new high mass particles led to detectors that were quite 
capable for the higher multiplicities expected in heavy ion collisions at the LHC, and 
which have excellent capabilities for high energy processes, similar to that needed for 
Higgs and searches beyond the Standard Model. As of late 2011, the LHC has completed 
its second lead ion run, as well as a first feasibility study for future proton lead collisions. 

In the following we will concentrate on results from Au+Au at the top RHIC energy 
of 100 GeV per nucleon. The transverse radius of a Au nucleus is approximately 6 fm, 
and the duration of a heavy ion event is r ~ (6 — 10) fm/c. The estimate for the lifetime 
comes from hydrodynamic simulations which we will describe in Sec. 3.5. The simplest 
observable in a heavy ion experiment, typically published very soon after a new machine 
becomes operational, is the total multiplicity of produced particles. In Au+Au collisions 
at 100 GeV per nucleon the total multiplicity is about 7000; see Sec. 3.4. Somewhat 
more detailed information is provided by the spectra dN/d 3 p of produced particles 
(e.g. Ref. [246, 247]). The momenta can be decomposed into a transverse momentum 
Pt = Px + Py an d a longitudinal momentum p z ; see Fig. 13. In the relativistic regime 
the natural variable to describe the motion in the z direction is the rapidity, 



At RHIC the energy of the colliding nuclei is 100+100 GeV per nucleon, and the 
separation in rapidity is Ay = 10.6. 

A simple picture of the initial state of the fireball created in the collision was 
suggested by Bjorken [249]. He proposed that the two highly Lorentz contracted nuclei 
pass through each other and create a longitudinally expanding fireball in which particles 




(35) 
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Figure 14. Glauber Monte Carlo calculation showing the collision of two gold nuclei at 
an impact parameter b = 6 fm head-on (left panel) and from a side view (right panel) . 
The Glauber model is a geometric model for the high energy scattering of composite 
quantum mechanical states. Spectator nucleons from the two nuclei are shown in pale 
red and blue, participants are shown as dark red and blue circles. Figure from [248]. 



are produced. In the original model the number of produced particles is independent 
of rapidity, and the subsequent evolution is invariant under boosts along the z axis. 
A simple model for the initial energy density in the transverse plane is the Glauber 
model [250, 248]. The Glauber model is based on the observation that high energy 
scattering can be described in the eikonal (geometric optics) approximation. The initial 
entropy density in the transverse plane is 



s(x ± , b) oc T A (x.± + b/2) 1 - exp (-a NN T a (x_l - b/2)) 

+ T A (x ± - b/2) [l - exp (-a NN T A (x ± + b/2)) 
where b is the impact parameter, 



dz p A (x) 



(36) 



(37) 



is the thickness function, x = (x,y,z) and x_|_ = (x,y), and o~nn(\^s) is the nucleon- 
nucleon cross section. We also define the nuclear density as p^(x). The idea behind 
the Glauber model is that the initial entropy density is proportional to the number of 
nucleons per unit area which experience an inelastic collision, which we call the number 
of participants, N part . Other variants exist. For instance, one can distribute the energy 
density according to the number of binary nucleon-nucleon collisions, iV co ii; see Ref. [251] 
for a comparison. A more sophisticated theory of the initial energy density is provided 
by the color glass condensate (CGC) [252, 253]. This model leads to somewhat steeper 
initial transverse energy density distributions. 
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Figure 15. Left: number of charged particles per unit pscudorapidity rj as as a 
function of 7/ for Au + Au collisions at several different energies. Pscudorapidity 
is defined like rapidity, but using the approximation E ~ \p\, and is more easily 
measured experimentally than the rapidity y. For light particles, such as pions, rj ~ y; 
note, however, that corrections to this relation are biggest at n = 0. Data from the 
PHOBOS collaboration at RHIC [185]. Right: multiplicity dN/drj per participant pair 
at midrapidity for both nucleon-nucleon and nucleus-nucleus collisions as a function 
of the collision energy Data and compilation of earlier results from the ALICE 
collaboration at the LHC; sec [254] for original references. 



It should also be noted that the simple Glauber model is typically not used by 
modern calculations, since it does not account for the large fluctuations which are 
currently attributed to the event-wise variations in the nucleon configuration coming 
from each nucleus and the collision process itself. Figure 14 illustrates this by means of 
a single event simulated by a Glauber Monte Carlo code, which counts participants and 
collisions by means of a simple prescription, that collisions occur when two nucleons are 
within d < \J &nn J ^ •> where ct^n is the total nucleon-nucleon inelastic cross section [248]. 

3.4- Particle Multiplicities 

Particle multiplicities are important both as first measurements at new machines, but 
also as a means of testing the current theoretical approaches in their ability to predict 
the degrees of freedom available at higher energies. 

From the theoretical perspective, multiplicities are argued to be a good proxy 
for the initial gluon density. This relies on one essential non-trivial argument, that 
subsequent scatterings in the QGP, after the initial phase where hard collisions occur 
and energy is deposited in the interaction region, do not increase the overall entropy, 
and thus do not increase the overall multiplicity. This is certainly not exactly true, as 
both thermalization and viscous effects during the hydrodynamic evolution will produce 
some amount of entropy, but there is evidence that the total entropy is dominated by 
initial particle production [255]. Assuming that entropy is conserved one can derive the 
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Bjorken estimate for the initial entropy density, 



so = 



3.6 (dN 
ttR 2 t \ dy 



) 



(38) 



where To is the thermalization time, R is the nuclear radius and dN/dy is the total 
number of particles per unit rapidity. This estimate is based on a model of the 
hydrodynamic evolution which we will discuss in more detail in Sec. 3.5. The left panel in 
Fig. 15 shows experimental results for dN/dy of charged particles at RHIC. Assuming 
that A?" (all) = 1.5A(charged) (most particles are pions) we get dN/dy\ y= a ~ 975 for 
Au+Au collisions as 100 GeV/nucleon. A conservative estimate for the equilibration 
time is r = 1 fm/c; see Sec. 5.2. This corresponds to an initial entropy density 
So — 30 fm -3 . For a weakly interacting QGP this implies an initial temperature To ~ 230 
MeV, and an initial energy density eo = 5 GeV/fm 3 . This number is significantly larger 
than the critical energy density for forming a QGP, e cr i t — 1 GeV/fm 3 [256]. 

Experiments typically focus on studying how multiplicities scale with beam energy 
and collision geometry. Empirically the multiplicity per unit rapidity scales as a small 
power of the beam energy; see the right panel of Fig. 15. In nucleus-nucleus collisions 
data ranging from the fixed target SPS program all the way to the LHC are described 
by dN/dy\ y= Q ~ s a , where s is the square of the center of mass energy in the nucleon- 
nucleon system and a ~ 0.15. In nucleon-nucleon collisions the exponent is somewhat 
smaller, a ~ 0.11. 

As a first approximation, it is natural to expect that heavy ion collisions can be 
built by a linear superposition of proton-proton or, more accurately, nucleon-nucleon 
collisions. This is the basis of the wounded nucleon model of the 1970's which was 
invoked to explain how multiplicities in proton-nucleus collisions tended to scale as 
^ch A = Nch P xN pa _n/2 [257, 258]. In this expression, A part is the number of participating 
nucleons estimated either using a Glauber model, or by counting the number of slow 
proton tracks, e.g. in emulsion. In this picture, soft particle production was expected 
to result from the excitation and subsequent independent decay of individual nucleons. 
Hard processes, in which single partons (quarks and gluons) from one nucleon scattered 
off partons of another nucleon, were also discovered in the 1970's. As these were not 
simply decay products of a soft excitation but were short-range elastic scatterings, they 
were expected to scale linearly with the number of binary collisions, also predictable 
using a Glauber model. 

In order to compare to experimental results, calculations typically scale the 
produced multiplicity by the number of participants, since the multiplicity is dominated 
by soft particles, with transverse momenta under 2 GeV. The experiments use the 
observed distribution of multiplicities at a fixed beam energy to estimate the centrality 
of each collision, usually expressed as a percentage of the total inelastic cross section. 
The 10% highest multiplicity events are the 0-10% centrality bin, the next highest are 
the 10-20% and so on. These bins are then compared to the events with the 10% smallest 
impact parameter, or 10% highest multiplicity events, in either a simple Glauber model 
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or one which includes some model of the experimental fluctuations [248]. 

The first comparisons of theory and experiment at RHIC using the multiplicity in 
the most central events were surprising. While many models predicted a combination 
of soft and semi-hard physics processes, it was the models that were predominantly 
composed of soft processes that came closest to the experimental data. This included 
particularly models based on parton saturation, which were preferred by the early data. 
The situation is a bit different at the LHC, as described briefly in the contribution by 
Steinberg [186]. At the higher energies, models which include a combination of hard 
and soft processes seem to do better at predicting the multiplicities in the most central 
collisions. However, the centrality dependence, i.e., the relationship between the number 
of measured charged particles and the number of participating nucleons (or effective 
nucleon-nucleon collisions), is essentially the same at RHIC and the LHC [259]. At 
present it is not clear how to reconcile this observation with the apparent increase in 
the fraction of bulk particle production stemming from hard processes. 

3.5. Hydrodynamic Flow 

The particle density measured by the early experiments is not only valuable to compare 
with different production models. It is also an essential empirical input that can be 
used to determine the initial conditions for hydrodynamical models of the subsequent 
evolution. The hydrodynamic description is valid if the system is in local thermal 
equilibrium, and if the thermodynamic variables are varying smoothly. If the system 
has a kinetic description in terms of quasiparticles this implies that the ratio of the 
mean free path over the characteristic length scale of the system, known as the Knudsen 
number, must be small. There is evidence that this condition is satisfied during the early 
stages of the evolution, but it must eventually break down during the late, dilute, stage 
of the evolution, since the mean free path increases with diluteness. The time when the 
mean free path becomes larger than the system size, more precisely, the expansion time 
multiplied by the mean velocity, is called freezeout. 

There are several experimental observations that provide evidence for the 
assumption that heavy ion collisions create a thermalized state. The first observation is 
that the overall abundances of produced particles, including very rare ones, is described 
by a simple thermal model that contains only two parameters, the freezeout temperature 
T and baryon chemical potential [i [183] . The second observations is that for momenta 
less than about 2 GeV the spectra dN/d 3 p of produced particles follow a Boltzmann 
distribution characterized by the freezeout temperature and a collective radial expansion 
velocity [182]. This radial flow is clearly seen from the fact that the spectra of heavy 
particles, which receive a larger momentum boost from the collective velocity field, have 
a larger apparent temperature than the spectra of light particles. 

The third piece of evidence in support of not only thermalization, but early 
thermalization, that means equilibration significantly before freezeout, is the observation 
of strong azimuthal anisotropics, typically called elliptic flow, in non-central heavy ion 
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collisions. Elliptic flow represents the collective response of the system to pressure 
gradients in the initial state. The analogous phenomenon in ultracold atomic gases is 
discussed in Sec. 2.1. At finite impact parameter the initial state has the shape of an 
ellipse, with the short axis along the x-direction, and the long axis along the ^-direction; 
see Fig. 13. This implies that pressure gradients along the a;- axis are larger than along 
the y-axis. Hydrodynamic evolution converts the initial pressure gradients to velocity 
gradients in the final state. Elliptic flow is sensitive to early thermalization because the 
initial anisotropy that drives elliptic flow disappears with time. Elliptic flow is quantified 
in terms of the second Fourier coefficient v 2 of the particle distribution in the transverse 
plane, 

v (pt) (l + 2 Vl (pr) cos(0) + 2v 2 {p T ) cos(20) + ...), (39) 



dN 
d 3 p 



Pz=0 



where <fi is the angle between the momentum vector and the x-axis. Odd Fourier 
moments like v\ and V3 arise mainly from fluctuations in the initial state [260, 261]. 
Techniques for measuring v 2 and systematic trends in the results are described in some 
detail in the contribution by Snellings [187]. Snellings emphasizes that measurements 
of v 2 constrain the equation of state and the transport properties of the QGP. A 
similar conclusion is reached in the contribution of Lisa et al. [188] to this Focus Issue. 
These authors focus on direct measurements of the shape of the final state using HBT 
(Hanbury-Brown-Twiss) interferometry. 

In this section we will explain how elliptic flow can be used to constrain the ratio 
rj/s introduced in Sec. 3.2. In a relativistic fluid the equations of energy and momentum 
conservation can be written as a single equation 

<9 M T^ = 0, (40) 

where T^ v is the energy momentum tensor. In ideal fluid dynamics the form of T^ v is 
completely fixed by Lorentz invariance, 

= (e + P)vPu v + Prf v , (41) 

where w M is the fluid velocity (u 2 = —1) and rf v = diag(— 1, 1, 1, 1) is the metric tensor. 
The hydrodynamic equations have to be supplemented by an equation of state P = P(e). 
The four equations given in Eq. (40) can be split into longitudinal and transverse 
components relative to the fluid velocity. The longitudinal equation is equivalent to 
entropy conservation, 

d lt (su») = 0, (42) 
and the transverse equation is the relativistic Euler equation, 



^=-r4p V M P > (43) 



where D = u^d^ and = <9 M — u^D. We observe that the inertia of a relativistic 
fluid is governed by e + P. In the non-relativistic limit we can have ~ (1 , rT) and 
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e + P ~ p (energy density is dominated by rest mass energy) . We also find D ~ d t + v ■ d 
and V" 1 ~ 9. These approximations lead to the usual Euler equation. We also note 
that V^P = c^V^e, where c s is the speed of sound. This implies that for a given 
initial energy density gradient the resulting acceleration is determined by the speed of 
sound. Ideal fluid dynamics corresponds to the limit that variations in the hydrodynamic 
variables occur on scales much larger than the mean free path. Dissipation arises from 
the leading gradient terms in the energy momentum tensor. In the rest frame of the 
fluid these corrections have the same form as in non-relativistic fluids. We have 

5T i3 = -7] f <9V + &u l - jjcf J <9 • u^j - (5 ij d ■ u , (44) 

where rj and ( are the shear and bulk viscosity. 

The application of hydrodynamics to relativistic heavy ion collisions goes back to 
the work of Landau [262] and Bjorken [249]. Bjorken discussed a simple scaling solution 
of the equations of fluid dynamics that corresponds to the space-time picture discussed 
in Sec. 3.3. This solution provides a natural starting point for more detailed studies 
in the ultra-relativistic domain. In the Bjorken solution the initial entropy density is 
independent of rapidity, and the subsequent evolution is invariant under boosts along 
the z axis. The evolution in proper time is the same for all comoving observers. The 
flow velocity is 

u ll = 1 {l,0,0,v z ) = {t/r,0,0,z/r), (45) 

where 7 = a/1 — v\ is the boost factor and r = \ft 2 — z 1 is the proper time. The 
velocity field in Eq. (45) solves the relativistic Euler equation (43). In particular, there 
is no longitudinal acceleration. The remaining hydrodynamic variables are determined 
by entropy conservation. Eq. (42) gives 

A \ts(t)\ = (46) 

and s(r) = sqTq/t. For an ideal relativistic gas s ~ T 3 and T ~ 1/r 1 / 3 . We saw in 
the previous section that the initial entropy density is constrained by the final state 
multiplicity. Typical parameters at RHIC are To ~ (0.6 — 1.6) fm and To ~ (300 — 425) 
MeV. We note that the initial temperature is significantly larger than the critical 
temperature for the QCD phase transition. The temperature drops as a function of 
t and eventually the system becomes too dilute for the hydrodynamic evolution to 
make sense. At this point, the hydrodynamic description is matched to kinetic theory, 
generally using the formalism described by Cooper and Frye in the early 1970's [263], 
and the spectra of produced particles are computed. 

In order to quantitatively describe the observed particle distributions several 
improvements of the simple Bjorken model are necessary. First, one has to include the 
transverse expansion of the system [265]. Second, one has to include deviations from 
boost invariance in the longitudinal directions. Rapidity distributions at RHIC most 
likely result from physics somewhere in between the Bjorken scenario, which assumes 
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Figure 16. Differential elliptic flow V2(pt) of identified hadrons from minimum bias 
Au+Au collisions at 200 GeV per nucleon compared to predictions from ideal (non- 
dissipative) hydrodynamics, from [264]. 



boost invariance, and the Landau picture, which assumes complete stopping of the initial 
nuclei [266]. One also has to include realistic equations of state and take into account 
the geometry of the initial state. Results for ^(pt) obtained from a calculation in ideal 
hydrodynamics are shown in Fig. 16. This calculation focuses on the central rapidity 
regime and maintains the assumption of boost invariance. The figure shows the result 
for a given centrality class, corresponding to a specific range of impact parameters. 
Hydrodynamic calculations show that the elliptic flow response Vi is approximately 
linear in the spatial anisotropy (y 2 — x 2 ) / (x 2 + y 2 ) of the initial state. We observe that 
ideal hydrodynamics provides an excellent fit to the RHIC data for p T < (1.0 — 1.5) 
GeV, depending on the particle species. This includes the observed hierarchy in the pt 
dependence of t>2 for different species at low px- The mass splitting of ^(pr) reflects 
an approximate transverse energy = + m 2 ) scaling of the particle spectra in 
hydrodynamics. 

Having established a baseline description of the spectra using ideal hydrodynamics 
we can now discuss the role of dissipative effects. We begin with the effect of shear and 
bulk viscosity on the Bjorken solution. The scaling flow given in Eq. (45) is a solution 
of the relativistic Navier-Stokes equation. If the transverse expansion of the system is 
neglected viscosity does not affect the flow profile but it does generates entropy. We 
find 

Ids _ 1/ Iv + C 
s dr t \ sTr 

The applicability of the Navier-Stokes equation requires that the viscous correction is 



-H 1 -^'- (47) 
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small [234], 



- + 3 J-«- A (Tr). (At 
s 4 s 4 



For the Bjorken solution Tr ~ r 2 / 3 grows with time, and this condition is most restrictive 
during the early stages of the evolution. Using tq = 1 fm and T = 300 MeV gives 
7]/s < 0.6. This result implies that hydrodynamics cannot be used in relativistic heavy 
ions collisions unless the QGP is strongly coupled and the shear viscosity is small. 

It is instructive to study the viscous contribution to the stress tensor in more detail. 
Neglecting bulk viscosity the stresses in the central rapidity slice are given by 

T Z z P ^XX Tyy P ~\~ ~ • (49) 

This means that shear viscosity decreases the longitudinal pressure and increases the 
transverse one. In the Bjorken scenario there is no acceleration, but if pressure gradients 
are taken into account shear viscosity will tend to increase radial flow. At finite impact 
parameter shear viscosity reduces the pressure along the x-direction, and increases the 
pressure in the ^-direction. As a consequence there is less acceleration in the x-direction, 
and elliptic flow is suppressed. 

Viscosity modifies the stress tensor, and via the matching to kinetic theory at 
freezeout, this modification changes the distribution functions f p of produced particles. 
In Ref. [267] a simple quadratic ansatz for the leading correction 5f to the distribution 
function was proposed, 

Sf P = ^~/o(l ± MPcPpdW , (50) 

where fo is the Bose-Einstein/Fermi-Dirac distribution and d^ a u^ is the symmetric 
traceless tensor that appears in Eq. (44). This ansatz is a very good approximation 
to the result of a more involved calculation using kinetic theory [6]. The modified 
distribution function leads to a modification of the single particle spectrum. For a 
simple Bjorken expansion and at large px we find 



SjdN) 1 r] fpr_ 
dN ~ 3T f T f b \T~ f 



(51) 



where dN is the number of particles produced in ideal hydrodynamics, S(dN) is the 
dissipative correction, and Tf is the freezeout time. There is an analogous formula for 
the second Fourier moment of the spectrum, related to v-i [267]. We observe that the 
dissipative correction to the spectrum is controlled by the same parameter 7]/{stT) that 
appeared in the entropy equation. We also note that the viscous term grows with px- 
These results are in agreement with experiment: Deviations from ideal hydrodynamics 
grow with px, and they are larger in smaller systems (which freeze out earlier). More 
detailed analyses can be found in the contributions by Snellings and Nagle, Bearden 
and Zajc [187, 189]. A conservative bound for rj/s at RHIC is rj/s < 0.4, but the best 
fits tend to give values that are even smaller rj/s ~ (0.1 — 0.2). 
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Figure 17. Nuclear suppression factor Raa for a variety of hadron species, photons 
and electrons measured in PHENIX. Note that the dominant source of electrons is 
the decay of heavy quarks (bottom and charm). References to data are indicated on 
the figure itself. The definition of Raa is discussed in the text. The value Raa = 1 
implies that for this particular probe a nucleus-nucleus collision behaves like a simple 
superposition of nucleon-nucleon collisions. Figure provided courtesy of Carla Vale, 
BNL. 



3.6. Jet Quenching 

Another way to directly probe the density of gluons present in the initial state is to 
scatter fast partons, formed in the initial collisions, which can be detected in the final 
state as hadronic jets. This idea was proposed as far back as 1982 by Bjorken, placed 
on firmer theoretical footing in the early 1990's, and finally discovered experimentally 
by RHIC experiments in the early part of the 2000's [268]. Figure 17 shows the nuclear 
suppression factor, defined as 

_ 1 dN/dp^ 

" "dN/W ' { } 

i.e. the yield of a particle species in heavy ion collisions, typically measured as a function 
of its transverse momentum, divided by the similar yield in proton-proton collisions, 
divided by the number of binary collisions calculated from a Glauber model. This 
ratio essentially tests the hypothesis that each individual binary collision has an equal 
probability to induce a hard process identical to that found in proton-proton collisions. 
The latter acts as a reference system in which hard processes can usually be calculated 
perturbatively [269]. As can be seen in Fig. 17, only the production of direct photons 
appears to be unmodified, and even then only below 12-13 GeV, above which isospin 




Figure 18. Nuclear suppression factor Raa for a variety of hadron species, photons 
and Z° bosons measured in Pb-Pb collisions at the LHC, from [190]. The plot 
shows results from the CMS and ALICE collaboration compared to various theoretical 
models. Sec [271, 272] for the original data and [190] for references to the theory 
predictions. 

effects are expected [270]. Conversely, light hadrons like n° and r\ are suppressed by 
about a factor of five above 5 GeV and perhaps slowly rise at high p?. Electrons from 
charmed hadron decays are unsuppressed at low px but quickly fall to a level similar to 
that found for light hadrons, as discussed in more detail in the next section. 

Significant effort has been devoted to the development of a perturbative QCD-based 
formalism that describes the energy loss of fast partons moving through a hot, dense 
medium. The main transport parameter that appears in the theory is the transverse 
momentum diffusion constant q, [274] 

q = p J dqlql^, (53) 

which determines the mean transverse momentum kick of the fast parton per unit length 
traveled. Here, dajdq\ is the differential cross section for scattering of the parton off 
the constituents of the dense medium, and p is the density of the medium. Recent work 
has focused on non-perturbative definitions of q that make no reference to the structure 
of the dense medium; see [184] for a review. The energy loss per unit length scales 
as dE/dz ~ qL for short path lengths (L < L c ), and as dE/dz ~ \J~qE for long path 
lengths. The characteristic length scale is L c ~ <J ' E /q. 
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Figure 19. Jet central-to-peripheral ratio Rcp for different centralities in Pb-Pb 
collisions at the LHC, measured by the ATLAS collaboration [273]. The plots show 
the ratio of the jet yield in central (0 — 10)% and semi-peripheral (50 — 60)% events 
relative to the yield in very peripheral collisions as a function of the transverse energy 
E T of the jet. 



Extracting q from experimental data on jet quenching at RHIC has proven to 
be difficult. Based on an analysis of Raa using a Monte-Carlo implementation of the 
theory of perturbative energy loss [274, 275] the PHENIX collaboration has reported q = 
13.2±^GeV 2 /fhi (the 2a errors are q = 13.2±^GeV 2 /fm) [276]. This number is large 
compared to expectations for a perturbative QCD plasma, q ~ (1 — 2)GeV 2 /fm [277]. 
However, there are significant uncertainties associated with different implementations 
of energy loss. Bass et al. consider three different methods and find that equally good 
descriptions for Raa can be obtained for values of q ranging from q = 2.5 GeV 2 /fm to 
q = 10GeV 2 /fm [278]. 

The LHC offers significantly larger jet production cross sections, a much increased 
p T range, and the capability of detailed studies of not only leading particles but identified 
jets and jet shapes. First results from the LHC are shown in Figs. 18 and 19. Figure 18 
shows the observable Raa introduced above. The low px behavior agrees with results 
at RHIC, despite the much larger collision energy. At high p? the suppression is 
smaller, Raa ~ 0.5, but there is no hint of Raa approaching 1. This is seen even 
more dramatically in Fig. 19, which shows the suppression factor Rcp for identified 
jets, not just high px hadrons. The ratio Rcp is defined relative to the yield in very 
peripheral nucleus- nucleus collisions. We observe that R C p is approximately equal to 
0.5 for jets with transverse energies as high as 300 GeV. These results have important 
implications for the dependence of energy loss on the energy density of the medium, but 
the theoretical analysis of the LHC data is still very much in progress. A discussion of 
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Figure 20. (a) Raa of heavy-flavor electrons in 0-10% central collisions [282] 
compared with ir° data [283] and model calculations (curves I [284], II [285], and 
III [286]). Model III is the Langevin calculation discussed in the text for two different 
values of the diffusion constant. Model II contains some extra processes, in particular 
final state coalescence. Model I is based on energy loss. The box to the far right at 
Raa = 1 shows the uncertainty in Taa- (b) v^ F of heavy- flavor electrons in minimum 
bias collisions compared with 7r° data [287] and the same models. 

the first jet data from the LHC can be found in the contribution in this Focus Issue by 
Steinberg [186]. 

There are some important connections between q and other transport parameters. 
In a quasiparticle description the cross section in Eq. (53) is the same cross section that 
governs momentum diffusion of approximately thermal particles, i.e., shear viscosity. 
This implies that rj/s oc T 3 /q, where the constant of proportionality is roughly 1. In the 
context of kinetic theory a small shear viscosity therefore implies strong jet quenching 
[279]. In detail the situation is more complicated. One can show that the viscous 
correction to ^(pt) at moderate px is proportional to l/\/q, not 1/q [280]. Also, elliptic 
flow at large p? is mainly a probe of the path length dependence of energy loss [281]. In 
the next section we will see that there are important connections between energy loss 
of light quarks, and energy loss of heavy quarks. 

3. 1. Heavy Quarks 

An important diagnostic of the properties of the QGP is the drag force on a heavy 
quark initially produced in a hard collision during the pre-thermal-equilibrium stage of 
the collision. We will see below that the drag force can be related to the heavy quark 
diffusion constant. Combining experimental constraints on heavy quark diffusion and 
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shear viscosity, which is related to momentum diffusion, provides additional information 
about transport properties of the nearly perfect fluid created in heavy ion collisions. 

The diffusion of heavy quarks in a QCD plasma was first studied by Svetitsky [288]. 
In the following we will follow the arguments presented by Teaney and Moore [286]. 
Consider a small density uq of heavy quarks inside a QGP. In a heavy ion collision, 
heavy charm and bottom quarks are produced in hard collisions between quarks and 
gluons during the pre-equilibrium stage of the reaction. If there is enough time, and if the 
interaction is sufficiently strong, then the heavy quarks can reach thermal equilibrium.* 
In this case the the time evolution of the density can be described by a diffusion equation 

^ = DV*n Q , (54) 

where D is the diffusion constant. We can relate D to the drag force on the quark by 
using a stochastic (Langevin) equation. 

^ = -VdP + UW)) = n5 l3 5(t - t'). (55) 

Here, p is the momentum of the particle, r\z> is the drag coefficient, and £ (t) is a stochastic 
force. In a kinetic theory picture the stochastic force models collisions with quarks 
and gluons in the plasma. The coefficient k is related to the mean square momentum 
change per unit time, 3k = ((Ap) 2 )/(Ai). The Langevin equation can be integrated to 
determine the mean squared momentum. In the long time limit (t ^> r/^ 1 ) the particle 
thermalizes and we expect that (p 2 ) = 3mT. This requirement leads to the Einstein 
relation 

^ = 2^T- (56) 
The relation between t\d and the diffusion constant can be determined from the mean 
square displacement. At late times ([Ax(t)] 2 ) = 6D\t\ and 

T 2T 2 

D — — — ——• (57) 
mrjD k 

The diffusion constant for heavy quarks in a QGP can be determined by computing the 
mean square momentum transfer per unit time. At weak coupling and for approximately 
thermal heavy quarks the diffusion constant is dominated by heavy quark scattering on 
light quarks and gluons, qQ —> qQ and gQ — > gQ. These processes are similar to the 
processes that determine the shear viscosity, except that heavy quarks move slowly and 
the interaction only involves the color Coulomb interactions, whereas shear viscosity is 
sensitive to both electric and magnetic interactions. The leading order result in QCD 
with three light flavors is [288, 286] 

67T 

= g*riog(2T/m D ) ' (58) 

* The charm and bottom quark masses are m c ~ 1.3 GeV and m& — 4.2 GeV, respectively. This 
implies that even if thermal equilibrium is reached, the density of heavy quarks is suppressed by large 
Boltzmann factors cxp(toq/T). We will see below that the expected thermalization times are larger 
than those of light quarks and gluons by a factor mg/T. 
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Comparing this result with Eq. (31) we observe that heavy quark and momentum 
diffusion are related. In the relevant range of coupling constants, and keeping terms 
beyond the leading logarithm, one finds DT ~ 6(r)/s). We note that the heavy quark 
relaxation time rf^ contains an extra factor mq/T compared to the hydrodynamic 
relaxation time t]/(sT). For charm quarks at T = 200 MeV this factor is m c /T ~ 7. 
This implies that even if hydrodynamic behavior is reached very quickly, t]/(sT) ~ 0.2 
fm, charms quarks have barely enough time to equilibrate, m c D/T ~ 8 fm. 

This simple relation between DT and rj/s is broken in the strong coupling limit of 
the SUSY Yang-Mills plasma. Using holographic duality one finds [289, 290, 291] 

D = h^ (59) 

This result implies that in the strong coupling limit there is no bound on the diffusion 
constant, and that even very heavy quarks can possibly equilibrate. It is not clear how 
A should be chosen in order to compare Eq. (59) to experiments in QCD. Gubser has 
advocated a value A ~ 6ir and concludes that D ~ l/(27rT) [292]. 

Experimental information on the diffusion constant comes from the observation 
of single electrons from charm and bottom decays, as shown in Fig. 20 [282]. There 
is data on both the nuclear modification factor Raa = dN Au+Aul 1 ((TAA)da p+ p) , where 
dN^u+Au is the differential yield in Au+Au, (Taa) is the nuclear overlap, and da p+p is 
the differential cross section in p+p, as well as on the elliptic flow parameter V2. The 
nuclear modification factor is sensitive to energy loss, which in turn is governed by the 
drag force. Drag also implies that the elliptic flow of light quarks and gluons induces 
a non-vanishing t> 2 parameter for heavy quarks. Different models for the drag force are 
shown as the dashed and dash-dotted lines in Fig. 20. We observe that these models 
can account qualitatively for the behavior of Raa but fail to reproduce the observed 
flow. This indicates that charm quarks show some degree of equilibration. A Langevin 
simulation with (2irT)D = (4 — 6) provides a qualitative description of both Raa and 
t> 2 . It is interesting to note that this result, combined with the estimate DT ~ 6(r]/s), 
gives (4ir)ri/s ~ (1.3 — 2.0), which is close to the result obtained from the elliptic flow 
of light particles [282]. 

There are a number of uncertainties in this analysis that will be addressed in 
the future. Experiments have not yet been able to determine the flavor of the heavy 
quark, and there are significant theoretical uncertainties in predictions of the charm 
and bottom spectra. Future experiments will be able to identify the flavor of the heavy 
quark. Two contributions in this Focus Issue address transport properties of heavy 
quarks. Meyer [293] presents a lattice study of the Euclidean chromo-electric field 
correlation function. This correlation function is related by a Kubo formula to the 
diffusion constant. Rapp and Riek [294] perform a non-perturbative T-matrix analysis 
of the charm quark diffusion constant and the charmonium spectral function. These and 
other studies, together with ongoing efforts to measure the spectra of identified heavy 
quarks, will shed light on the relation between the different transport properties of the 
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plasma such as diffusion, energy loss, and shear viscosity, and help to characterize the 
degree to which the initial state in a heavy ion collision thermalizes. 

4. Holographic Duality 

Holography is a duality relating quantum field theory (QFT) and gravity. Roughly 
speaking, holographic duality maps the quantum physics of strongly correlated many- 
body systems to the classical dynamics of black hole horizons in one higher dimension, 
replacing quasiparticles with geometry as the salient degrees of freedom. As such, 
holography is very much in the tradition of emergent critical phenomena: when the 
system is strongly coupled, new weakly-coupled degrees of freedom dynamically emerge. 
The novelty is that the emergent fields live in a dynamical spacetime with an extra 
spatial dimension. This extra dimension plays the role of an energy scale in the QFT, 
with motion along the extra dimension representing a change of scale, or renormalization 
group (RG) flow, in the QFT. Holography thus translates problems in quantum many- 
body physics, such as thermodynamics and transport physics, into equivalent problems 
in classical gravity. 

The power of this reorganization is that various phenomena which are easy to 
see, or perhaps even universal, in one presentation may be very surprising in the dual 
presentation, but just as universal. For example, it is a classic result that the physics 
of black hole horizons in general relativity (GR) is largely independent of the details of 
the black hole [295, 296, 297, 298, 299], with small fluctuations of the horizon obeying 
the equations of viscous hydrodynamics [300]. The dissipation of waves in this fluid 
encodes the absorption of energy by the black hole, with the shear viscosity taking a 
universal value determined by basic properties of the Einstein-Hilbert action [300, 301]. 
Holographically, this is a remarkable fact: it tells us that any generic strongly-interacting 
quantum many-body system at finite temperature and density and with a sufficiently 
large number of degrees of freedom per unit volume (necessary for the application of 
holographic duality) can be expected to behave, at low energies, like a nearly-perfect 
liquid, and not like a gas of long lived quasiparticles as one might naively imagine. 

More generally, holography gives us an entirely new way to define quantum field 
theories. Explicitly, holography provides a recipe [302, 303, 304] for using classical 
gravity in (c?+l)-dimensions to compute quantum amplitudes that manifestly satisfy the 
consistency conditions of a d- dimensional QFT (locality, causality, etc.). s In general we 
do not know how to identify the corresponding QFT in terms of more familiar tools, for 
example by specifying a Hamiltonian governing the interactions of a set of well-defined 
quasiparticles. Moreover, as we shall explore in some detail below, the holographic 
description is generally reliable precisely in situations when most traditional techniques 
are not: when A, the typical coupling in the QFT, is strong, A 3> 1, and when N, 



E This recipe is colloquially referred to as the holographic dictionary, a sketch of which is presented in 
Table (4.2.3). 
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the number of degrees of freedom per unit volume, is large, N ^> 1." In such cases 
we can simply take the gravitational description as a constructive definition of the 
QFT. Such holographic QFTs thus define a special subset of the space of well-defined 
QFTs which does not depend on any quasiparticle picture or conventional perturbation 
theory. Importantly, some of the most theoretically and experimentally interesting 
real-world systems manifestly do not have any well-defined quasiparticles upon which 
to base a standard QFT. Holography thus provides an entirely new way to construct 
consistent models of these strongly correlated quantum many-body systems, replacing 
quasiparticles with geometry as the central organizing principle. 

An important corollary is that it is typically more fruitful to use holographic QFTs 
as windows on a general class of phenomena, or to study general concepts in the space 
of QFTs, than to try to exactly reproduce or solve a specific QFT of previous interest. 
Indeed, building "the holographic dual" of one's favorite QFT is generally quixotic, as 
the regimes of validity of holographic QFTs typically exclude the theories of previous 
interest. The classic example is SU(N) gauge theory: QCD is given by N = 3, while 
the holographic dual is classical only when N ^> 1 and weakly-coupled only when the 
't Hooft coupling is also large, A> 1. It is thus futile to try to reproduce QCD exactly. 
Where holography has proven useful, rather, is in studying general properties of SU(N) 
gauge theories at high temperature and densities where many interesting phenomena 
arise which appear to be relatively insensitive to the precise value of N. For example, 
holographic models suggest that any strongly-coupled QGP should behave as a liquid 
whose viscosity is very low and relatively insensitive to the precise value of the coupling, 
in sharp contrast to the large viscosity and strong coupling-dependence predicted by 
weakly-coupled QCD. The low-viscosity of the QGP observed at RHIC (see Sec. 3.5) 
thus suggests that the RHIC fireball is indeed an extremely strongly-coupled quantum 
liquid. 

To be sure, the simplest holographic models are in many ways very different from 
QCD - in the most well-understood example there are no quarks and no mass gap! 
What is remarkable from this point of view is that many of the striking features of 
the simplest holographic models persist even after the inclusion of quarks, a mass gap, 
and other phenomenologically important ingredients. For example, in more realistic 
holographic models of QCD which contain fundamental quarks and display confinement, 
the physics at low energy is again governed by hydrodynamics with an exceptionally 
low and coupling-insensitive viscosity. Such holographic models may thus be treated of 
as computationally tractable toy models which exhibit a rich set of behaviors analogous 
to those observed in the lab, the study of which can reveal qualitative, and sometimes 
even quantitative, general properties of the larger class of strongly-coupled QFTs. 

" More generally, holographic duality relates QFT in <i-dimensions to quantum gravity in (d + 1)- 
dimensions, with the gravitational description becoming classical when the QFT is strongly-coupled, 
as discussed in Sec. (4.2.3). We will focus on regimes where the QFT is strongly coupled and the 
gravitational dual classical, but we emphasize that holography remains true even when quantum 
gravitational (string theoretic) effects become important. 



Strongly Correlated Quantum Fluids 



62 



It is important to stress that while the original discovery of holographic duality [302] 
required esoteric tools such as supersymmetry and conformal invariance,^ the duality 
itself does not depend on supersymmetry or conformality. + In fact, most current work 
on holographic duality involves non-supersymmetric, non-conformal QFTs. Indeed, we 
now have a long list of examples in which each of these constraints is weakened or 
removed (see e.g. [305, 306, 307, 308] and references therein). To date there is no 
example of a violation of a sharp holographic duality. From the current point of view, 
holographic duality is simply a true, if as yet unproven, fact about quantum field theories 
and quantum gravity. 

While our understanding of the duality remains in various ways incomplete, 
holographic model building has already generated novel insights, including predictions 
for the anomalously low viscosity of cold atom gasses, which follow from the same 
universal horizon physics as that of strongly-coupled SU(N) plasmas, and lessons about 
jet quenching and rapid thermalization in the RHIC fireball [184]. Beyond providing new 
tools with which to model specific phenomena, holography has sometimes suggested 
new organizing principles for strongly quantum dynamics, as seen for example in 
semi-holographic models of non- Fermi liquids [309, 310, 311], or in the application 
of classical numerical relativity to study far-from-equilibrium dynamics in extreme 
quantum liquids [312]. 

The goal of this section is to provide a brief introduction to the basic structure 
of holography, first heuristically and then more precisely, as well as its origins and a 
few of its applications. We will also review a few major themes in the holographic 
study of extreme quantum matter and give a broad, if selective view of some of the 
key open questions in the field. Of necessity, we omit many important topics and 
references. For a more detailed introduction, several excellent reviews are available, 
including [184, 310, 313, 314, 315, 316, 317, 20, 318, 319, 320, 321], from which we have 
drawn heavily in preparing various parts of this section. 

4-1. Why Should Holography Be True? Two Heuristic Pictures. 

Why should a quantum field theory (QFT) in <i-dimensions have anything to do with 
gravity in d+1, or vice versa? To give some intuition we describe two heuristic arguments 
that motivate such a holographic correspondence. We first consider a system that 
includes gravity, i.e., objects falling into a black hole, and argue that it should admit 

^ Indeed, holographic duality is often referred to as the AdS-CFT correspondence and also as gauge- 
gravity duality, among many other names. 

+ The role of all this structure is simple: given enough symmetry, it becomes possible to compute 
quantities in an interacting QFT at both weak and strong coupling. Holographic duality was discovered 
when Maldacena pointed out [302] that scattering amplitudes in the maximally supcrsymmctric 4d QFT 
were in fact identical to amplitudes for the low-energy scattering of closed strings off a maximally 
symmetric 5d black hole in string theory. Conformal symmetry still plays an important role in 
holographic QFTs, but no more so than it does in any interacting QFT: conformal symmetry arises 
dynamically at fixed points of the RG. 
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an equivalent description without gravity in one lower dimension. We then examine 
a strictly non-gravitational system, a QFT on a lattice, and argue that it should be 
related to an equivalent theory with gravity in one higher dimension. 

4-1.1. Starting with Gravity: Falling into a Black Hole Imagine standing far from a 
black hole* in (d+l)-dimensions and sending a robotic probe on a straight line trajectory 
directly toward the black hole. To follow the trajectory of the robot, we have attached 
to it a strobe light which flashes once a second. According to the robot, it will itself 
reach the horizon in finite proper time, as measured for example by a clock onboard the 
robot, and will thus have flashed its strobe only a finite number of times before crossing 
the horizon into the black hole interior. 

However, according to you, the distant observer, the story looks very different. As 
the robot falls into the gravitational well of the black hole, the light emitted by the 
flashing strobe must do work to escape the gravitational potential, and is thus red- 
shifted. By the constancy of the speed of light, this means that the flashes arrive to you 
at ever greater intervals. As the robot approaches the horizon, the observed red-shift in 
fact diverges, so that the final flash emitted at the moment of crossing the horizon takes 
an infinite time to arrive at the distant observer; light from flashes after the robot has 
crossed the horizon will never make it out to you. Thus, what you actually see is not 
the robot falling into the black hole, but rather the robot approaching the black hole 
and gradually slowing down and compressing into a thin membrane on the surface of 
the black hole. Indeed, to you, it appears as if the robot has stopped falling toward the 
massive black hole entirely, as if it had stopped responding to gravity at all. 

Now imagine sending in a swarm of robots in an isotropic shell around the black 
hole, so as to increase the mass of the black hole by a small amount corresponding 
to the matter and energy in the collapsing shell. Again, according to you, the distant 
observer, the shell will never actually appear to fall into the black hole. Rather, you 
will see it approach the black hole, slow down, and effectively stop, forming a non- 
gravitating membrane which stretches over the true horizon. As each new bit of matter 
reaches this stretched horizon, it generates a disturbance in the membrane which spreads 
in waves through the stretched horizon as if through a fluid. Remarkably, all of the 
observable behavior of this system can be precisely captured by a non-gravitational 
viscous hydrodynamical model for a fluid on the stretched horizon which lives in only 
d, rather than (d + 1) spacetime dimensions. 

If we describe this same process from the perspective of one of the robotic probes, 
the physics looks very different. In this frame, nothing special happens as the robots 
approach and pass the horizon. On the contrary, the swarm of robots continues to fall 
towards the black hole, interacting with each other and with the gravity of the black 

* For example, in a spaceship following a distant circular orbit of the black hole. We emphasize that 
everything in Sec. (4.1.1) follows from standard results about black holes in GR. For a clear and concise 
review of the physics of black holes in GR; see [322] and references therein. More complete references 
can be found in [323]. 
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Figure 21. Two views of a swarm of robots falling toward a black hole in an isotropic 
shell formation. [Left] A distant observer sees the robots slow down and spread over the 
horizon. According to her, the proper description of the robots is as a non- gravitating 
fluid stretched over the horizon. [Right] An observer falling along with the robots would 
see them continue steadily towards and beyond the horizon, responding as one would 
expect to the black hole's gravitational pull. According to the in-falling observer, the 
appropriate description is general relativity in the bulk. 



hole in the full volume of {d + l)-dimensions. These dynamics can be modeled in a 
straighforward manner with GR in (d + l)-dimensions coupled to the matter and energy 
of the robots and the black hole. 

This leaves us with two very different descriptions of our probes as they approach 
the horizon, one involving the gravitational dynamics of probes falling into a (d + 1)- 
dimensional black hole, the other involving the hydrodynamic response to our probes of 
a d- dimensional fluid. Both descriptions accurately capture the behavior of our probes 
as viewed by two different observers. But there cannot be two facts of the matter 
about the physics of our probes before they cross the horizon: two observers cannot 
accurately observe contradictory events unfold." The gravitational dynamics of a (d+l)- 
dimensional black hole must thus be, in some deeply non-local way, equivalent to the 
hydrodynamics of a d- dimensional fluid. 

4-1-2. Starting without Gravity: Taking the Renormalization Group Literally. Rather 
than start with gravity, let's start with a familiar non-gravitational field theory. Consider 
a system on a lattice with lattice spacing a and Hamiltonian, 



Here, x labels the sites in the lattice, i labels the various operators O l (x) defined at each 
site, and the Ji(x) are coupling constants / sources for the operators^ O l . Note that 

' Explicitly, so long as both observers remain outside the black hole, it remains physically possible for 
them to exchange information and compare their observations. 

n We write sources with index down and operators with index up, for later notational convenience. 




(60) 
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the sources will in general depend on both space and time; we take x to stand in for 
both as is typical spacetime notation in relativity. Given this setup, what we generally 
want to compute is the physics of the ground state and the low energy excitations over 
this ground state as a function of the microscopic coupling constants at the lattice scale. 
In general, however, exactly diagonalizing the Hamiltonian is intractably difficult. 

Kadanoff and Wilson taught us a beautiful approach to this problem: the 
Renormalization Group (RG) [324, 325, 326]. The basic move is to iteratively coarse 
grain the lattice, making the lattice spacing larger with each step, so that at each step 
a single site represents the average of multiple sites in the previous lattice. We then 
tune the couplings so as to preserve the physics of the ground state and the low-energy 
excitations over it, replacing the fixed coupling Ji(x) with a scale-dependent coupling 
Ji(x,u), where u denotes the length scale at which we probe the system. By iterating 
this procedure, we eventually arrive at an effective description of long wavelength modes 
of the system. As explained by Kadanoff and Wilson, the resulting flow of the couplings 
with scale u can be encoded in a beta function which is, remarkably, local in energy 
scale, 

d 

u—Ji(x, u) = 0i(Jj(x, u), u) (61) 

When the beta function can be determined, for example in perturbation theory, this 
renormalisation group approach is extremely powerful. Indeed, even finding the fixed 
points of the RG flow, corresponding to scale-invariant or conformal points governed 
by conformal field theories (CFTs), can be enormously revealing, so much so that we 
often define scale-dependent QFTs with non-trivial RG flows by starting with a well- 
understood CFT and turning on a relevant deformation to generate the desired RG 
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Figure 22. Coarse-graining spins on a lattice a-ld Kadanoff and Wilson. At each step 
we coarse-grain by replacing the degrees of freedom on a block of sites by an average 
value on a single site, rescaling the lattice spacing a and tuning the couplings Ji(x,a) 
in the Hamiltonian H so that the physics of the ground state and low-lying excitations 
remain invariant under the scaling operation. 
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flow. 



However, in many complex, strongly coupled systems, including many strongly 
correlated systems governed by interesting conformal fixed points, the beta functions 
cannot be straightforwardly derived. It is thus tempting to search for a reorganization 
of the RG which determines the correct RG flow of the couplings without requiring an 
explicit calculation of the /3-functions. 

To this end, consider the following recasting of the Kadanoff- Wilson picture. As 
before, consider our series of coarse-grained lattices with coarse grained Hamiltonians 
and suitably tuned couplings, Ji(x,u). Now imagine arranging each coarse-graining of 
the lattice into a stack ordered by scale, so that at each subsequent level of the stack 
the lattice spacing grows, u={a, 2a, 4a,. . . }. By construction, motion down the stack 
reproduces RG flow in the original lattice theory. In this one-higher-dimensional hyper- 
lattice, however, the couplings at each scale, Ji(x,u), look a lot like fluctuating fields in 
a one-higher-dimensional lattice, varying both in space-time and in the scale direction. 

This raises an interesting question. We know that the true solution Ji(x,u) solves 
the /3-function equations for our original lattice theory, but we do not know what the 
correct /3-functions for our lattice model are. Might there be some simple equation 
of motion on this one- higher- dimensional lattice whose solution is also Ji(x,u)7 More 
generally, might there be a conventional local field theory^" whose dynamics encode the 
full, and unknown, /3-functions of the original lattice theory? 

Consider the full RG-stack of lattices labeled by a new RG-coordinate, r, which 
runs from the original ultraviolet (UV) lattice cutoff, r = a, to the deep infrared (IR), 
r — > oo. We seek a simple QFT defined on this one- higher- dimensional space whose 
(rf + l)-dimensional bulk fields, §i(x, r), are in 1-to-l correspondence with the couplings, 
Ji(x), of the underlying lattice theory, with the values of the bulk fields at the UV end 
of the stack, r = a, determined by the microscopic couplings, 



As a consequence, the bulk field $(x,r) must have the same charges, tensor structure 
and other quantum numbers as the corresponding coupling,-'- Ji(x). Thus to each scalar 
operator 0(x) is associated a bulk scalar field Q(x,r) such that a)0(x) is a scalar 
operator we can add to our Hamiltonian, to each current operator J^(x) is associated a 
bulk vector field, A^x, r) such that A^(x, a)J' ll (x) is a scalar, and so forth. In particular, 
to the canonical stress-tensor T^ v (x) is associated a canonical spin-two field in the bulk, 
9^u(x,r). 

' By a conventional local QFT we mean a QFT with a finite number of local fields governed by a 
Lagrangian with canonical quadratic kinetic terms and local interactions. One can certainly imagine 
other possibilities, but this structure turns out to be particularly useful in what follows. 
$ Recall that our couplings are defined as the coefficients of the associated operators in the Hamiltonian, 
and so have definite quantum numbers. Note too that, at fixed points of the j3- functions where the 
system is scale-invariant, the dimension of an operator is also a quantum number. This is also true in 
perturbations around such conformal fixed points. 
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Figure 23. Holography promotes the tower of coarse-grained lattices to a one-higher- 
dimensional lattice, with the RG scale recast as a spatial dimension and the running 
couplings Ji replaced by dynamical fields <&i which asymptotically approach the UV 
couplings at the UV top of the stack and which satisfies a gravitational equation of 
motion in the bulk of the stack. 



What can the Lagrangian of this bulk QFT be? It is not the original lattice 
Hamiltonian - the fields of the bulk QFT encode the dynamics of the couplings in 
RG space, not of the operators on the original lattice. A priori, the natural thing 
to do is to write down the most general effective field theory allowed by the fields 
and symmetries of the system. However, a simple argument greatly constrains the 
possibilities. Any QFT defined by perturbation from a fixed point CFT as discussed 
above includes among its operators a canonical stress-energy tensor, T^ u , arising as 
the Noether current corresponding to the space-time symmetry group. The bulk QFT 
should thus include a corresponding canonical spin-2 field, But by the Weinberg 
and Weinberg- Witten theorems, s and assuming the existence of a Lorentz-invariant 
continuum limit, any spin-2 field must either decouple at low-energy, which would imply 
that sources of momentum or energy in the original lattice theory do not affect the 
system, or it must couple universally according to the equivalence principle. In other 
words, the bulk field theory must be either topological or be a theory of gravity which 

s Weinberg [327] argued that any Lorentz-invariant theory of a spin-2 field must respect the equivalence 
principle or suffer IR pathologies; the Weinberg- Witten theorem [328] further states that any massless 
spin-two field must either be the graviton or suffer one of a host of pathologies, for example not having 
a conserved stress tensor. For an illuminating discussion of these constraints; see [329] and references 
therein. 
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reduces, at long wavelength and low energy, to GR. 

If all this can be done consistently, we would appear to have two distinct descriptions 
of our system: the original lattice theory with lattice spacing a and sources Ji(x) which 
satisfy first-order /3-function equations; and a gravitational description living in a one- 
higher-dimensional bulk whose fields, $j(x,r), are in 1-to-l correspondence with the 
sources Ji(x) but which satisfy conventional second-order local field equations. In other 
words, we have two ways to compute correlation functions in our system: we may 
either work in the lattice theory and solve the QFT problem; or we may work in the 
gravitational description and solve a GR problem. 

In order for these two descriptions to have any chance of being equivalent" the bulk 
theory must satisfy two among many rather unusual properties, both of which turn out to 
be basic properties of black holes. First, the entropy as computed in the two descriptions 
must be the same. In the original lattice QFT, the entropy will generically be extensive, 
with the entropy Sqft of the rf-dimensional QFT in a region R<i scaling linearly with 
its spatial volume, Sqp^{Rd) oc V(Rd). As a result, the entropy as measured in the 
one-greater-dimensional bulk gravity, Sqr must be sitfr-extensive, scaling as the area 
bounding a region, Sgr(R<i+i) oc A(Rd+i) oc V(dRd+i)- This is a basic feature of any 
theory of gravity:' the entropy in a volume is bounded by the entropy of a black hole 
that fits inside that volume [330]. Since the entropy of a black hole is proportional to 
the surface area of its horizon, 5bh = zAff, this tells us that the entropy of any region 
in a theory of gravity is bounded by its surface area, Sgr(-Rcz) < \A(R d ). Indeed, the 
fact that the entropy of a black hole can be entirely associated with its surface, and that 
as a result the entropy in gravity is necessarily sub-extensive, is the origin of the term 
holography [331, 332]. + 

Second, the /3-function equations in the boundary lattice theory are first-order in 
the scale r, while the gravitational equations of motion in the bulk are second-order in 
r. Physically this means that specifying the sources in the lattice QFT at one scale 
completely determines the coupling at all other scales, i.e., completely fixes the RG 
trajectory. By contrast, in the gravitational system we have two solutions to our second- 
order equation, so we must further specify the derivative of the couplings with scale to 
uniquely fix the RG trajectory. How can these be equivalent? Somehow it must be the 
case that the gravitational system automatically selects one of the two solutions of the 

" By equivalent we mean that their partition functions are exactly equal, so that both descriptions 
contain precisely the same information, as we shall spell out in detail in Sec. 4.2.3. Note that this 
implies their free energies and other thermodynamic data must similarly, and identically, match. 
" Various other theories, including topological theories, share the property of sub-extensivity. However, 
by Weinberg- Witten, to have both an area law and an interacting spin-2 field, gravity appears to be 
the only option. 

+ The term holography is used by analogy with familiar holograms, and is meant to convey the 
surprising fact that the information about the (d + l)-dimensional bulk spacetime is being reliably 
encoded in the data of the QFT living on the d-dimensional boundary. The analogy is in various ways 
quite apt. For example, information which is local in the (c?+l)-dimensional bulk is encoded non-locally 
in the d-dimensional boundary. However, one should not take the analogy too literally; for example, 
the encoding here involves quantum-mechanics and gravity. 
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bulk equations as physical. 

It is again the physics of black holes which explains how this can happen. Imagine 
solving a simple wave equation in the neighborhood of a black hole horizon. In principle 
there are two solutions to the second-order wave equation. However, since things can 
fall into a black hole horizon but nothing can escape, it is natural to organize the 
solutions near the horizon into in-going and out-going waves. Moreover, when working 
in Euclidean signature as is done for example in computing thermodynamic data, only 
the in-going solution is regular at the horizon; the out-going solution is always irregular 
at the horizon. Thus, while the gravitational field equations are indeed second-order, 
the presence of a black hole horizon effectively adds a second boundary condition, so 
that we again need only specify a single boundary condition to completely determine 
the solution. Since the location of the horizon encodes the thermodynamic variables 
(temperature, entropy) of the system, the physical meaning of this second boundary 
condition is to specify the state of the QFT. 

These heuristics suggest a connection between quantum field theory in d spacetime 
dimensions and quantum gravity in d + 1 spacetime dimensions, with the fields in the 
gravitational system exactly paired with sources in the dual QFT and with the extra 
spatial coordinate in the gravitational bulk playing the role of an RG scale. The next 
sections will describe a precise and computationally effective theory of this connection 
known as holographic duality, and review some of the lessons gleaned from its study to 
date. 

4-1.3. Coda: The Surprise of Locality It is worth pausing to contemplate just how 
remarkable this proposal truly is. Consider a lattice theory in ^-dimensions which has, 
in its Hamiltonian, a real parameter q. You might think of q as a global charge, or as 
a momentum, or something else entirely, which can in principle be varied adiabatically. 
All correlation functions in the theory are thus functions of <i-dimensional spacetime, 
(x,i), as well as of our parameter, q. 

We might be tempted to take a stack of such lattices, each with a different value of 
q, call the "space" labeled by (x, t, q) a (d+ l)-dimensional spacetime, and claim that 
the resulting tensor product of lattice QFTs is well described by a (d + l)-dimensional 
QFT living on this (d + l)-dimensional spacetime. 

This would be wrong. Being (d + l)-dimensional means that all observables 
should transform appropriately under translations, rotations, accelerations, and general 
coordinate transformations amongst x, t and q. They should also respect (d + 1)- 
dimensional causality and locality. But that is absurd! What does it mean to rotate 
between the x direction and the q parameter? Meanwhile, there is absolutely no reason 
for correlators of two operators which live on lattices with very different values of q 
to vanish outside some q — t light-cone, or for operators from distant g-slices to have 
vanishing equal-time commutators in whatever state you choose to call the ground state. 
Locality and causality are almost sure to be violated, and badly. This is not a (d + 1)- 
dimensional system. 
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To summarize, saying that the system is well-described by a local, causal, unitary 
(d + l)-dimensional QFT implies a long list of constraints on the (d + l)-dimensional 
correlation functions. These translate into constraints on the original lattice theory and 
on the g-dependence of its correlation functions. It is by no means obvious that any 
such stack of lattices can ever have a causal, local (d+ l)-dimensional description. From 
this point of view, holographic duality, in which the role of q is played by the energy 
scale and in which the (d + l)-dimensional description is in fact local, causal, and even 
classical when the d- dimensional QFT is very strongly interacting, is nothing short of 
astonishing. Precisely why, when, and how this locality arises is perhaps the deepest 
question in holography* 

4-2. Essential Holography 

Holographic duality [302, 303, 304] is a precise equivalence between certain <i-dimensional 
quantum field theories and (d + l)-dimensional gravitational theories which provides a 
sharp realization of the heuristics described in Sec. 4.1. In particular, all the basic 
features noted above will reappear below: the fields in the bulk correspond to the 
couplings in the QFT; the RG flow of the QFT is encoded in the radial evolution of the 
gravitational theory along an extra dimension; and black hole horizons play a key role. 
Of particular importance is the precise geometry of the emergent spacetime, as well as 
the precise relationship between observables in the QFT and those of the gravitational 
dual. We will explore a few simple examples along the way. 

For clarity of presentation, we will focus our discussion on the simplest possible 
QFTs from the point of view of the RG, i.e., Lorentz-invariant conformal field theories 
(CFTs) which are fixed-points of the RG. Since the RG maps to radial evolution in the 
gravitational dual, QFTs corresponding to RG fixed-points must be dual to geometries 
which are translationally-invariant along the emergent radial direction up to overall 
rescalings; adding Lorentz-invariance then entirely fixes the dual geometry to be Anti- 
de Sitter space (AdS), whose geometry we will now describe." 

4-2.1. Gravity and Matter in Anti-de Sitter Space AdS^+i is a homogeneous, isotropic 
spacetime with d+1 spacetime dimensions and constant negative curvature whose metric 
can be conveniently written as,^ 

L 2 

ds 2 = — [-dt 2 + dx 2 + dr 2 ] , (63) 

* See for example [333] for a discussion of locality in holographic duality, and [334] for an interesting 
spin of the problem. 

' The first examples of holographic duality involved such AdS/CFT dual pairs, and hence holographic 
duality is often referred to as AdS/CFT. Just as we can always construct a general QFT by perturbing 
away from a CFT fixed point, wc can extend such AdS/CFT dualities to more general QFTs and 
geometries by perturbing both sides in corresponding ways, so we do not lose too much by focusing on 
these examples. 

" For the remainder of our discussion of holography wc will work in natural units in which h = 1 and 
c=l. 
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where x are (d-1) spatial coordinates, t is a timelike coordinate, and r is a final "radial" 
spatial coordinate. The constant parameter L, known as the AdS-radius, sets the radius 
of curvature, with the Ricci curvature scalar given by the constant 

R = -^A. (64) 

The space is thus weakly curved when L is large and strongly curved when L is small. 
At fixed values of the radial coordinate r*, the metric reduces to the d- dimensional 
Minkowski metric on flat space-time rescaled by L 2 /r 2 , 

dsl = ^ [-dt 2 + dx 2 ] , (65) 

It is thus useful to think of AdS as a stack of slices of flat space, with r controlling the 
physical scale on that slice. The slice at r — > is referred to as the boundary. Note 
that the volume element of the boundary slice at r — > diverges, while the proper 
distance to the boundary also diverges. Thus the boundary at coordinate r — >■ should 
be understood as the set of points at spatial infinity. ^ 

Importantly, the full geometry is invariant under not just translations, rotations 
and Lorentz boosts along each slice, but also under general scaling transformations, 

(r, x, t) — > (ar, ax, at) (66) 

as well as rotations and boosts that mix r with the other dimensions. The full isometry 
group of AdSrf + i is in fact identical to the conformal group in d space-time dimensions.* 

It is often useful to think of AdS as a harmonic trap for gravity. More precisely, 
the constant negative curvature of AdS acts as a harmonic trap such that if you sit 
inside AdS at an arbitrary point and throw a ball, it will inevitably fall back to you in 
finite time. Indeed, if you fire out a photon, that photon will run away to the boundary 
at r = and return, again in finite observed time. s And yet while nothing, not even 
gravity, gets out of the AdS box, the spacetime remains completely homogenous and 
isotropic. AdS thus acts as a peculiarly graceful IR regulator for gravity which breaks 
none of the symmetries of flat space. 

All of this implies an intimate connection between the radial coordinate in the bulk 
of AdS and spatial scales along the boundary: probing short distances (or high energies) 

' In computations, it will often be useful to regulate this divergence by introducing a cutoff at r = a, 
with a <C L, removing the regulator, a — > 0, only at the end of all computations. 

* Note that our regulator at r = a explicitly breaks conformal symmetry but does not break space-time 
rotational or translational invariance in x and t, as a more conventional UV regulator such as a hard 
momentum cutoff would. Matching regulators between bulk and boundary turns out to be quite subtle, 
c.f. [335]. 

s Unlike the electromagnetic harmonic traps used in trapping cold atoms, AdS is homogeneous and 
isotropic: there is no center to which all objects return. Rather, the curvature ensures that any two 
initially parallel trajectories always curve towards each other, no matter where they begin, and then 
oscillate. 
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along the boundary corresponds to probing the bulk only near the boundary at r = 0, 
while probing long distances (or low energies) along the boundary corresponds to data 
deep in the interior of AdS. Roughly, then, we can think of the region near the boundary 
of AdS as associated with the UV physics of the boundary, and of the deep interior of 
AdS as associated to the IR physics of the boundary 




r 



Figure 24. Hanging rope in AdS. When the ends are close together in the spatial 
direction, x, so that the rope probes short distances and high energies in the QFT, the 
stretches only slightly into the bulk of the AdS spacetime. When the ends are far apart, 
probing long distances and low energies in the QFT, the rope stretches deep into the 
AdS bulk. Thus is the boundary of AdS associated to the UV of the QFT, while the 
deep bulk is associated with the IR. 

To visualize this association, consider a fixed tension rope whose ends are glued to 
the boundary at r = 0. Due to the constant negative curvature of AdS, the bulk of the 
rope is drawn into the bulk of the AdS spacetime, forming an arc drooping away from 
the boundary. If the ends of rope are held close together - probing short distances or 
high energy in the 4d boundary - the rope barely hangs into bulk. If we instead pull the 
ends of the rope far apart - probing long distances or low energy in the 4d boundary - 
the rope dips much further into the bulk. In this way, probing the IR of the boundary 
corresponds to probing the deep interior of the AdS, r 4 oo, while probing the UV 
of the boundary corresponds to focusing on the near-boundary region of the geometry, 
r — > 0. Note that the near-boundary cutoff at r = a <C 1 thus acts as a UV regulator, 
while the horizon at r = th provides a natural IR regulator. 

Importantly, AdS is a solution to the equations of motion of a generic Wilsonian 
action for the metric 

/Gravity = I (-2A + R + c 2 R 2 + c 3 R 3 + ...). (67) 

Here, Gn is the Newton constant, g is the determinant of the spacetime metric, 
g = det(g^), R is the Ricci curvature scalar built out of two derivatives of the metric, 
R ~ ddg , A is a cosmological constant (a.k.a. the tension of the vacuum), and the 
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. . . represent all other scalars one can build out of the metric and its derivatives. The 
R term in the action plays the role of a two- derivative kinetic term for the metric, 
g^ v . GR is defined by retaining only the lowest dimension kinetic term, R, giving the 
Einstein- Hilbert action whose equation of motion is the Einstein equation, 

R^u - ^g^R = -J^g^v (68) 

The c n R n + . . . terms then represent higher- derivative corrections to GR. The reader 
can verify that the AdS metric solves the Einstein equation (68) with the AdS-radius L 
determined by the cosmological constant, A, as A = — ■ O ne can then show [336] 
that AdS continues to be a solution even when we include generic higher-curvature 
corrections: all that changes is the precise relationship between A and L. 

Importantly, the Ricci scalar in AdS scales as R ~ j? (c.f. (64)). As a result, 
we may neglect higher-curvature corrections in AdS when L 2 is large compared to 
the appropriate powers of the dimensionful couplings c n . For example, if the higher- 
derivative operators are generated by quantum gravity effects, the dimensional couplings 
are controlled by the Planck length, £ p . Quantum corrections to the metric can thus be 
neglected in AdS when L is large compared to the Planck length, j- 3> 1. Similarly, 
GR receives corrections from string theory even at the classical level. Such corrections 
become important at a characteristic length scale known as the string scale, i s . In AdS, 
these stringy corrections can be neglected so long as j- ^> 1. Our generic Wilsonian 
action for quantum gravity in AdS is thus well-approximated by classical GR when the 
AdS-radius is large compared to all scales in the problem, 

L L 

F »l, 7T» 1 ' ( 69 ) 

where the . . . denote other sources of corrections to GR which may kick in at other 
length scales. 



Black Holes in AdS If pure empty AdS is a ground state for gravity, finite-temperature 
states correspond to black holes inside AdS. The simplest such asymptotically- AdS black 
hole is the AdS-Schwarzchild black brane, 



ds 2 = ^r 



-f{r) dt 2 + dx 2 



-dr' 



with emblackening factor 



fir) 



i - 



(70) 



(71) 



H 



Near the asymptotic boundary at r — >■ 0, / — > 1, so this metric is asymptotically AdS^+i. 
At r = rn, however, / — > 0, signaling the presence of a black hole horizon. This horizon 
in turn shields us from a physical singularity at r — > oo by ensuring that nothing which 
is inside the horizon, and thus sensitive to the singularity, can ever escape to influence 
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events in the rest of the spacetime. Since the entire solution is translationally invariant 
in the (d — 1) spatial directions, x, this black brane is not a compact object, but rather 
extended in all directions other than r. 

Importantly, a classical black hole is a thermodynamic object with a definite 
temperature, energy and entropy, as shown by Bekenstein and Hawking (see [330] 
and references therein). If this were not the case, we could violate the third law 
of thermodynamics by throwing all of our waste heat into a black hole and thus 
build perfectly efficient Carnot engines. Black holes in AdS similarly carry definite 
temperatures, energies and entropies, though the details are a bit different from the flat 
space case [337]. In the case at hand, the Hawking temperature T, energy density e and 
entropy density s of the black brane are then, 

rp d = d-l ( L\ dl = _J_ (L 

"47rr H ' " 167rr£ \£ P J ' * Ar^ 1 \£ p 

where £ p is the Planck length defined by Gn = The ratio f- thus measures the 

AdS-scale in Planck units. 

Physically, these thermodynamic relations are telling us that if we increase e by 
throwing additional mass or energy into the black brane, the horizon swells outward 
toward the asymptotic boundary (rjj — > 0) and the black brane heats up. The specific 
heat of this black brane is thus positive. This means that black holes in AdS can come 
to equilibrium. Again, this can be traced to AdS playing the role of a harmonic trap 
for gravity: any radiation that the black hole evaporates away will, in finite time, fall 
back into the horizon." 




Charged Black Holes in AdS More generally, we can add charge to our black brane by 
adding a non-trivial Maxwell field, 



- f(r)dt 2 + dx 2 + -rhdr 2 
ds 2 = L 2 ——— f M , A = A t (r)dt , (73) 



with emblackening factor 

/ = 1 - M r d + Q 2 r 2 ^ , (74) 
and electromagnetic scalar potential^ 

(75) 




" This is very different from black holes in asymptotically flat space, whose specific heat is negative 
and which never come to equilibrium in the absence of external forces. Black holes in AdS thus behave 
like familiar systems in the thermodynamic limit, while black holes in flat space behave more like finite 
volume sub-systems. 

^ Note that we are working with units in which h = 1 and c = 1. For example, M and Q scale as energy 
and charge densities on the (non-compact) horizon, M ~ E/Vd-i ~ L~~ d and Q ~ l/Vd-i ~ L~( d ~ x ' , 
as required by Eq (74). Similarly, both the d+l-dimensional scalar potential, A t , and the d-dimcnsional 
chemical potential, fi, scale as 
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where ^ = f r H d ~ 2 and C = y^l. 

This metric and gauge field together extremize the Einstein-Maxwell action, 

1 



I ME 



16nG 



N 



drxy/=g 



I 2 

-2A + R- —F 2 

4e 2 



(76) 



again with A = ~ d f I y ) . In this geometry, the horizon lies at the radial position r = r^ 
implicitly defined as the value of r where f(r) vanishes. M and Q then determine the 
Hawking temperature of the horizon, 

T - " ^«- 2 ) • < 77 > 

as well as its energy, entropy and charge densities, 

d-lfLY" 1 1 {LY- 1 ^d-lfLY" 1 



e = M wUJ ' ,= i^UJ • ' =0 &5rUJ • (78) 

It is then straightforward to check that these variables satisfy the first law of 
thermodynamics, de = Tds + \i dp, with \i playing the role of the chemical potential. 

A curious feature of these charged black holes is that the zero-temperature limit 
has finite entropy density. To see this, note that we can tune Q to make T vanish, 
Q — > Q* = J rfZf^r- I n this limit, r H , and thus the entropy density, remains non-zero. 
This finite zero-temperature entropy represents an enormous degeneracy of the ground 
state of this black hole,* and correspondingly a large number of potential instabilities 
of these black hole solutions; these turn out to play an interesting role in holographic 
descriptions of quantum phase transitions. For our present purposes, let it suffice to say 
that this is an interesting feature of these black branes which will crop up from time to 
time in the following. 



Matter in AdS: The Effective Action as a Functional of Boundary Conditions The fact 
that AdS behaves like a homogeneous harmonic trap has important consequences when 
we consider matter fields in AdS. Consider for simplicity a scalar field <3> with mass m 2 
in the background of an uncharged AdS black brane (70) with classical action, 

hoc J d d+ \ y/=g (-~ (<9$) 2 - ^$ 2 + . . .) , (79) 

where the . . . denote interactions with other matter fields which we will for the moment 
neglect. Working for a moment with a simple plane wave with momentum k in the 

+ Here = d^A^ is the totally antisymmetric electromagnetic field strength tensor (E l = —F - 1 
and B i = -e ijk F jk ) and F 2 = F^F^ = (E 2 -B 2 ) is the Maxwell kinetic term. 

* The precise microcanonical counting of this ground-state entropy was a long-standing puzzle in 
quantum gravity which was first solved, for a specific set of black holes in string theory, by Stromingcr 
and Vafa in 1996 [338]. A similar counting can now be performed for a much larger set of black holes, 
and in fact played a very important role in the discovery of holographic duality. 
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boundary dimensions, <3>(x, r) 

_ $( r ) e «-* j the wave equation following from the above 

action becomes 

r 2 f$"{r) -r[r/'-(d-l)/]$'(r) - [k 2 r 2 +m 2 L 2 ] $(r) = . (80) 

This equation has two simple degenerate points, one at the boundary where r = and 
the other at the horizon where / = 0. Let's study a general solution in the neighborhood 
of each. 

Near the boundary at r — > 0, where / — > 1, the general solution of this equation 
reduces to, 

$ ~ ^ d -A(k) r d ~ A + (fi A (k) r A + . . . , (81) 

where the scaling dimension A is determined by" A(A — d) = m 2 L 2 . So long as 
m 2 > —d 2 /4L 2 , the scaling dimension is real. This Breitenloner-Freedman (BF) bound 
tells us that a small negative m 2 does not lead to an instability in AdS as it would in flat 
space - instead, the would-be instability is lifted by the harmonic potential generated by 
the AdS curvature. Note that in the flat space limit, L — > oo, this window of "allowed 
tachyons" disappears. So long as the BF bound is satisfied, specifying a solution reduces 
to specifying the two radial integration constants, (fid- a and 0a- 

A key fact about this asymptotic behavior is that, so long as m 2 L 2 > 1 — the 
mode (fid- a is non-normalizable according to the natural norm on a constant-time spatial 
slice, £ t , 

($!, $ 2 ) = -i f dzdx^T~gg u ($^$ 2 - $^$0 . (82) 

Varying this non-normalizable mode thus corresponds to a large (divergent) change in 
the action. We must thus fix the non-normalizable mode to make the variational problem 
in the bulk well-posed, for example by specifying a boundary condition for the bulk field 
$ near the boundary, r — > 0. This leaves us with a one-parameter family of solutions 
labeled by the remaining integration constant, a , which can be chosen arbitrarily. 

We thus draw a general conclusion for studying matter fields in AdS: in order for 
the least-action principle to be well-defined in AdS, we must fix the values of the non- 
normalizable modes of all fields. This is typically done by fixing boundary conditions 
for the non-normalizable modes of all fields at the AdS boundary, r — > 0. 

Near the horizon at r#, where / ~ — r) — > 0, the equation of motion again 
degenerates. Here, however, the physics is rather different. Physically, this equation 
has two kinds of solutions: one which is regular, so that the degenerating /<£>" term 
is negligible; and one which is irregular at the horizon, so that the /<&" term is not 
negligible. But why are half the solutions irregular? The reason was alluded to above: 
near the horizon, all waves can be expressed as a superposition of ingoing and outgoing 
waves. To see that this is precisely what is going on, let's study our scalar equation at 
non-zero frequency. Recalling that the magnitude of the d-momentum k in the metric 

" Similar scaling behavior is obtained for any set of matter fields, with the their resulting scaling 
dimensions A determined by their precise spins, masses and interactions. 
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(70) is k 2 = -jW 2 + r 2 k 2 , the dominant terms in the scalar equation near the horizon 
are 

fd r (fd r $) +w 2 $ ~ 0. (83) 

It is useful to switch to the near-horizon coordinate p defined by fd r = d p , i.e., 
e~P ~ (1 — r) 1 " 11 , in terms of which the horizon lies at p — > oo. In terms of this 
coordinate, the general solution can be expressed as a superposition of ingoing and 
outgoing waves as 

$ = $ in (p)e-^ (t - rf + $out(p)e-^ ( ' +p) . (84) 

For future reference, note that the in(out)-going solutions satisfy the first-order 
constraint at rjj, 

f <9 r $ in = iou<S> m , / <9 r $ out = -iu$ out • (85) 

Note that, in terms of the original radial coordinate, r, the phases of both of these 
time- dependent solutions accumulate near the horizon. One can check that, in the limit 
of zero-frequency, the in-going wave is regular at the horizon while the out-going wave 
is irregular, matching our zero-frequency analysis above. Note, too, that this suggests a 
natural prescription for analytic continuation to Euclidean time: we should choose the 
continuation so that the ingoing wave is regular at the Euclidean horizon, e luJp — > e~ UJEp , 
i.e., we should define u = iuje, where oje > 0. 

Again, this turns out to be a general result: in the presence of a black hole horizon, 
the Laplacian in any tensor representation takes the form □ ~ r 2 fd 2 + "^j- . . . , so 
we must impose a boundary condition on our solutions such that they correspond to 
regular, in-falling modes at the horizon. 

Thus, while our equation of motion is second order, there is only one linearly 
independent solution which is regular everywhere in the bulk and in-going at the 
horizon. Our gravity problem in AdS, in the presence of a black hole horizon, becomes 
effectively a first-order problem.^ Explicitly, fixing the non-normalizable mode (fid-A 
at the boundary and imposing in-falling boundary conditions at the horizon completely 
determines our solution for $ throughout the bulk of AdS, and as a result determines 
0a- However, in order to actually compute a given (fid- a we must solve our problem 
not just at the boundary but all the way through the bulk to the horizon, where we 
impose the appropriate boundary conditions. The relation between (fid- a and 0a is thus 
also determined by IR physics, not just UV physics. 

It is illuminating to see this work in detail in an analytically solvable example. 
Consider again our scalar field $ but now in pure AdS, i.e., with no black hole. This 
corresponds to studying our dual QFT at zero temperature and zero density. Expanding 
in plane waves, e"'^'"' 1 ' , the general solution to the bulk equation of motion can be 
found in closed form, 

(86) 

tt A similar story holds for fermions, though the details differ because the Dirac equation is naturally 
first-order [339]. 
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where v = A — |, k = \/ uj 2 + A; 2 , _?r„ and J v are modified Bessel functions, and reg 
and 0i rrcg are the two integration constants. Since I u {nr) ~ e Kr near the "AdS horizon" 
at r — > oo, regularity requires that 0i rrcg = 0. We thus have just a single integration 
constant, reg , which fixes the overall normalization of the in-falling mode in $. Near 
the boundary, the two asymptotic integration constants (frd-A and 0a i n (81) are not 
independent. We can see how they are related by setting 0i rreg = in (86) and expanding 
the solution near the boundary. After a little algebra, this gives, 

^ = 2 4r ( "A A> f) ( " 2 + * 2)A "^ A - (87) 

Thus, once we impose regularity on our bulk solutions, we need only one boundary 
condition in the UV at r — > to fully specify a solution of the 2 nd order bulk equations 
of motion. Meanwhile, we know we need to fix the non-normalizable mode <pd-A a t the 
boundary so that our variational problem is well posed. Once we impose regularity in 
the bulk, i.e., in-falling at the horizon, the value of a is completely determined. 

This teaches us an important lesson. Suppose we want to study the partition 
function for our gravitational theory in AdS, integrating over all bulk fields, <3>(x,r). 
To make the variational problem well-posed, we must specify a boundary condition for 
each of the bulk fields at the boundary of AdS, as discussed above, 

<pd-A{x) = limr A ~ d $(x,r) , (88) 

r— >0 

where the factor of r A ~ d should be thought of as the appropriate wavefunction 
renormalization of the bulk field $(x, r). This means that the partition function in 
AdS, and thus the effective action I\ds — — hi^AdsL is a functional of the boundary 
conditions for all bulk fields, 

Z A ds[<f>d-A(x)} = Z AdS [$[<^_ A (z)]] . (89) 



What One Means by 11 A Theory of Gravity in AdS" In discussing holography, we will 
regularly refer to classical gravity, quantum gravity and perturbative gravity in AdS, so 
it's worth taking a moment to define our terms here. 

By quantum gravity in AdS we mean a complete quantum description of the 
gravitational system, i.e., string theory in AdS. In general, string theory is a rich and 
complicated quantum theory which has no simple perturbative classical description. 
However, in some cases closed strings have a tractable semiclassical low-energy expansion 
involving a metric, g^, minimally coupled to a host of dynamical matter fields, 
governed by an effective action of the form, 

Igr = ~T7^~Fi / d d+1 x (-2A + R + ... + £ mattcr ($ 4 )) . (90) 

167rGAr J 

Here C mattCT (^i) is the Lagrangian density for the matter fields and R is the Ricci 
curvature scalar, which is roughly two derivatives of the metric. ^[g~R plays the role of 
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a kinetic term for the metric in GR. The . . . then represent an infinite tower of higher- 
derivative terms involving the metric, for example R 2 , i? 4 , etc, which represent closed 
string corrections to GR. Since these higher-derivative terms are also higher-dimension 
irrelevant operators, they must appear with dimensionful couplings. In general these 
couplings are controlled by two independent length scales: £ p , the Planck length, which 
controls quantum corrections to the dynamics of classical strings, and £ s , which controls 
classical "stringy" corrections to the dynamics of point particles. 

When expanding the theory around AdS with AdS-length L (and thus with Ricci 
curvature R ~ L~ 2 ), we thus have two dimensionless parameters controlling the 
various possible corrections to GR: which tells us whether quantum corrections are 
important, and ^, which tells us whether stringy corrections are important. Note that 
ensuring that the curvature is weak in Planck units, <C 1, does ensure that the system 
can be treated classically, but does not tell us that the gravity must be exactly GR: the 
curvature of our AdS may still be large compared to the string scale, so higher-curvature 
corrections to the Lagrangian of GR may not be negligible. For pure classical GR to be 
a good approximation, two independent conditions must thus hold, 

^«1 and (91) 

Lj Lj 

i.e., the curvature must be small compared both to the Planck scale and to the string 
scale. 



4-2.2. Field Theories in Flat Space Our next job is to define the d- dimensional QFT 
we want to study. A canonical way to do so begins by specifying a list of local operators 
Oi labeled by their Lorentz structure, their charges g; and their scaling dimensions Aj, 
all defined at some UV fixed point. To generate an RG flow of interest, we perturb away 
from this fixed point by turning on appropriate sources J;. The basic observables of the 
theory are then correlation functions of products of local operators, 

(0 1 (x 1 )...O n (x n )). (92) 

These can be conveniently encoded in terms of the quantum generating functional, 

ZqftIMx)} = ( e /<^(^M), (93) 

where the Ji(x) represent a set of sources and couplings for the operators Oi. Note that 
the scaling dimensions and tensor structure of the sources Jj are completely determined 
by those of O^. This allows us to express correlation functions as derivatives of the 
partition function, 

(O,(x0...O,(x0) = ^, h f q "'^ )1 , , , (94) 

where the restriction |j. =0 means evaluate the final result with all remaining sources 
turned off. Solving the theory then requires one to compute the partition function, 
Z QFT [Ji(x)}. 
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Canonical Example: SU(N) Yang-Mills at Large N The classic example of a field 
theory with a well-understood and controlled holographic dual, to which we will appeal 
below, is a special version of SU (N) Yang-Mills gauge theory in flat 4d spacetime called 
the M = 4 theory. Here M refers to the amount of supersymmetry enjoyed by the theory; 
M = 4 is the most one can have in 4d without including gravity. For our purposes, the 
role of supersymmetry is nothing more than a way of turning off quantum mechanics 
without totally trivializing the theory: with M = 4 supersymmetry, some quantities 
do not receive quantum corrections beyond one-loop, and so can be computed at weak 
coupling and reliably extrapolated to strong coupling. In particular, the 6-function of 
the theory can be computed exactly and is identically zero for all values of the coupling. ^ 
The basic ingredients of the Af = 4 theory are a gauge group, G = SU(N), a 
gauge field, A^x), transforming in the adjoint of G, six scalars also transforming 

in the adjoint of G and further enjoying a global 5*0(6) symmetry, and large set of 
fermions. We can safely neglect the fermions for now, as their only role is to ensure 
supersymmetry. The Lagrangian is 

£ym = ~4^ ( Tr ^ + Tr |jD0/|2 + • • •) ' (95) 

where g\ M is the gauge coupling, is the gauge-covariant derivative, . . . denotes 
fermionic terms we neglect, and all fields are canonically normalized with an overall 
factor of — -j-? in front of their kinetic terms. 

What are the good operators of this theory? We may suppose that any good 
gauge-invariant combination of the fundamental fields is a perfectly good observable. 
However, since we have a global symmetry, it is useful to write this list in terms of 
50(6) irreducible representations. For example, the following three Lorentz scalars, 

O = ^Tr0?, I j = Tr<j> {I </>j ) , 0» = TrJ», (96) 

I 

transform as a scalar, a symmetric 2-tensor and a vector, respectively, under the global 
50(6), with the latter corresponding to a charge current operator. In the weak-coupling 
limit, their scaling dimensions A,; are just their engineering dimensions which can be 
read off the Lagrangian. * 

t Our ability to compute various quantities at both strong and weak coupling in the Af = 4 theory was 
key to the original discovery of holographic duality, and is the reason this example is both illuminating 
and canonical. However, that is the extent of the role of supersymmetry, to provide unusually simple 
and tractable examples. It is no more a necessary feature of holography than spherical symmetry is 
a necessary property of hydrogen. The maximally-supcrsymmetric Af = 4 theory is an illuminating 
example to study, so we turn to it now, but we stress that the general structure of holographic duality 
thus revealed is more general, and does not depend on supersymmetry. 

* Here, engineering dimension refers to the naive dimension of an operator as determined by 
dimensional analysis on the kinetic terms in the classical Lagrangian, while the true scaling dimension 
is determined by the scale-dependence of correlation functions containing the operator in the full 
interacting quantum theory. For a review of engineering and general scaling dimensions see e.g. the 
text by Sachdev [340]. 
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Now, as pointed out by 't Hooft [341] (see also [342]), anytime you have a gauge 
theory of N x iV matrices, there is a natural way to organize the Feynman diagrams in 
terms of underlying genus-g surfaces. s Explicitly, since every field comes with two gauge 
indices, and since all indices in a gauge invariant observable must be contracted, every 
Feynman diagram can be written as a ribbon diagram in which each propagator is a 
ribbon with a gauge running along each side of the ribbon. One can then show that any 
given Feynman diagram, with any number of underlying loops, can be drawn without 
the index lines crossing only on a genus-g Riemann surface, where the specific genus 
g depends on precisely how the index lines are contracted to make the given Feynman 
diagram. Explicitly, a diagram with genus g = can be drawn on a piece of paper with 
no lines crossing, a diagram of genus g — 1 must be drawn on a torus to have no lines 
cross, etc. 

A remarkable result is that if we reorganize the loop expansion in g\ M as a double 
expansion in A" and A = gy M N, then every observable can be expressed as a power 
series expansion in A" where the powers of A" for a given diagram is determined only by 
the genus g on which that diagram can be drawn with no crossings. For example, the 
free energy J 7 = In Z takes the form 

T = N 2 fo(X) + /i(A) + ^/ 2 (A) + • • • = N 2 ' ( 9? ) 

9 

where f g (X) is the sum of all diagrams arising at genus g and is a function of only A, 
independent of N. A similar result is obtained for every possible observable, with the 
only difference being the overall factor of N, 

^ = Nm U(A) + ^i(A) + ^A(A) + ...)= i\T £ ^ (98) 

^ '9 

This double expansion implies a remarkable simplification at large A^: if A^ 3> 1, the only 
term which matters in the genus expansion is the leading planar term; any diagram that 
cannot be drawn on a piece of paper without crossings, while possibly large in numerical 
value, is dwarfed by the much larger terms coming from the planar diagrams. 

Another way of expressing this simplification at large A" is on terms of large N 
factorization. Consider a set of single-trace operators" of the form Oi = Tr(. . . ). At 
large N, correlation functions of single trace operators factorize, 

(0 1 (x 1 )0 2 (x 2 )) = (0 1 (x 1 ))(0 2 (x 2 )) + O (j^) . (99) 

This follows because the disconnected parts of the diagrams necessarily involve more 
closed index loops than connected parts, and thus additional powers of N. This does 

s Here the genus refers to the number of handles on a 2-d surface: a sphere or plane has genus zero, a 
donut has genus 1. a donut with 2 holes has genus 2, etc. 

II The meaning of single-trace operators is slightly subtle, but can roughly be understood as the set of 
operators which have local, extensive classical limits. 
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not imply that the large N theory is free, for while these correlators do factorize, 
the individual one-point functions remain non-trivial. For example, the anomalous 
dimensions of all operators are controlled by the 't Hooft coupling, A. 

We thus see that large N gauge theories have two natural coupling constants: 
jj, which controls the genus expansion and factorization; and A, which controls the 
perturbative corrections to each term in the genus expansion and the anomalous 
dimensions of large N factorized operators. Suggestively, this is precisely the structure 
of closed string perturbation theory, where every amplitude is expressed as a sum of 
contributions from various genera, with the loop-counting parameter given by the string 
coupling g s and the amplitude of a given genus determined by an auxiliary quantum 
field theory whose interactions are determined by the string tension a! via the relation 

A = g n s (A(«0 + <?Mi(«0 + giMa?) + ■■■)= 9 n s E 9 2 s 9 A g (a') . (100) 

g 

This analogy led various people to speculate that large N gauge theories should be 
captured by some theory of closed strings. This turned out to be wrong, but only just 
barely: at large N, these theories are dual to closed sting theories, but the closed strings 
live in one higher dimension! 

As an aside, while not every theory is a gauge theory of N x N matrices, it is useful 
to keep in mind this example when thinking about more general examples. What N 2 is 
really measuring is the number of degrees of freedom per unit volume, or per lattice site 
if we put the theory on a lattice. Similarly, A is measuring the strength of the dominant 
interactions when the number of degrees of freedom grows large. Many theories which 
are not simple gauge theories nonetheless have a useful notion of N and A. Exactly 
when one obtains such a controlled double expansion remains in general a key open 
question in holography. 

4-2.3. The Holographic Dictionary The basic claim of holography is that every d- 
dimensional QFT defined as above can be exactly reorganized into a (d+ l)-dimensional 
quantum theory of gravity and matter propagating in AdS^+i- The precise relationship 
between these two theories is known as the holographic dictionary, a sketch of which is 
outlined in Table (4.2.3). The rest of this subsection will develop various of the entries 
in the dictionary, including a few canonical examples of computations performed via the 
holographic dictionary. 

Operators and Fields The first entry in the holographic dictionary relates the operators 
of the QFT to the matter fields in the bulk gravity: to every local operator Oi(x) with 
quantum numbers (or charges) qi and scaling dimension Aj in our QFT is associated a 
field <&i(x, r) in AdS which carries the Lorentz structure and quantum numbers required 
to act as a coupling for Oi(x) in the Hamiltonian, Note that the dimension of the 
operator, A i; is determined by the mass of the bulk field, as we saw above in (81).' 

™ More generally, A depends on the Lorentz structure of the field $i and on the full set of its bulk 
interactions. As a simple example, the relation for a scalar operator corresponding to a scalar field in 
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the bulk is m 2 L 2 = A(A - d). 
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This means that the mass and other couplings of the bulk gravity theory generally do 
not correspond to interactions of the boundary theory in any simple way, but rather 
determine which boundary QFT we are studying. 

An Equivalence of Partition Functions Given this map, the central claim of holographic 
duality [302, 303, 304] can be succinctly expressed as 



where Zqft[^] is the partition function of the QFT as a function of the sources 
Ji for each operator Oi, while Zqg[$[Ji]] is the quantum partition function of 
the gravitational theory described in Sec. (4.2.1) computed in an AdS spacetime 
background." 1 " As discussed, the quantum gravity partition function must be evaluated 
on field configurations $j which asymptote at the boundary to the sources Jj of the 
QFT, hence the notation $[J]. This recalls the heuristic derivation above: the bulk 
fields are precisely the coupling constants of the QFT promoted to dynamical fields on 
the full RG-extended spacetime in which the RG scale becomes a physical coordinate. 

An important ingredient in this recipe is a proper definition of the boundary value 
of the bulk field. As we saw above for a bulk scalar <3> of mass m, dual to a scalar 
operator O of dimension A given by A(d — A) = m 2 L 2 , the boundary value of the bulk 
field is in general divergent. To specify the boundary condition we must rescale the 
bulk field by an appropriate wavefunction renormalization which picks off the leading 
divergence of the bulk field near the boundary, 



We will return to this renormalization when we discuss the computation of one- and 
two-point functions. 

Strong Coupling and Perturbative Gravity In general, both sides of Eq. (101) are 
complicated and computationally intractable objects. However, when one side or the 
other can be evaluated semi-classically, we can use Eq. (101) to generate a controlled 
strong-coupling expansion via the dual description. Precisely whether, and when, such 
a limit is obtained is a delicate question that must be studied in detail on a case by 
case basis. To get a sense for the general structure, let's consider the canonical example, 
the large N limit of the 4d M = 4 gauge theory with gauge group G = SU (N) and 

+ In general, Zqq means the partition function of closed string theory expanded around the AdS 
geometry. In practice we can only compute this explicitly in special cases, and then only in regimes in 
which the string coupling and curvature are sufficiently weak. Thus we may use this duality either to 
deduce the LHS of Eq. (101) when the gravity is classical and weakly curved, or to define the RHS of 
Eq. (101) when the QFT is computable. 




(101) 




(102) 
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't Hooft coupling A. The dual theory turns out to be a closed string theory on AdS 5 
with AdS-radius L. The holographic dictionary then tells us that, 



where £ p is the Planck length determining the scale at which quantum effects occur in 
the bulk, and l s is the string length controlling the scale of stringy higher-curvature 
corrections to the bulk gravitational action, as discussed in Sec. 4.2.1. The first relation 
in Eq. (103) thus tells us that the bulk description becomes classical (L S> £ p ) only in 
the large N limit,* and that the gravitational description is deeply quantum mechanical 
when N is small. In other words, the loop counting parameter in the bulk is jr. 

With an eye toward the general class of holographic theories, it is worthwhile 
unpacking Eq. (103) in more detail. In the case at hand, i.e. SU(N) SYM, N 2 
measures the number of degrees of freedom per point in the QFT. In particular, all 
extensive thermodynamic quantities (entropy, energy, etc.) must scale as N 2 . In this 
light, what Eq. (103) shows is that the quantum corrections in the bulk are negligible 
when the number of degrees of freedom per point in the dual QFT is large, and vice- 
versa. Indeed, in a host of controlled models in string theory, the appropriate version of 
the first equation in Eq. (103) takes the form, N b = (L/£ p ) d_1 , where N is a conserved 
charge and b is some real parameter [343]. Similarly, it is not always possible to interpret 
A as a t Hooft parameter in the dual QFT. More generally, the role of A is played by 
the typical scale of anomalous dimensions in the QFT [333]. We can understand this as 
follows: if the QFT is weakly-interacting, the effects of renormalization will generally 
be weak and most observables will be well-approximated by their classical cousins, i.e., 
quantum anomalous dimensions will be small; if interactions are strong, quantum effects 
will generally drive anomalous dimensions to be large. 

While N 3> 1 ensures that the Planck scale is negligible, so that the gravitational 
interactions can be treated with classical effective field theory, it does not tell us that 
the system is well approximated by pure GR. Indeed, in principle there will be a host 
of higher-curvature corrections to the classical GR action suppressed by powers of ^f. 
To ensure that these corrections are in fact negligible, i.e., that L £ s , we must 
furthermore take the 't Hooft coupling to be large, A ^> 1. Conversely, when the QFT 
is weakly coupled, A 1, the bulk geometry is sufficiently strongly curved that GR is 
swamped by higher curvature corrections to the stringy effective action. 

The main lesson here is that when the QFT is perturbative, the bulk gravity is 
out of control, and that when the QFT is strongly interacting and has a large entropy 
density, the dual gravitational description is simple semi-classical gravity. To see why 
such a weakly-curved limit is interesting, let's consider our canonical example at large 
N, N ^> 1, and with strong 't Hooft coupling, A ^> 1. Let's also analytically continue 
to imaginary time, u = iloe with ue > 0. In this limit, the gravity sector reduces 

* This provides an out from the Weinberg- Witten argument, which assumes that there arc a finite 
number of local degrees of freedom. 




(103) 
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to classical Euclidean GR. In terms of our fundamental relation, the RHS can thus be 
expressed as a sum over saddles. Focusing on the dominant saddle, we get 

Z QFT [A - e~ Jo *™, (104) 

where 

IgrMJ]] = 77TTT- [d d+1 x^(R-2A + C m m (105) 

167TGAT J 

is the classical gravitational action expanded about the dominant classical saddle, 
subject again to the condition that the normalizable modes of all matter fields are 
determined by the sources of the dual QFT. Equation (105) thus defines a simple classical 
functional of the sources from which we can compute quantum correlation functions of 
the dual QFT via the expressions above as, 



<5Ji(xi) . . .6J n (x r 



Ji=0 



(106) 



The upshot is that, when such a dual semi-classical limit exists and is reliable, we 
can compute the quantum correlation functions of our QFT by finding solutions to the 
classical field equations of the dual gravitational system, evaluating the action on-shell 
as a functional of the boundary values of the fields, and taking appropriate derivatives. 

Example: One- and Two-point Functions for a Scalar Operator As an example of this 
machinery at work, let's compute the one- and two-point functions for a boundary 
operator O dual to a free bulk scalar field <fi. For simplicity, we will work in Euclidean 
momentum space throughout. 

For the one-point function, the basic form of the computation is straightforward. 
Intuitively, our prescription (106) tells us that the one-point function corresponds to a 
variation of the classical gravitational action with respect to the boundary value of the 
bulk field. But a variation of the action with respect to the boundary value of the field 
is precisely the conjugate momentum, II, of the field $ in the direction normal to the 
boundary, 

n = - y /=jjg rr d r $. (107) 



Thus we should expect to find, with k = {oje, k 



^- SJ fm 1 ^^ mrh (108) 

This is almost correct. The trouble is that the on-shell classical action, 7gr[0]> an d 
the bulk field, $, both generically diverge as we approach the boundary, cf. Eq (81). 
To absorb these divergences we must (a) regulate all quantities by evaluating them not 
at the boundary at r = but at r = e C 1, (b) rescale the bulk field by an overall 
wavefunction normalization as in (88), and (c) add boundary counterterms to the action 
evaluated at r = e so that the renormalized on-shell action remains finite as we remove 
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the cutoff. In retrospect, this should not be surprising: the radial coordinate r is playing 
the role of a lattice cutoff, so we must work with properly renormalized quantities to 
avoid confusions. This holographic renormalization [344, 345, 346, 347, 335, 348] is 
simply the bulk realization of the UV renormalization of the boundary QFT. 

In fact, we have already seen this effect above. Recall that the relationship between 
the amplitude of the non-normalizable mode <pd-A and the boundary value of the bulk 
field $ also required a subtle wavefunction renormalization (88), so that the source is 
given by, 

J(k) = (j) d - A (k) = \imr A - d $(k,r) . (109) 

r->0 

Upon performing a similar renormalization of the bulk Euclidean action, one finds that 
the correct expression for the one-point function in terms of the renormalized radial 
momentum is" 

(OVk)) = lim r d ~ A U(k, r) . (110) 
Together with (107) and (81), this gives 

(0(k)) = ^ r ^Mk), (HI) 

where a is the coefficient of the sub-leading term in (81). Thus, just as the non- 
normalizable mode of the bulk field 4>d-A determines the source, J, for the boundary 
operator Oa, the subleading normalizable term a determines the response, (O). 

This make computing the linear-response Green functions easier. In linear response 
theory, an infinitesimal source J{x) generates a response which is linearly proportional 
to the source, with the ratio defining the linear-response Green function, 

G F (k) = MM> _ i imr 2(A-«)HM (m) 
Using our above expressions for the source, J, and response, (O), we thus find 

L (pd-A{k) 

Note that we could also have derived this result by taking two derivatives of the partition 
function. 

Computing the 1- and 2-point functions thus follows from determining 0a in terms 
of 0d-A- As we have seen, in AdS these two modes are not independent: requiring the 
bulk solution to be regular at the Euclidean horizon imposes a relation which we can 
use to determine the response 0a in terms of the source <pd-A- I n practice, finding the 
precise relation involves solving the bulk equation of motion subject to the boundary 
conditions, i.e., to solving a set of 2 nd -order elliptic partial differential equations. For 
simple examples this can be done analytically; more generally one is forced to use some 

' For a detailed discussion of this calculation sec Appendix C of [184]. The basic strategy can be found 
in [349, 350, 351, 237, 339]. 
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form of matched asymptotic expansion or numerical integration to determine 0a m 
terms of 4>d-A- 

The key point is that the computation of the quantum correlation function of the 
boundary QFT has been reduced to solving a set of classical partial differential equations 
with fixed boundary conditions. 

An example where this can be done analytically is our example of a massive scalar 
in pure AdS, for which we have in fact already determined this relation in momentum 
space, Eq. (87). Plugging this in and again working in momentum space gives, for the 
Euclidean Green function, 

' l 2 Mr(A - i ) y ' ' 

This is the exactly the form for the two point function of a scalar operator with scaling 
dimension A in a CFT. 



Real-time Response and Retarded Green Functions So far we have focused on Euclidean 
correlation functions. In fact, these holographic techniques can be readily extended to 
intrinsically Lorentzian computations needed for real-time response. 

The simplest thing to do would be to simply compute the Euclidean Green function 
and analytically continue. In practice this is often not tractable. For example, we will 
sometimes be interested in particles moving on light-like trajectories which are not easily 
described in Euclidean continuation. Instead, one can construct an intrinsically real-time 
holographic prescription, as was first proposed by Son and Starinets [352] by essentially 
analytically continuing the Euclidean prescription. Their results have since been justified 
by developing a full Holographic Schwinger-Keldysh formalism [353, 354, 355]. Here we 
will forego the formalities. We assume that such a justification can be made and proceed 
to the prescription. 

The prescription requires repeating our Euclidean prescription step by step in 
the real-time Lorentzian geometry. The key difference is that we must choose an 
appropriate boundary condition at the horizon. The appropriate choice depends on 
which Green function we wish to compute.^ Intuitively, one might expect that the 
retarded Green function, Gr, should correspond to imposing causal in-falling boundary 
conditions at the horizon, while the advanced Green function, Ga, should involve 
acausal out-going boundary conditions. This turns out to be precisely correct. Once we 
have a solution satisfying the appropriate boundary conditions, we again compute the 
properly renormalized on-shell Lorentzian action, identify source and response from the 
asymptotic behavior of the solution near the boundary, and compute the Green function 

^ Since we are only studying infinitesimal sources and responses, there is no reason to restrict to finite 
action modes, i.e., to solutions which are regular at the horizon. This is a simplification afforded by 
linear response. More generally, when considering non-trivial bulk field configurations, we should again 
impose regularity and in-falling boundary conditions. 
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via linear response. The final result takes an analogous form to Eq. (113), 




(115) 



where oj is the Lorentzian frequency and </>a, <j>d-A ar e found by solving the linearized 
Lorentzian equations in the bulk with in-falling boundary conditions at the horizon. The 
advanced Green function is computed analogously with out-going boundary conditions 
enforced at the horizon. 

Cautionary Notes on the Semi- Classical Limit While it is true that we have simplified 
our job by focusing on limits where the gravitational partition function Zqr can be 
computed by saddle point approximation, this is not the same as expanding Zqft semi- 
classically; in the present saddle, the classical description is of classical gravity in one 
higher dimension. Thus, rather than studying a semi-classical limit of the QFT governed 
by quasiparticle perturbation theory, we are focusing on a very different emergent limit 
in which quasiparticles are replaced by geometry. 

Note, too, a practical consistency condition: it cannot be the case that both Zqft 
and Zqq are simultaneously perturbative, for if that were the case we would learn that 
every classical field theory can be reorganized into classical gravity - but this is explicitly 
forbidden by the Weinberg- Witten theorem. It is thus clear that this holographic duality 
had to be a strong-weak duality, with one side becoming strongly quantum mechanical 
whenever the dual becomes semi-classical. 

Finally, we could just as easily have inverted our logic and studied a limit in which 
the QFT partition function became semi-classical. While this is unlikely to teach us 
very much about the QFT, since perturbation theory is already well understood, it may 
teach us a great deal about the strongly quantum gravitational dual. In particular, this 
approach provides the only non-perturbative definition of 4d quantum gravity currently 
known, and can be used to quickly and convincingly argue for the preservation of 
unitarity by the formation and evaporation of black holes. 

Thermodynamics Encoded by Black Holes The final entry in the holographic dictionary 
sketched in Table 4.2.3 relates the thermodynamic state and ensemble in which we place 
our QFT to the precise geometry in which we study the bulk gravity [356]. When 
the QFT is exactly conformal, the bulk geometry is pure AdS. Recalling that the UV 
physics of the boundary is associated with the asymptotic near-boundary region of 
AdS. Deforming the theory away from this conformal UV fixed point to generate a QFT 
with non-trivial RG flow corresponds to studying a geometry which is asymptotically 
AdS near the boundary but flows away from pure AdS as we run into the bulk. 
For example, turning on a finite temperature, T, and chemical potentials, /Xj, in this 
CFT corresponds to studying asymptotically-AdS black hole spacetimes with Hawking 
temperature Tjj = T and charges Qi = /jj. More generally, the thermodynamic data 
of the QFT is entirely encoded in the thermodynamics of the black hole in the dual 
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geometry (See Table 4.2.3 for the mapping between Thermodynamic parameters of the 
QFT and the quantum numbers of the dual black hole). 

The classic example of the power of these thermodynamic relations is the 
computation of the thermodynamic properties of a large N Ad Yang-Mills gauge theory 
at strong 't Hooft coupling. The holographic dual of this system is a planar AdS black 
brane, as discussed in Sec. (4.2.1). In the large N, large-A limit, the gravitational 
description reduces to GR and the black brane entropy becomes 

Sbh = ^ (116) 

where Ah is the area of the horizon. Since the system is translationally invariant, this 
diverges, so it's convenient to consider instead the entropy per unit area in the field 
theory, 

_ Sbh A h 
S ~ Aq ~ AGnAq ' [LU) 
where Aq is the area of the AdS boundary. Plugging in factors, this gives a strong- 
coupling computation of the entropy density and energy density as, 

7T 2 3VT 2 

<? - __/\r 2 T 3 e - —N 2 T 4 (1181 

■^strong 2 ' c strong 8 V / 

These strong coupling results are extremely close to the results in the free theory as 
computed in perturbation theory: 

5free = ^iV 2 T 3 , 6 frec = y iV 2 T 4 . (119) 

Remarkably, despite running from weak to strong coupling, A = to A — > oo, all that 
happened to the thermodynamics is a mild renormalization by a factor of |: 

■Sstrong _ ^strong _ 3 (120) 
"Sfree ^free 4 

Studying the full A-dependence of this ratio reveals that the energy and entropy densities 
flow smoothly and monotonically as we run from weak to strong coupling. Analogous 
results are obtained for a wide variety of gauge groups, matter contents, and even 
dimensions, with the ratio SBtIong typically within 10% or so of the Af = 4 value, §. 

•Sfree 4 

This holographic result already tells us an important fact: thermodynamic 
quantities like the free energy and the entropy density are not good probes of the strength 
of interaction of a quantum liquid. Similar results apply to other thermodynamic 
quantities, for which the strong/weak ratio is typically 0(1) as wehV 

' This can of course change if there is a quantum phase transition as parameters in the system are 
varied. 
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Holography is far more general than the refined context in which it was discovered. 
In particular, using holography as a tool does not require supersymmetry, string theory, 
d-branes, or 5-dimensional AdS, although all those ideas were central to its discovery 
and can provide useful guidance to the applied holographer. It is useful to address this 
directly. 

In its modern form, holography was first discovered in string theory while 
studying very special toy systems with as much symmetry as possible, in particular, 
maximally supersymmetric 4d A/"=4 SU(N) Yang-Mills theory at large N. The reason 
supersymmetry was important in the discovery of holography is that SUSY turns off 
many quantum effects so that a semi-classical analysis can give, for at least some special 
quantities, exact results. This allowed the computation of a host of correlation functions 
at weak coupling which could then be reliably extrapolated to strong coupling. For 
example, thermodynamic quantities like the free energy and entropy density, as well as 
transport properties such as conductivities, viscosities, susceptibilities and sound mode 
dispersion relations can all be calculated analytically, something one would never expect 
in a less constrained theory such as QCD. Of particular interest were grey-body factors 
of supersymmetric black holes, corresponding to the probability with which a particle 
would be absorbed by a black hole. These could be computed at weak coupling via 
classical gravitational perturbation theory, then reliably extrapolated to strong coupling, 
giving a precise computation of a small set of very special absorption cross-sections of 
deeply quantum mechanical black holes. Thus, in the late 90's a number of researchers 
intensively computed strong-coupling results for the string theorist's spherical cow. 

The results were astonishing. As first pointed out by Juan Maldacena in the fall of 
1997 [302], while the weakly coupled M = 4 field theory behaved much as familiar 4d 
gauge theories, at strong coupling the correlation functions of the theory reproduced, 
rather miraculously, classical scattering off black holes in one higher dimension! More 
precisely, Maldacena [302] (and, shortly thereafter, Gubser-Klebanov-Polyakov [303], 
and Witten [304]) gave a precise relationship between the 4-dimensional M = 4 theory 
and a particular theory of gravity in 5-dimensional AdS. Crucially, the 't Hooft coupling 
A and number of colors N of the 4d supersymmetric Yang-Mills theory were related to 
the AdScurvature in units of the Planck and string lengths precisely as reviewed above 
in Eq. (103). Thus, when the 4d theory is weakly coupled, A < 1, the 5d theory is 
strongly curved and cannot be treated as a theory of perturbative gravity. Conversely, 
when the 4d theory is strongly coupled, A ^> 1, and at large N, the 5d theory is weakly 
curved, so while the system cannot be treated as a weakly interacting gauge theory in 
4d, it can be treated as a weakly coupled theory of gravity in 5d. 

The point is this: while the original discovery of holographic duality required 
supersymmetry and conformal invariance, the duality itself does not depend on 
supersymmetry or conformality. Indeed, we now have a long list of examples in which 
each of these constraints is weakened or removed, with every consistency check passed. 
To date there is no example of a violation of a sharp holographic duality. From the 
stringy point of view, holographic duality is simply a true, if as yet unproven, fact 



Strongly Correlated Quantum Fluids 

about quantum field theories and quantum gravity. 



92 



4-3. Applied Holography 

Applied holography is the application of holographic dualities to construct 
computationally tractable toy models of behavior which has proven difficult to model 
with traditional tools. Key to this endeavor is the fact that various features of the physics 
which appear universal, or even trivial, on one side of the duality may be significant, or 
simply obscure, in the dual. Most current work in applied holography has used universal 
features of gravitational models, and in particular, the universality of black hole horizon 
physics, to make new discoveries about QFTs, though using universal features of weakly 
coupled field theories to learn about extreme phases of gravity is of interest as well. 
Given the aim of this Focus Issue, we will emphasize the first direction in what follows. 

The simplest application of holography involves studying QFTs for which the precise 
gravitational dual is explicitly known. In such cases we can use the perturbative QFT 
to compute any quantities of interest at weak coupling, A — > 0, then use the weakly 
coupled gravitational dual to compute the corresponding quantities at strong coupling, 
A — > oo. This strategy has very limited applications because the precise duals of most 
QFTs are not yet known. At the moment we only know the precise duality for a special 
set of non-generic quantum field theories, although they are fairly close cousins of the 
theories studied in particle physics. For example, we know the precise gravitational 
dual of the N = 4 theory discussed in detail above and have used it to learn a great 
deal about the structure and strong-coupling physics of this intricate and exotic theory 
(which we emphasize does not appear to be realized anywhere in nature). By contrast, 
we do not know the dual of the Hubbard model. 

This is not to say progress hasn't been made, just that it is difficult. An example 
of this point is the recent conjecture of a duality between the Id Ising model and 
pure Einstein gravity in AdS3. For this relation to hold the gravitational side must 
be treated not only quantum mechanically, but fully non-perturbatively, summing over 
an infinite number of topologically inequivalent locally-AdS3 spacetimes, including a 
series of quantum corrections around each classical saddle, something that can only be 
done in closed form in the absolutely simplest possible gravitational systems, namely 
pure Einstein gravity in 2 + 1 dimensions. This conjecture underscores the difficulty of 
constructing the precise holographic dual of a garden- variety QFT. 

If the precise dual is not known, a second approach to applied holography is to 
make use of a known exact duality with a class of QFTs that are structurally similar 
to the QFT of interest. Universal properties of these sibling QFTs can then be inferred 
to apply to the original QFT, too. An example of this strategy is the construction 
by Kachru et al. of holographic dimer models [357, 358]. In this construction, the 
degrees of freedom of interest for the dimer problem form a small subset of the full set 
of holographic degrees of freedom. At low energy, this subset forms a relatively generic 
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dimer model. By exploiting the holographic description, much can be derived about the 
low-energy physics, including a novel phase transition between Fermi-liquid and non- 
Fermi liquid behavior. The most powerful example of this approach so far has been the 
study of strongly-coupled quark-gluon liquids in large N toy models of QCD, for which 
a menagerie of models have been considered, each emphasizing and reproducing various 
features of the observed phenomena. As we shall see in detail in Sections (4.3.1) and 
(4.3.2), this approach has been used to motivate the anomalously small viscosity of the 
RHIC fireball, to explain why the rapid thermalization of the RHIC droplet is not so 
surprising, and to estimate the jet quenching parameter. 

However, perhaps the most interesting use of holography is to provide an entirely 
new way of defining consistent QFTs. Traditionally, QFTs in dimensions greater than 
2 are defined by specifying some Lagrangian or Hamiltonian governing the interactions 
of a set of well-defined quasiparticles. Yet many systems, including some of the 
most theoretically and experimentally interesting systems, manifestly do not have any 
well-defined quasiparticles upon which to base a standard QFT. From this point of 
view, holography looks interesting because it provides a recipe for computing quantum 
amplitudes that manifestly satisfy the conditions of a good QFT (local, causal, etc.) but 
which is defined without any appeal to a quasiparticle picture or any sort of conventional 
perturbation theory. Holography thus provides an entirely new way to construct 
consistent models of many-body physics, replacing quasiparticles with geometry as the 
central organizing principle. 

4-3.1. Viscous Hydrodynamics from Gravity Suppose we are given a QFT at finite 
temperature and density and want to know whether it is weakly or strongly interacting. 
A natural observable to query is the shear viscosity, rj, which measures the efficiency of 
momentum transport across a velocity gradient. Explicitly, the shear viscosity 77 of a 
fluid determines the frictional force F induced on plate of area A moving at a velocity 
v a distance L over a fixed surface, 



This relation follows from the hydrodynamic stress tensor given in Eq. (44). As discussed 
in Sec. 3.2 shear viscosity has dimensions hn, where n is a density. The natural 
thermodynamic quantity to which one can compare the viscosity is thus the entropy 
density, s, with dimensions of fc^n. rj/s thus gives a simple dimensionless measure of 
the efficiency of transverse momentum transport. 

^From kinetic theory, the transport of momentum proceeds by scattering 
quasiparticles between layers of the fluid, and is thus proportional to the mean free 
path between scattering events between fluid quasiparticles. Since the mean free path 
decreases as the scattering cross section increases, the shear viscosity should be inversely 
related and thus strongly sensitive to the strength of interactions. However, if we drive 
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up the strength of interaction until the mean free path becomes of order the Compton 
wavelength of the would-be quasiparticle, the kinetic-theory relation between interaction 
and viscosity should no longer apply. Meanwhile, in this post-quasiparticle setting, the 
mean free path is much smaller than any typical velocity gradient, so a hydrodynamic 
description should set in. Thus, as the interactions become strong, we expect the shear 
viscosity to become small, and furthermore expect coupling-dependence of the viscosity 
to qualitatively change, since the kinetic behavior must eventually cease. But what 
values should rj/s take, and how will it vary as the coupling grows large? A priori, the 
answers to these questions would appear to be strongly system-dependent. 

To address this question, let's compute rj/s holographically.-'' This will give us a 
controlled strong-coupling estimate of t]/s. However, before proceeding, let's pause to 
consider, from the holographic perspective, why there must be any viscosity in the first 
place. Recall that the holographic dual of a strongly-coupled relativistic large iV gauge 
theory at finite temperature is a black hole with non-zero mass. But anytime we have a 
black hole, we will have dissipation: if we scatter a wave off the black hole, some fraction 
of its amplitude, and thus of the momentum it carried, will be absorbed into the black 
hole. So the fact that there is always a generic source of dissipation is very much built 
in to holographic duality. 

A useful way to compute rj/s begins with the Kubo formula for 77, 



which relates r/ to the retarded Green function Gr for shear modes of the stress tensor, 



Note that we take oj — > only after sending k — > 0. The Kubo formula follows from 
matching linear response theory to the expected low frequency, low momentum limit of 
the correlation function in hydrodynamics; see for example [317]. To obtain rj we need 
to compute G R (t,x). In general, this is difficult in a strongly coupled QFT. However, 
holographically this calculation turns out to be relatively straightforward. According 
to the holographic dictionary, the operator T xy is dual to the corresponding shear mode 
of the bulk metric, g xy . We thus need the response of the classical gravitational action, 
/Adstfl 1 ], to variations of this mode, Sg xy = h xy . So long as we are expanding around 
pure AdS, h xy is governed by a simple massless scalar action of the form in Eq. (79), 
but with the normalization of the kinetic term determined by the normalization of the 
gravitational kinetic term, 1fi 1 r , 



— lim 



lmG R (uj,k) 



(122) 



G R (t, x) = (T xy (t, x)T xy (0, 0)> 0(t) . 



(123) 




(124) 



* Sec [359] for a summary of the current state of such computations, and [237, 339] for a more detailed 
exposition of the computation we outline, and for references to the original literature. 
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Since m = 0, we have A(d — A) = 0, so the dimension of our dual operator is A = d. 
Let's use this action and the formalism developed above to compute Gr. With A = d, 
the retarded Green function takes the form 

G R {k) = hm— — . (125) 
<S> ixi (u,k,r) 

It is readily checked (cf. 114) that the real part of the ratio above scales, after sending 
k —> 0, as u 3 , and thus vanishes in our limit. We can thus replace the taking of the 
imaginary part in Eq. (123) by 4 to get 

7? = - lim lim "frfofo . (126) 
w,fc->or-»o z W $ in (w, k,r) 

At this point something remarkable happens. In the limit that uj, k — > 0, the radial 
evolution of II and u& trivialize. This can be seen by inspection of the definition of II 
and the Hamiltonian equation of motion, d r (u§) ~ ull and <9 r II ~ k^k^d^dn^. Thus, 
in the limit u, k — > 0, both II and u& become independent of r. 

We can thus evaluate the ratio appearing in our Kubo formula at any radial 
coordinate we like. A convenient choice is the horizon, where the solution must satisfy 
in-falling boundary conditions. From our earlier discussion of the in-falling solutions 
(85), which implied that /9 r $i n = iouQ- m at the horizon, and using the definition 
nin = ~ wnG N V~99 rr d r ^m, we have that 

lim n in = -— — 3 iu$ in {r H ) ■ (127) 



The factor of r H can be nicely interpreted as the area of the horizon in boundary units, 

Ah 
Ag 



J j l . Plugging all of this into our Kubo formula, Eq. (126) then becomes 



IQtiG^Aq ' 



128) 



where Ah is the area of the horizon and Aq is the area of the boundary. From Eq. (117), 
^ = AGns. Putting these together, we see that the viscosity-to-entropy-density ratio 
takes the simple form now known as the KSS form after the authors of [11] (though see 
also e.g. [235, 236, 301]), 

5 = i (129) 
s Att 

which is orders of magnitude smaller than that of typical liquids. In fact, the liquids 
which are believed to come close to this level of perfection are the QGP at RHIC and 
fermions at unitarity, as shown in Fig. 2. 



Many features of this result are remarkable. First, as motivated heuristically above, 
rj/s is indeed a good diagnostic of strong coupling in any many-body system. Near weak 
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coupling, QFT perturbation theory gives 

7} A 



s A 2 In(VX) 



(130) 



where the constant of proportionality is system-dependent. In particular, weakly 
coupled QFTs generally have ^ > 1. Near strong coupling, by contrast, a more careful 
gravitational analysis gives 

- = 1z + 0(4e), (131) 



4vr \X 3 /\ 

almost independently of the details of the system. Thus, in any such system, the 
behavior of ^ changes dramatically as we run from weak to strong coupling. 

Second, Eq. (129) is remarkably generic in holographic systems. Indeed, in any 
QFT which is holographically dual to GR, whether the gauge group is SO(N) or 
SP(2N) or something more complicated, whether enjoying maximal supersymmetry 
or no supersymmetry, we obtain precisely - = j-. To see why this result is so general, 
note that the radial evolution dropped out of our computation of rj, so that all that we 
needed was the local physics near the horizon. The rest of the holographic description 
(including the UV physics near the boundary) was moot. Meanwhile, s is explicitly 
proportional to the horizon area, so again is completely determined by the physics of 
the horizon. The generality of this holographic result then follows from the exceptional 
universality of the physics of horizons in GR. 

Observations such as these led Kovtun, Son and Starinets (KSS) to make the 
conjecture that r\ > j- is a universal lower bound on the shear viscosity of any quantum 
many-body system, with saturation occurring only for those strongly coupled theories 
which are holographically dual to pure Einstein gravity. We now understand that this 
result is not universal. Explicitly, Eq. (129) holds when the spacetime geometry is 
AdS, the gravitational action is Einstein-Hilbert, and the fluid is time- independent, 
homogeneous and isotropic. However, when higher-derivative terms are present in the 
gravitational action (corresponding to finite- N and finite- A effects in the dual QFT), 
or when the geometry is deformed away from AdS such that the radial evolution does 
not factor out (as occurs for example when the geometry is time-dependent, spatially 
disordered or sufficiently anisotropic), - can depart from the KSS value, and indeed can 
be lower than j-. 

The cleanest example of such an effect involves [360, 10, 361] adding the leading 
irrelevant operator to the bulk gravitational action, the 5d Gauss-Bonnet term, a 
quadratic scalar built out of the curvature tensors, 

Igb = y^g^ j ^ ^ (~ 2A + R + XgbR2gb) ■ (132) 

Such terms arise in known string theories where we know and have control over both 
sides of the duality. E As it turns out, the holographic analysis including Xqb ^ is 

s In the QFT, the Gauss-Bonnet coupling \qb is related to the Euler and Weyl anomaly coefficients. 
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straightforward. All that materially changes is the normalization of the action for the 
scalar mode h xy , which picks up an extra factor of (1 — 4Xgb), yielding 

77 _ 1 - 4X GB 
s 4ir 

When Xqb > 0, the KSS conjecture is violated." 

It is worth understanding why this happens. The value of t]/s is controlled by the 
normalization of the kinetic term in the scalar action for our tensor mode h xy . Typically, 
such coefficients are best understood as coupling constants, since scaling them away to 
make the kinetic term canonical shifts the coefficient into the interactions.^ When 
we make the coupling stronger, as in the case of Gauss-Bonnet with \ GB > 0, the 
viscosity goes down. This is reminiscent of the dependence of viscosity on coupling 
in kinetic theory: increasing the strength of quasiparticle scattering inhibits transport 
and decreases the shear viscosity. At strong coupling, there are no longer any well- 
defined quasiparticles to which to apply kinetic theory; instead, the effective degrees 
of freedom mediating low-energy momentum transport are weakly coupled gravitons in 
the holographic spacetime. 

Interestingly, while this evidence rules out the KSS conjecture, there exist 
compelling arguments for a weakened version of the KSS bound. Physically, the 
arguments for such a refined bound involve demonstrating that violations of these 
bounds would lead, in a wide class of theories, either to violations of causality in the 
field theory [10, 362] or to violations of the positivity of energy [363]; in fact these 
two pathologies are intimately related. In the case of Gauss-Bonnet gravity which is 
holographically dual to 4d QFTs with M = 1 supersymmetry, both arguments lead to 
the same constraint, ^ < Xqb < 155, leading to the refined lower bound, 

77 > 1 16 
~s ~ 4^25 

In more general models, similar bounds obtain, but the precise values differ. The 
meaning and reliability of such refined bounds remains a matter of active debate. 
For example, it has been shown [364, 365] that sufficiently contrived examples may 
be concocted which violate any such bound and in fact push 2 arbitrarily close to zero, 
and even negative! However, it is not at all clear that these models are themselves well 
defined; indeed, one can immediately show [366] that all such models with 2 negative 

While it cannot be tuned arbitrarily in these examples while remaining within the regime of validity 
of the classical gravitational description, it can certainly be taken away from zero in a controlled 
fashion [361]. 

II Note that Xqb < \ is a consistency condition for GB gravity. When Xqb — j, the kinetic term for 
the graviton degenerates and the theory becomes ill-defined. 

^ The universality of 77/s in GR thus physically derives from the equivalence principle: in Einstein 
gravity, the metric couples universally with a single overall coefficient to all forms of stress-energy. 
When we muck with this universality, for example by breaking symmetries or by modifying the way 
the metric interacts with itself via higher-curvature corrections, we make our chosen tensor mode of 
the metric couple more, or less, strongly. 



(133) 



(134) 
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or vanishing are pathological, with causality and self-consistency enforcing a minimal 
value for 2 . Meanwhile, whether any such bound should be expected when we give up 
time-independence, homogeneity or isotropy is a topic of much current research. To 
summarize, there are well-defined lower bounds on - in various classes of theories, but 

' s ' 

each such bound has been violated by considering a more general class of models. It 
remains possible that there is a simple, completely universal lower bound on but the 
evidence for such a bound is increasingly tenuous. 

Worrying about universality of a particular bound, however, misses the point. The 
lasting import of Eq (129) is that it establishes a link between ultra-low shear viscosity 
and strongly interacting many body systems. Said differently, if one finds a fluid with 
an extremely small shear viscosity, one has very strong evidence that it enjoys strong 
quantum correlations. Conversely, extreme quantum fluids should be expected to have 
ultra- low viscosity of order j-. Whatever the final story about a hard lower limit, this 
connection is an important lesson from holography. 



4-3.2. Diffusion, Jet Quenching and Dynamics in the Holographic Quark-Gluon Plasma 
A related transport coefficient which can be fruitfully studied holographically is the 
diffusion constant for a heavy quark in a strongly-coupled holographic QGP. It is an 
interesting example for us for several reasons. First, the structure of the holographic 
computation is very different from the computations described above, and in fact owes 
a great deal to the string-theoretic origin of holographic duality. Second, the result is 
very different from the simple universality of the viscosity computation. Specifically, 
the diffusion constant for a heavy quark scales as 

D ~ ^7% (135) 

and is thus very strongly sensitive to the 't Hooft coupling, A. In particular, this suggests 
that even very heavy quarks should rapidly equilibrate in the QGP, in sharp contrast 
to perturbative estimates. Note that this also differs sharply from the perturbative YM 
result + 

6 / 1 



c =2^U;l • < 136) 



for while Eqs. (135) and (136) look similar, the latter is only valid for a s <l, while the 
former is only valid for A ^> 1. 

The basic structure of the holographic computation is described in [289, 291]; 
see [286] for a perturbative analysis. Suppose we place a heavy quark in the QGP and 
treat it in a Born-Oppenheimer approximation as a slow, heavy, pointlike variable which 
sources the gauge field. As we run the RG, this charge is conserved. Holographically, 
this means that we should have a point-like defect not just on the UV boundary, but 

+ This result corresponds to a specific numerical estimate of the logarithmic terms in Eq. (58); see 
Eq. (6.1) of [286]. 



Strongly Correlated Quantum Fluids 



99 



also on every slice of the holographic bulk. The heavy quark on the boundary is thus 
the endpoint of a string in the bulk which hangs down from the boundary into the AdS 
horizon, hanging along a straight line according to the gravitational tug of AdS. The 
string can be thought of as the holographic image of the color flux of the fundamental 
quark in the 3 + 1 boundary. Explicitly, the dynamics of the string are determined by 
a simple least action principle, where the action is the area of the surface in space-time 
swept out by the string in AdS. For a static string, the extremal string is simply the 
string hanging straight down from the boundary quark. 

To measure the diffusion constant, we need to give our quark a kick. The easiest 
way to do this is to drag the quark through the QGP with a fixed velocity and then 
determine the drag force. In equilibrium, we must have 

f = -VDP, (137) 

where rjjj is the drag coefficient. Holographically, the string hanging off the quark into 
the bulk will now dangle behind the moving quark, imparting a drag force to the quark. 
Using the action principle we just defined, the drag force is computed to be 

f = -l^XT 2 — . (138) 

For heavy quarks, for which the thermal corrections to the kinetic mass are negligible, 
this gives 

Vd = ^V\T 2 . (139) 
2m 

Using this in the relation D = T/rjom derived in Eq. (57) yields Eq. (135). 

It is worth commenting on how this computation is similar to the viscosity 
calculation, and how it differs. Perhaps the most important similarity is that in both 
cases we computed deeply quantum-mechanical transport quantities in our strongly- 
coupled QGP by solving simple, classical differential equations in the dual gravitational 
system. This is a basic hallmark of holographic duality. Secondly, the resulting 
computation was surprisingly independent of the details of the holographic setup, 
depending only on the gravitational dynamics of a string hanging into AdS from a 
moving point on the boundary. This suggests that a heavy fundamental in any strongly 
coupled QGP should have similar qualitative behavior, namely, efficient diffusion and 
rapid equilibration. By contrast, the final results differ in significant ways. For example, 
while r)/s was remarkably insensitive to the coupling, D depends sensitively, and is thus 
a sensitive probe of precisely how strongly coupled the theory is. 

A similar computation leads to a holographic estimate of the jet-quenching 
parameter, q, introduced in Sec. 3.6. Holographically, the role of the QGP fireball 
is played by the dual black hole, while the hard jet can be modeled as a heavy quark 
at very high energy, and thus moving on an effectively light-like trajectory. Our job is 
thus to compute the transverse momentum gained by a light-like quark moving through 
the QGP. This can be reduced [367, 368, 369] to computing a light-like Wilson loop 
in the QGP. Holographically, the expectation value of this Wilson loop can again be 
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computed by studying the surface swept out by the trailing string, determined as before 
by extremizing the string action. Since the Wilson loop is light-like, the computation 
must be performed in Lorentzian signature. The resulting real-time holographic analysis, 
while somewhat intricate, leads to a simple prediction for the jet-quenching parameter 
in W=4 SYM [369, 370], 

<k=4 = v/AT 3 . (140) 

Recalling that A = g 2 N c and g 2 = 47ra s , and plugging in T ~ 300MeV^ and a s = |, 
which are reasonable estimates for the values at RHIC, this gives ~ 4.5GeV 2 /fm, 
which is comparable to results from RHIC; see Sec. 3.6. Note that this holographic 
result leads to a semi-strongly coupled treatment of energy loss, treating transverse 
momentum diffusion in strong coupling via a holographic estimate of q, while the actual 
energy loss (longitudinal drag) is done in perturbation theory, as discussed in Sec. 3.6. 
For fully strongly-coupled treatments see e.g. [371, 372]. 

The important lesson of this holographic computation however is not the specific 
value of q we find but rather how q depends on the physical parameters of the system. 
In particular, it turns out that in any holographic QGP dual to pure AdS, and thus 
governed by a dual CFT, one finds [370] 



/ScFT (141) 



qu=A V S Af=4 

where s is the entropy density of the QGP. The energy lost by a hard parton plowing 
through the QGP is thus not proportional to the entropy density, and so cannot be well 
modeled by scattering off a gas of persistent quasiparticles. 

This holographic approach to studying the physics of a generic SU(N) QGP to glean 
insight into the dynamics of e.g. the physical QGP at RHIC has been followed in many 
directions, including several of the papers in this Focus Issue. For example, Hubeny [373] 
uses a simple but remarkably illuminating gravitational model for an accelerated heavy 
quark in large N QCD to study the collimation of synchrotron radiation emitted by 
a circling quark. Like the diffusion and jet-quenching computations sketched above, 
Hubeny's calculation involves tracing the path of the bulk string which trails behind 
the accelerating quark into the bulk of AdS. The result of the acceleration is that this 
bulk string becomes coiled, leading to the radiation emitted by the string remaining 
tightly collimated. Curiously, this result at first appears to be in tension with standard 
holographic intuition. Hubeny resolves this apparent conflict through a careful study of 
the gravitational backreaction induced by the coiled string. 

4-3.3. Gravitational Engineering and Holographic Superconductors An interesting 
example of defining the field theory by its holographic dual is the holographic 
superconductor. The idea, first outlined by Gubser [374] and developed in detail by 
Hartnoll, Herzog and Horowitz [375, 376], goes roughly as follows. Suppose we have 
a QFT which enters a superconducting phase at low temperatures by developing a 
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non-zero condensate for some scalar operator, O, spontaneously breaking a global U(l) 
symmetry. What would that look like in the holographic context? By appeal to the 
holographic dictionary in Table 4.2.3, the global current in the boundary is dual to a 
gauge field in the bulk, while the charged scalar order parameter in the boundary is dual 
to a scalar field in the bulk which is charged under the bulk gauge field. Spontaneously 
breaking (or Higgsing) the global U(l) symmetry on the boundary by giving a vacuum 
expectation value to O is then dual to spontaneously breaking the bulk gauge U(l) via 
a non-trivial profile for the bulk scalar field, $(r). 

Our goal is thus to build a gravitational system which is dual to some QFT which 
exhibits this symmetry breaking pattern - not a specific QFT, but any QFT with this 
structure. By adjusting the gravitational theory we will then be able to explore a whole 
space of strongly-coupled QFTs with low-energy superfluid phases. Procedurally, we will 
take the gravitational system as a definition of the dual QFT, computing the correlation 
functions etc of the QFT entirely by appeal to the holographic dual. 

The holographic dictionary thus tells us that, whatever else exists in the 
gravitational dual of a superconducting QFT, we will need a vector field dual to 
the boundary current J m and a scalar field $ dual to the scalar order parameter O. 
The dictionary, however, does not tell us the Lagrangian, so we must choose it. A 
minimal first guess is, 



1 1 



-F 2 + \D<$>\ 2 + m 2 \<S>\ 
4 11 11 



(142) 



We can further simplify our lives by considering the probe limit, g — > oo, in which 
case the matter field stress tensor (which scales as -4-) is negligible so that the Einstein 
equation, and thus the metric, is not modified from its original AdS form. 

Our job is now to find non-trivial solutions of the bulk equations of motion which 
follow from our trial action expanded around an AdS black brane geometry, so as to put 
the system at finite temperature. Non-trivial means such that the boundary operator 
O acquiring a spontaneous vacuum expectation value; from our dictionary, this requires 
$ to have non-trivial profile, Higgsing the vector, as expected. Once we find such a 
solution, we can compute correlation functions in the resulting system by evaluating the 
bulk action on-shell as a functional of the non-normalizable modes and taking derivatives 
with respect to those modes, as discussed in Sec. 4.2.3. 

For simplicity, let's focus on homogeneous configurations of the form A = A t (r)dt 
and <3> = <3>(r). The bulk equations of motion in our AdS black brane take the form 

r 2 fA'{-2\$\ 2 A t = 0, r 2 f&' + (rf'-2f)r&+l—±-m 2 \<l> = 0. (143) 

Near the boundary, the general solution to the bulk field equations takes the form, 

A t ~ n + pr + ... , $ ~ <p d -Ar d ~ A + <pAr A + ... , (144) 

where m 2 L 2 = A(A — d). A convenient choice for the mass is m 2 L 2 = —2, which 
gives A = 1 or A = 2. Let's focus on A = 2. The holographic dictionary tells us to 
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Figure 25. Condensates and conductivity of a holographic superconductor. Results 
adapted from code by C.Herzog [377] presented the original work in [375]. 

interpret n as the source for the J t component of the current, so n plays the role of the 
chemical potential for the boundary charge density, p, being the subleading piece of A t 
near the boundary, then represents the expectation value of the operator conjugate to 
p, p = {Jt), so we can identify p as the charge density induced by p. Similarly, <pd-A, as 
the leading term in $ near the boundary, is interpreted as the source for the boundary 
operator O, and a as the vacuum expectation value, a = (O). 

To construct a spontaneous condensate, we should turn off the source, 4>d-A = 0, 
fix a reference value for the chemical potential, p = p Q , and look for a solution of the 
bulk equations of motion which have, below some critical temperature T c , a non-zero 
value for (O) = 0a- This requires that we solve Eqs. (143) subject to the boundary 
conditions that $(0) = (no source) and A t (0) = p Q . Generally this must be done 
numerically, as Eqs. (143) are nonlinearly coupled. By evaluating the on-shell action 
and using the holographic dictionary in Table 4.2.3, we can compute physical quantities. 
For example, we find the critical temperature by searching for the minimum T of the 
background black brane such that a non-trivial solution with these boundary conditions 
exists. We find the AC conductivities by looking at linear response for fluctuations of A x 
linearized around a non-trivial solution for A t and $, etc. The precise behavior we find 
depends on the details of the gravitational model chosen. By adding additional fields 
and interactions, we can use these techniques to engineer a wide range of dynamical 
phenomena in the dual QFT. 

The point to emphasize here is that while we are explicitly computing correlation 
functions in a strongly-coupled QFT, and indeed studying the detailed physics of 
transport, the QFT under investigation is defined purely through its gravitational dual. 
We thus call these superconductors holographic superconductors, and more generally 
refer to such QFTs as holographic QFTs. 

Importantly, these models can be used to study much more than just the 
translationally invariant ground states of these strongly interacting superfluids. For 
example, they provide a simple and powerful framework for studying spatially ordered 
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phases and solitons. For example, the paper by Keranen et al. in this Focus Issue [164] 
uses solitons in these holographic superconductors as precision probes of the superfluid, 
including in particular the excitation spectrum [378, 379] around a dark soliton at 
unitarity in the BCS-BEC crossover [162, 165, 166, 167]. As the solution of this problem 
reduces to a detailed study of coupled PDEs, considerable progress can be made with 
common computational techniques towards a concrete prediction observable in present 
experiments on ultracold quantum gases. 

Notably, this holographic mechanism for spontaneously breaking a symmetry can 
be applied to any boundary QFT with symmetry breaking. For example, the paper by 
Basu et al. in this Focus Issue [380] uses the holographic superconductor paradigm to 
study color superconductivity in holographic models of the QGP. As good effective field 
theorists, Basu et al begin by writing down a gravitational effective theory with the 
minimal ingredients to mock up a color superconductor, then examine what constraints 
they must impose on the parameters of this effective description to reproduce the 
expected physics of a true color superconductor. 

4-3.4- Non-Relativistic Holography and Cold Atoms Strongly correlated quantum 
liquids of cold atoms provide an excellent target for applied holography. From a 
holographic point of view, cold atoms at unitarity are very similar to J\f — 4 models of the 
RHIC fireball: they provide examples of strongly-correlated liquids at finite temperature 
and chemical potential that are governed by a conformal field theory. The key difference 
is that these systems are non-relativistic, with dispersion relations scaling like* 

w ~ k z , (145) 

so that the RG fixed points which govern their dynamics are non-relativistic conformal 
field theories (NRCFTs) [381, 382, 383]. Furthermore, in the case z = 2 arising in cold 
fermions at unitarity, the symmetry algebra is enlarged [383] and includes a central 
Number operator, N, such that every operator has, in addition to a dimension, A, a 
conserved particle number eigenvalue, N. 

Holographically, the non-relativistic scaling has a dramatic effect. In the case of 
the M = 4 theory, relativistic conformal invariance was enough to determine the dual 
space-time, with the conformal symmetry group of the QFT, 5*0(4,2), mapping to the 
isometry group of the dual space-time; this fixes the geometry dual to the J\f = 4 theory 
to be AdSs. More generally, the fact that relativistic QFTs flow to Lorentz-invariant 
RG fixed-points in the UV implies that their holographic duals should asymptotically 
approach pure AdS near the boundary. If our QFT does not possess Lorentz invariance, 
the dual space-time cannot be simple AdS. To build the holographic dual of an NRCFT, 
then, we must start from the beginning and identify the right space-time. 

* The scaling exponent z is called the dynamical exponent, and can be extracted either from the 
dispersion relation for low-lying modes or from the scaling of the characteristic timescale with the 
correlation length, r ~ £ z , as we approach a critical point [340]. 
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The simplest geometry with these properties is the Lifschitz metric [308] 




(146) 



In this geometry, every slice at fixed r is again a copy of flat space, but the scaling of 
the time-like coordinate and the spatial coordinates under a shift in r is inhomogenous, 
with scaling uj ~ k z . At z — 1, this reduces precisely to the AdS case, a good check. 
However, at z = 2 there is no enhancement of the isometry group of the geometry, so 
it's natural to look for another geometry which realizes the full enhanced z = 2 non- 
relativistic conformal symmetry group. Following this logic leads to the Schrddinger 
metric [306, 307], 



where £ is a new periodic variable. Again we find that scaling implies u ~ k z . Now, 
however, the isometry group of the manifold is in fact enhanced at z — 2. Moreover, 
since £ is compact, we can expand every field in eigenmodes of N = i<9g, so that bulk 
fields, and thus all boundary operators to which they are dual, are labeled by both a 
dimension, A, and also a ^-momentum, N. We thus have a geometry which enjoys both 
of the peculiar features of fermions at unitarity identified above. 

Repeating the full construction in Sec. 4.2 in this non-relativistic setting remains 
an open problem. However, considerable progress has been made in the context of 
both Lifschitz and Schrddinger systems, including the construction of charged black 
holes corresponding to NRCFTs at finite temperature and chemical potential [384, 385, 
386, 387] and the construction of some toy non-relativistic superfluids (see e.g. [388] 
in this Focus Issue). Many of the basic universal results from the relativistic setting 
have been replicated. For example, in the simplest models we again find - = Note 
that, while we have suppressed factors of c throughout, there are no such factors in 
this result for the viscosity ratio so we get a non-trivial prediction for non-relativistic 
systems. Furthermore, while - is a natural observable in relativistic QFTs even if we 
don't know about holography (it is the dimensionless momentum diffusion constant in 
linearized hydrodynamics), the fact that ^ is relevant in non-relativistic systems is an 
interesting feature of holography; the momentum diffusion constant in non-relativistic 
hydrodynamics is -. In the context of Lifschitz scaling, considerably more progress has 
been made, including interesting results for models of strange metals [389]." 

" It is worth emphasizing that it is straightforward to build holographic models in which the dynamical 
exponent changes as the theory flows from the UV to IR. A particularly illuminating example involves 
the fate of charged black holes with large ground-state degeneracies. In at least some holographic 
models it has been shown that the ground state degeneracy is lifted by a dynamical instability which 
changes the near-horizon geometry from AdS to Lifschitz geometries which have no, or at least vastly 
diminished, ground state degeneracy. Since the near-horizon region encodes the IR physics of the 
dual QFT, such holographic models describe QFTs where the dynamical exponent of the low-energy 
degrees of freedom arises due to strong interactions amongst constituents with no simple quasiparticlc 
description, and is thus effectively invisible in the UV regime. The dynamical scaling is in this sense 
emergent. 
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For models of cold atoms, on the other hand, our non-relativistic holographic models 
remain rather distant from the systems studied in actual experiments. Indeed, it seems 
likely that fermions at unitarity, like QCD, do not have weakly coupled holographic 
duals. However, unlike QCD, where we have already built holographic toy models 
which capture a great deal of the strong-coupling phenomena, non-relativistic CFTs 
with dynamical exponent z = 2 and Schrodinger symmetry have proven difficult to 
model holographically. In particular, the thermodynamics and phase structure of current 
holographic models are not obviously related to any known systems, if not outright 
pathological [384, 385, 386]. For example, we expect the finite-density ground state 
of these theories to form a superfluid which spontaneously breaks the U(l) symmetry 
generated by the number operator, N. Holographically, N is mapped to d^, the generator 
of translations in the 2 nd holographic direction. Such a superfluid ground state must 
thus break translation invariance in the ^-direction. Unfortunately, all of our current 
holographic models of z = 2 Schrodinger CFTs feature manifest unbroken ^-translation 
invariance. 

This problem can be traced [390, 391, 392] to the fact that all our current models 
were derived from highly symmetric stringy systems which enforce this symmetry. As 
a result, these limitations do not appear to be intrinsic to the holographic framework 
but rather to the very specific models which have so far been studied. It would be 
of great interest to develop more realistic models (see eg [392] for recent work in this 
direction). A step in this direction is described in the paper of Adams and Wang in this 
Focus Issue [388], in which a superconducting ground state for such a non-relativistic 
holographic CFT is constructed by extending the "holographic superconductor" strategy 
from AdS to these non-relativistic geometries. This leads to various interesting effects, 
including a surprising quantum phase transition as the density of the fluid is varied and 
a multicritical point where the transition to the superfluid state switches from 2 nd to I s * 
order. This does not solve the underlying problem, as the geometry remains unmodified, 
but is at least a proof of principle that such holographic ground states do exist. 

4-3.5. Holographic Non- Fermi Liquids and Critical Phenomena Fermi surfaces are 
familiar beasts, and are generally well- described by standard quasiparticle perturbation 
theory. This is true despite the fact that the electrons living on the Fermi surface 
may interact through the infinite-range Coulomb potential. The key point is that the 
phase space available for two electrons to scatter is extremely limited, as both the initial 
electron states and the final electron states must all lie on the Fermi surface. Folding 
in this phase space constraint, the effective interaction for (dressed) electrons scattering 
on the spherical Fermi surface mediated by the electromagnetic interaction is small, and 
hence the applicability of perturbation theory. The corresponding state is called a fermi 
liquid.^ 

tt We discuss the electron gas as an example of a Fermi liquid, but we emphasize that the arguments 
given here are only true for purely Coulombic forces. Current-current interactions, for example, are 
unscreened and destroy quasiparticle behavior [393] . 
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Experimentally, many systems with a fermi surface are manifestly not well 
described by quasiparticle perturbation theory, and are thus referred to as non- 
fermi liquids. Perhaps the most famous such examples are High-T c superconductors, 
whose spectral function (as measured e.g. by ARPES experiments) do not contain 
any sharp quasiparticle peaks. While it is relatively straightforward to write down 
phenomenological models of such non-fermi-liquids, it has proven somewhat non- 
trivial to build controlled microscopic models which generate non-fermi-liquid behavior 
dynamically, and famously intractable to derive such non-fermi-liquid behavior from 
realistic models of the crystal lattices which lead, experimentally, to such behavior. 

Since holographic QFTs neither require nor generally admit any quasiparticle 
description, it is natural to ask what happens when we put a holographic QFT at 
finite fermion number density [394, 395, 396, 397]. The simplest way to do so is study a 
QFT with a fermion carrying a conserved charge under some global U(l) for which we 
turn on a chemical potential. Holographically, this corresponds to studying a fermionic 
field in AdS in the presence of a charged black hole. The charge of the black hole 
induces a chemical potential in the QFT, while the boundary value of the bulk fermion 
specifies a source for the fermion in the QFT with which to probe the spectral function. 
To compute the spectral function, which is equal to the imaginary part of the retarded 
Green function for the fermion in the QFT, we turn on a known source and measure the 
response. Holographically, this maps to solving the Dirac equation in the background of 
a charged black hole in AdS subject to a Dirichlet boundary condition at the boundary 
and infalling boundary conditions at the horizon. This turns out to be a straightforward 
numerical computation. The resulting spectral functions show clear evidence of a sharp 
Fermi surface, but, as expected, enjoy no quasiparticle poles, and so the holographic 
QFT describes a non-fermi-liquid. 

Much work has been done to clarify the physics of these holographic fermi liquids. 
For example, it has been shown that incorporating gravitational effects can change the 
results dramatically at low energy, with the black hole in some cases vomiting forth 
its charge and disappearing, leaving behind an electron star which is extremely dense 
but completely smooth, much like a neutron star but carrying a net charge [398, 399]. 
Meanwhile, completely different models of fermi-surfaces in holographic QFTs have been 
concocted [389]. Despite using very different ingredients, and having in some ways very 
different physics, all such models share the important property of having no simple 
quasiparticle poles in their spectral functions. 

One of the key lessons of this study has been the role of emergent conformal 
symmetry in the near-horizon region [397, 400]. Recall that the near-horizon region 
corresponds to the deep IR of the dual QFT. Saying that this region has an emergent 
conformal symmetry then implies that the QFT has, at very low energy, an emergent 
scaling symmetry. Importantly, the dynamical exponent of this conformal symmetry is 
generally larger than 1, and indeed often infinite, corresponding to ultra-local physics 
with a finite ground state degeneracy. As usual, this enormous ground-state degeneracy 
is a marker of instability, and these systems are quick to decay into a host of other 
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phases. As such, these holographical critical points correspond to quantum critical 
points of the dual QFT. 

All of these models of holographic Fermi surfaces and holographic critical 
phenomena share a generic two-component structure: first, a set of degrees of freedom, 
X, which correspond to the naive fermions living on the Fermi surface; second, operators 
ip in an emergent IR CFT which couple to the x fields. Integrating over the V'-fields 
then generates new, strong interactions amongst the original x fields, which is the 
origin of the strong interactions on the Fermi surface in these holographic models. The 
natural question then becomes, if all we care about is the low energy physics, which 
is entirely encoded in the near-horizon region, why bother with the entire holographic 
framework? Why not just build a truncated model which only includes the fermions 
of interest coupled to the emergent CFT which controls the interactions? Doing so 
has come to be known as semi-holography [309, 311], and has proven quite useful (see 
e.g. [401, 400, 402]). 

An example of this semi-holographic approach is the paper by Nickel and Son in 
this Focus Issue, [311]. In this paper, the authors argue that the transport of conserved 
charges is governed by an effective theory for a set of goldstone modes coupled to an 
emergent near-horizon CFT. The fact that the physical modes are goldstone modes 
powerfully constrains the dynamics and enforces much of the structure of holographic 
transport. The origin of the Goldstone bosons is quite beautiful. From the holographic 
dictionary, we know that the conserved {7(1) of the boundary QFT maps into a gauged 
£7(1) in the bulk, and thus to a dynamical emergent vector field in the bulk. Importantly, 
this vector has two boundary values, one at the boundary of AdS and one as we approach 
the horizon. Now imagine we run the RG, integrating out the high energy modes 
(which live near the boundary) shell by shell. If we integrate out the entire bulk, we 
are eventually left with a thin shell just outside the horizon, on of the boundaries of 
which we have a dynamical {7(1) gauge field, and thus a total symmetry group which 
is {7(1) x i7(l). However, this is a fake - there is only one {7(1) in the system - 
so integrating out the degrees of freedom in the bulk between these two boundaries 
must spontaneously break {7(1) x {7(1) to a single {7(1). As Nickel and Son show, 
this indeed happens, leaving behind a goldstone mode which remains coupled to the 
surviving emergent CFT. A Wilsonian analysis of the most general action allowed by 
the symmetries of the system then efficiently reproduces many classic results about 
transport in holographic models. 

These sorts of holographic approaches to quantum critical phenomena have also 
been applied to myriad problems in the physics of condensed matter broadly construed. 
For example, the paper by Bayntun et al. in this focus issue [403] studies plateau 
transitions in quantum- Hall- like holographic models. Such phenomenological models 
are becoming increasingly sophisticated and, as they improve, increasingly compelling. 
In this case, the authors construct a model for which the scaling exponent for plateau 
transitions is |, in rough agreement with experimental values of ~ 0.42. 
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4-3.6. Far- from- Equilibrium Physics in Strongly Coupled Quantum Field Theories 
Perhaps the most exciting application of holographic duality is to the study of far-from- 
equilibrium physics in strongly-interacting quantum many body systems. The basic 
idea here is extremely simple. As we have seen, the equilibrium thermal physics of a 
holographic QFT is encoded in the thermodynamics of a stationary''' black hole horizon 
in the dual spacetime. To study non-equilibrium physics in the QFT, then, we should 
study non-equilibrium physics in the dual gravity - ie, time-dependent spacetimes whose 
black hole horizons evolve in time. 

It is easy to visualize this connection. Imagine the collision of two nuclei at RHIC 
or LHC. At some initial time we have two approximately non-thermal objects in a head- 
on collision. At some special moment they collide. Quite rapidly - in fact, amazingly 
rapidly according to experimental results - the ensuing fireball thermalizes, then slowly 
evolves according to near-equilibrium hydrodynamics. Now what would this process look 
like in a holographic dual? Since the initial temperature is zero, the initial geometry 
has no black hole horizon. Meanwhile, the initial non-thermal nuclei are basically just 
moving sources of stress-energy. According to the holographic dictionary, the source 
for stress-energy in the boundary QFT is equal to the boundary value of the bulk 
metric. Our moving packets of mass are thus represented holographically by traveling 
gravitational waves which extend from the boundary into the deep interior of AdS. The 
collision of the nuclei then corresponds to a collision of the dual gravitational waves. 
But we know what happens when you collide gravitational waves at sufficiently high 
energy - they form a black hole. 

However, since the waves are not perfectly thin and the process of black-hole 
formation is not perfectly efficient, the details of black-hole formation can be quite 
messy. In particular, the horizon does not form instantaneously but grows and expands 
rapidly during the collision - this black hole is far from equilibrium. Meanwhile, while 
the initial collision, and thus the initial horizon, will generically be inhomogeneous, 
anisotropic and just generally messy, it is a theorem of black hole mechanics that the 
horizon always settles down to a completely uniform spherical surface with a uniform 
radius and a uniform temperature - and indeed that it does so relatively rapidly, with 
fluctuations of the horizon damped by the viscosity we studied above. Thus our collision 
of nuclei and thermalization of the QGP fireball at RHIC look, in the holographic dual, 
like a collision of gravitational waves which form an initially far-from-equilibrium black 
hole which then rapidly rings down to a simple uniform sphere. 

Note the simplification afforded by this dual description: while the QFT description 
involves real-time quantum dynamics in a strongly-coupled many body system, which 
are extremely difficult to simulate, the holographic model involves solving a set of PDEs 
- the Einstein equations, together with a well-defined set of initial data and boundary 
conditions - which are considerably easier to simulate. 

This dramatic simplification was exploited by Chesler and Yaffe [312], who 

t In GR, stationary means, roughly, time-independent, while static means, again roughly, non-rotating. 
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simulated the formation of a holographic QGP by numerically colliding gravitational 
shock waves to form a black hole in AdS. An important lesson of the resulting numerical 
simulations is that these holographic QGPs thermalize exceedingly rapidly. Rapid can 
be quantified in several ways. One simple measure is that the thermalization time 
is of order the time required for light to cross the horizon of the final equilibrium 
black hole. Alternatively, we can take the numerical results and rescale variables 
to reproduce the real scales of QCD. Doing so [312] gives a very rough prediction 
of the RHIC thermalization time of 0.35 fm/c, which is safely below the upper 
bound estimated via hydrodynamic simulations from the RHIC data; see Sec. 5.2. 
Another important lesson of this and related computations is that hydrodynamics 
become an excellent approximation very rapidly, with the pressure and energy tracking 
hydrodynamic predictions as fast as, and sometimes apparently even faster than allowed 
by thermalization. 

Importantly, the holographic models studied thus far do not do nearly as good 
a job in other regards. For example, the most simple models fail to reproduce the 
observed energy dependence of the total multiplicities. This is perhaps not surprising 
given the simplicity of the models studied to this end, but is an important point to 
rectify. Further, while these models give us useful insight into the phenomena, they 
do not give us any guidance on how to understand the relation between such rapid 
thermalization and asymptotic freedom or deep inelastic scattering. Nonetheless, these 
are remarkable results. 

5. Conclusions 

We conclude this review by providing a list of open questions in the areas of ultracold 
quantum gases, QCD plasmas, and holographic dualities. This list is far from complete, 
and the selection of problems is guided by the topics discussed in the main body of 
the review. Many of these problems will benefit from the connections between the 
physics of ultracold quantum gases, the QGP, and holographic duality that are being 
explored in this Focus Issue. For example, improved studies of the hydrodynamic 
behavior of ultracold gases will benefit from what has been learned about second 
order hydrodynamics in the context of the QGP and the AdS/CFT correspondence. 
Holographic dualities will also continue to play an important role in studying the 
approach from the initial, far-from-equilibrium, stage of a heavy ion collisions to the 
hydrodynamic regime. Far-from-equilibrium states have also been explored in ultracold 
atomic gases. Finally, new experimental results on properties of the QGP obtained in 
the ongoing programs at RHIC and LHC, as well as new results from experiments with 
trapped gases may point to new, unexpected, connections between these fascinating 
systems. 
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5.1. Open Problems and Questions in Ultracold Quantum Gases 

(i) What is the structure of the unitary Fermi gas in both two and three dimensions? 
In particular, is it possible to understand properties of the system in terms of well 
defined quasiparticles? 

At high temperature the gas is always weakly correlated, independent of the 
strength of the interaction. In the weak coupling regime we know that the 
quasiparticles near T c are BCS-like gapped fermions, and that at very low 
temperature the quasi-particles are phonons. At unitarity it has been argued that 
the regime T c < T < 1.5T C can be understood in terms of a pseudogap [404, 405, 35]. 
This means that there is a gap in the single particle spectrum, but no long range 
coherence. Pseudogap behavior in three dimensions has been observed using radio 
frequency (RF) spectroscopy [406, 407], and it has also been seen in quantum Monte 
Carlo calculations [408]. Pseudogap behavior was also reported in experiments with 
two-dimensional systems [409]. However, it is not entirely clear to what extent 
thermodynamic and transport properties can be understood in terms of gapped 
quasiparticles. One quasiparticle that may turn out to be relevant is the polaron, 
a single spin up particle immersed in a sea of spin down fermions [410]. In highly 
polarized gases RF spectroscopy has been used to measure the dispersion relation 
of polarons. It was observed that the polaron description can be extended to more 
spin balanced systems. In particular, a simple estimate based on a gas of polarons 
provides a very good estimate for the critical polarization at the Chandrasekhar- 
Clogston point, the first order transition from an unpolarized superfluid to a 
polarized normal gas [411]. Recently, polarons have also been observed in two- 
dimensional Fermi gases [412, 413]. 

(ii) How do we make local measurements of transport coefficients like shear and bulk 
viscosity, thermal conductivity and the spin diffusion constant? 
Measurements of elliptic flow and collective mode damping provide trap averaged 
values of the shear viscosity (see Sec. 2.5) but there is no direct measurement of 
the local value of the shear viscosity r)(n,T). The status is roughly analogous 
to the situation after the first experiments aimed at measuring the equation of 
state. These experiments only determined the total energy and entropy [66]. 
Local measurements were made possible by new ideas like the Gibbs-Duhem 
method employed by the ENS group [68] and the compressibility thermometer 
introduced by the MIT group [36]. It is possible in principle to unfold the currently 
available data for trap averaged viscosities, but this method requires a very careful 
treatment of the dilute corona, and a more direct method would be desirable. 
One possible direction is the development of new experimental techniques that 
generate local shear flows and study their relaxation. Another option is to measure 
the frequency dependent shear viscosity via spectroscopic methods as suggested 
in [414]. Similar questions apply to other transport properties, like bulk viscosity, 
thermal conductivity and spin diffusion. Current methods are mostly sensitive to 
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average transport coefficients, but more local measurements are needed. 

(iii) Are there reliable approaches to transport theory beyond the Boltzmann equation? 
In the unitary limit, transport properties can be computed reliably using the 
Boltzmann equation in the limit T ^> Tp. For T ~ Tp the quasiparticles are 
strongly interacting, and the kinetic description in terms of binary scattering 
between atoms breaks down. Enss et al. have computed the shear viscosity using 
a T-matrix resummation [129]. A similar approach is discussed by Guo et al. 
in this Focus Issue [110]. These methods are quite successful in describing the 
thermodynamics, but there is no expansion parameter, and it is not clear if non- 
equilibrium properties are predicted as well as equilibrium features. In Sec. 2.2 we 
briefly discuss the application of 1/N and e expansions to equilibrium properties 
of the unitary gas, but except for the recent work of [415] these methods have 
not been extended to non-equilibrium properties. Similar remarks apply to other 
methods as well: There are detailed studies of equilibrium properties using the exact 
renormalization group [416], but no corresponding studies of transport properties. 
Finally, there now exist very accurate quantum Monte Carlo calculations of the 
equation of state, but there is only a single, pioneering, study of transport properties 
[5]. 

(iv) Are there other universal quantum gases that can be created by manipulating the 
dimensionality, the spin structure, the type of interaction, or the dispersion relation? 
The unitary gas exhibits a very high degree of universality: in the limit of infinite 
scattering length and zero effective range the many body system is completely 
characterized by the fermion mass, the temperature, and the chemical potential. 
The dependence on dimensionful combinations of these parameters is completely 
governed by dimensional analysis. Universality can be exploited to use the unitary 
gas as a model system for dilute neutron matter; see Sec. 3.1. The question 
is whether this idea can be extended to more complicated systems, for example 
three-species gases that can serve as a model for the quark matter /nuclear matter 
transition (see Sec. 2.7), or four-species gases that correspond to self-bound nuclear 
matter, which is a liquid of equal numbers of spin up and down protons and 
neutrons. Universality in systems with more than two degrees of freedom is 
more complicated, because the interaction also depends on at least one three-body 
parameter [417], and the system may require repulsive short range forces in order to 
prevent collapse. New universal states can also be created by combining Bose and 
Fermi gases in different dimensions [418], or by manipulating the dispersion relation 
of the atoms. Steps in this direction were recently taken by Salger et al., who 
demonstrated linear dispersion and Klein tunneling of a Bose-Einstein condensate 
in a one- dimensional optical lattice [419], and by Taruell et al, who created a two- 
dimensional optical lattice with the same honeycomb lattice structure as graphene, 
and therefore obtained low energy fermions with a linear dispersion relation [420]. 

(v) What is the physics of the crossover in the discrete context of optical lattices? 
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Will the predictions of holographic duality of a lower limit on the ratio of viscosity 
of entropy hold true in the lattice context, or not? What class of theories from 
holographic duality might map onto the Fermi Resonance Hamiltonian or other 
effective two-channel models? How can we extend such models to treat the 
imbalanced Fermi gas? How will a polaron gas be described in the lattice context? 

(vi) What will be the interplay between disorder and interactions for the unitary gas? 
The relative effects of disorder and interactions on both thermal and quantum 
phase transitions in quantum many body physics is an open question generally 
(see e.g. [421]). An essential feature of disorder is Anderson localization, in which 
scattering off many small defects conspires to localize all states above a certain 
threshold in 3D. This concept is in fact a wave physics or single particle quantum 
concept: it is so far unclear where and how interactions destroy or allow for 
Anderson localization in a generic sense. Disorder was treated in this Focus Issue 
for ultracold Fermi gases specifically by Han and Sa de Melo [172]. Here, too, we 
find potential connections to holographic duality, as explored in recent forays into 
disordered systems by Adams and Yaida [422, 423]. Although these explorations 
have not treated lattice physics particularly, nevertheless disorder is an essential 
feature of lattice systems arising in nature, and can be induced artificially in various 
ways for optical lattices [424]. Can holographic duality capture the underlying 
lattice causes of disorder? How will experiments in ultracold Fermi gases combine 
disordered optical lattices with unitary fermions at sufficiently low temperatures? 
Might studies in this direction gain us some fundamental insight into the interplay 
between disorder and interactions in strongly correlated systems? 

5.2. Open Problems and Questions in Quantum Chromodynamic Plasmas 

(i) Can we determine all the transport properties of the QGP? 

We would like to determine, with fully quantified experimental and theoretical 
uncertainties, the value of the shear viscosity to entropy density ratio rj/s, the heavy 
quark diffusion constant D, and the jet quenching parameter q. This program is 
probably closest to completion in the case of t]/s. In this case the basic method, 
second order relativistic viscous hydrodynamics, is now fairly well understood. 
Final state corrections can be handled by coupling to a hadronic cascade. There 
is a large set of data, including the energy, impact parameter, and system size 
dependence of elliptic flow, as well as the observation of higher harmonics that 
points to the validity of the hydrodynamic description. The internal consistency of 
hydrodynamics already implies a fairly strong bound rj/s ^ (0.5 — 1.0). The main 
problem is that we have no direct method for establishing the initial conditions 
for the hydrodynamic description. This problem may be addressed by focusing on 
a larger set of observables, in particular higher harmonics of the flow [425, 426]. 
There also remain theoretical uncertainties regarding the coupling of viscous fluid 
dynamics to a Boltzmann description, the role of bulk viscosity, and the equation 
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of state. A discussion of the known uncertainties in rj/s can be found in [427]. 

The main strategy for determining the heavy quark diffusion constant is to 
use a Langevin description of the diffusion process where the background collective 
expansion is taken from hydrodynamics. The most important problem at present 
is the fact that there are only data for single leptons, and not for identified heavy 
quarks. This problem will be solved by improved detection capabilities at RHIC 
and LHC. The main theoretical issue is how to handle the interplay between elastic 
and inelastic scattering. For very heavy quarks elastic scattering is dominant, 
but for charm quarks inelastic processes may be important. This also includes 
hadronization effects such as coalescence. 

In the case of jet energy loss there are still a number of theoretical issues that 
need to be resolved. The radiative energy loss of an energetic parton (E 3> T) in 
a partonic medium can be expressed in terms of a single medium parameter, the 
transverse diffusion constant q. It is not clear if this carries over to the spectra of jets 
or leading particles produced in the medium. Another source of uncertainty is that 
different implementations of the in-medium radiation in perturbative QCD lead to 
significantly different values of q when applied to the RHIC data [278]. There is a 
significant ongoing effort to understand these differences [428]. On the experimental 
side the LHC represents a large step towards higher jet energies E ~ 100 GeV, 
and significantly improved statistics in the range E < 20 GeV already explored at 
RHIC. Both at the LHC and RHIC experimentalists are studying not only spectra 
of hard particles, but also identified jets, jet shapes, and correlations between hard 
particles and the associated soft emission. 

(ii) Can we understand early thermalization? 

The success of the hydrodynamic description implies that the system must 
thermalize early. The question is how to quantify this statement experimentally, 
and how to understand early equilibration theoretically. The most reliable 
constraint on the thermalization time arises from the observation of elliptic flow. 
Elliptic flow is driven by the anisotropy of the initial state. Since ballistic expansion 
dilutes the initial anisotropy local equilibration must happen early for elliptic 
flow to develop. A very conservative constraint, discussed in [182], is tq < 2.5 
fm/c. Typical values used in hydrodynamic fits are t < 1 fm/c. Another 
hint of early thermalization comes from the observation of thermal photons with 
temperatures significantly larger than those seen in hadronic spectra. The PHENIX 
collaboration has recently reported a photon transverse momentum spectrum with 
a slope parameter T = 221 ± 19 ± 19 MeV [429]. The spectra can be reproduced 
in hydrodynamic calculations with thermalization times ranging from 7~o ^ 0.15 
fm/c to t < 0.6 fm/c, corresponding to initial temperatures T ~ (300 — 600) 
MeV. Establishing experimental constraints on the equilibration time is further 
complicated by the fact that not all observables require full thermalization. For 
example, for the development of radial and elliptic flow it is not necessary that 
longitudinal momentum distributions are fully equilibrated [430]. 
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Understanding thermalization theoretically is difficult, even if the early evolution 
is governed by weak coupling. The first attempt to understand thermalization 
starting from an initial color-glass state using 2—7-2 and 2—7-3 scattering is the 
bottom-up scenario described in [431]. It was later realized that this picture is 
modified by collective effects that lead to plasma instabilities [432, 433, 434]; see 
also the contribution by Dusling in this Focus Issue [435]. In all of these approaches 
it is hard to understand how thermalization at RHIC can occur on a time scale 
t < 1 fm/c. 

Attempts to understand thermalization at strong coupling have focused on 
holographic duality; see Sec. 4.3.6 and [436, 437] for recent overviews. Holographic 
duality provides an interesting geometric picture of thermalization: the initial 
state corresponds to two colliding shock waves, and thermalization is signaled by 
the formation of an event horizon. The hydrodynamic stage is described by the 
relaxation, or ringdown of the black hole. In holographic duality one can achieve fast 
equilibration [312], but it is hard to make contact with asymptotic freedom and the 
well-established theory and phenomenology of nuclear deep inelastic scattering, and 
it is difficult to understand the observed energy and impact parameter dependence 
of the total multiplicity [438] . 

(iii) Does the quark-gluon plasma have a quasiparticle description? 

We would like to understand how the properties of the hot and dense matter 
produced at RHIC and the LHC are related to its structure, in particular whether 
the matter can indeed be described as a plasma made of quarks and gluons. 
Holographic duality would seem to indicate that a strongly interacting quantum 
fluid near the viscosity bound does not have a quasiparticle description; see Sec. 3.2 
and Sec. 4. In Sec. 3.2 we discussed a number of observables that have been 
suggested as probes of the quasiparticle structure: fluctuations, the observed 
constituent quark scaling of flow, and the relation between the heavy quark diffusion 
constant and the shear viscosity. None of these observables provides a conclusive 
test of the quasiparticle picture, and new ideas would certainly be helpful. We also 
discussed the possibility to study quasiparticles by extracting the spectral functions 
of the correlators of conserved charges in lattice QCD, such as energy, momentum, 
flavor, etc. These studies, too, are somewhat indirect because on the lattice we 
only have direct access to imaginary time correlation functions. This means that 
the spectral functions have to be extracted via analytic continuation. Analytic 
continuation of numerical data is difficult, but there has been significant progress 
over the last couple of years [439]. 

(iv) Can we experimentally locate the phase transition? 

There is a continuing effort to establish the presence of a critical endpoint in the 
QCD phase diagram. The basic idea is to look for non-monotonic behavior of 
fluctuations or correlations as a function of the energy of the colliding system. The 
first data from the RHIC beam energy scan have recently been released [440]; see 
also the contribution by Mohanty in this Focus Issue [441]. The data contains new 
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information about the beam energy dependence of hydrodynamic flow which will 
help to establish the onset of nearly perfect fluidity. Gupta et al. have argued 
that data on fluctuations can be directly compared to lattice QCD, and that it can 
be used to set the scale for the phase transition at zero chemical potential [214]. 
There is not yet conclusive evidence for the existence of a critical point, but data for 
additional beam energies and with better statistics will become available in the near 
future. In heavy ion collisions the growth of the correlation length is limited by finite 
size and finite time effects, and it is important to identify observables that provide 
the best sensitivity to critical behavior. Observables that have been considered 
include fluctuations of conserved charges, like baryon number, other fluctuation 
observables, like the pion-to-kaon ratio, the mean px, or elliptic flow, and higher 
moments like skewness and curtosis [442]. There is a parallel effort to locate the 
critical endpoint in lattice QCD. This effort is complicated by the sign problem 
that affects lattice calculations for fermions at finite chemical potential [443]. 

The discovery of a critical point would establish the presence of a phase 
transition between hadronic matter and the QGP. At small chemical potential this 
transition is only a crossover. One might hope that it would be possible to detect 
the crossover through hydrodynamic effects associated with the minimum in the 
speed of sound near the crossover temperature [444]. However, explicit calculations 
seem to show that there is no direct link between the presence of a minimum in the 
speed of sound and the beam energy dependence of multiplicity and flow [182]. 

Finally, we would like to demonstrate the existence of quark matter at high 
density and low temperature. It is not known whether the cores of compact stars are 
sufficiently dense to contain quark matter or other exotic phases. The presence of 
quark matter can potentially be detected through its effect on the equation of state, 
which is reflected in the mass-radius relationship. There is a very active program to 
pin down this relationship using observations of compact stars [445, 446]. Once the 
presence of a high density phase has been established, more detailed information 
regarding its properties can be derived using the cooling and spin-down behavior 
of compact stars. 

5.3. Open Problems and Questions in Holographic Duality 

(i) What does it take to derive a holographic duality? 

We have no first-principles derivation of holographic duality. This is perhaps not 
surprising, since it is a strong-weak duality, but it is frustrating. We do have a 
number of heuristic arguments, such as appeared in Sec. 4.1 above. We also know 
how to derive a long but finite list of dual pairs as decoupling limits of particular 
string theories, e.g., the original construction, A^=4 Yang-Mills in 4d. But we do 
not know how to derive the duality directly within the bulk or boundary theories.* 
This is important for two main reasons. First, we do not have a thorough 

* This is not for lack of trying; see[333, 447, 448] for recent approaches. 
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understanding of the consistency conditions or regimes of validity of our holographic 
dual descriptions. We do have general rules of thumb, such as the scalings discussed 
in Sec. 4.2.3, and we can check for self-consistency in specific examples. But 
particularly when we are defining a QFT via the holographic prescription, it is 
not always obvious what the precise regime of validity is, and to what extent we 
should trust the chosen gravitational truncation. 

Second, if we want to generalize holographic duality, it would be very helpful 
to know why it is true in the more familiar examples. For instance, we would like 
to generalize to intrinsically anisotropic systems, to asymptotically flat or positive- 
curvature (dS) spacetimes, or even to understand why this is not possible. 

(ii) Given a QFT, what is the gravitational dual? 

Even if a general derivation is not possible, it would be useful to have a constructive 
map which allowed us to take an arbitrary QFT and identify the appropriate 
gravitational dual. For example, what is the exact dual of QCD, not QCD plus a 
collection of exotic fields, not of large iV QCD, not of some toy model, but pure 
QCD? Similarly, with an eye towards ultracold atoms, what is the gravitational 
dual of the Hubbard model, or of the Kondo problem? 

Again, the fact that holography is a strong-weak duality should give us pause. 
As we have seen, when the gravitational system is weakly coupled, the dual QFT 
generically does not have any simple quasiparticle interpretation. If we want to 
define the QFT using conventional Lagrangian methods, the best we will in general 
be able to do is say that it is the strong-coupling limit of some explicit perturbative 
theory. In the models we understand best, the coupling that is getting strong 
is tightly constrained by symmetries, e.g. gauge symmetries, flavor symmetries, or 
supersymmetry, so we know how to follow the RG flow well beyond the perturbative, 
quasiparticle regime. We do not know how to solve the more general case without 
these extraordinary symmetries and the constraints they imply. 

(iii) Can we derive gravity from field theory? 

The holographic graviton looks very much like an emergent vector boson in more 
familiar condensed matter systems. It arises as an emergent dynamical gauge 
field which reorganizes the strongly-coupled physics into a set of perturbative 
local fluctuations. However, in holography the graviton (a) lives in one higher 
dimension and (b) couples universally. Point (a) fits naturally with a basic feature 
of gravity: since gauge transformations are coordinate transformations, there can be 
no truly local degrees of freedom, so any emergent graviton must come along with 
an emergent spacetime dimension with which to fix its gauge redundancy, much as 
any emergent vector boson must come along with an extra Hubbard-Stratanovich 
scalar to fix its gauge redundancy. This was in fact the origin of the term holography 
in gravitational physics. Point (b) is more spectacular from the point of view of 
emergence, but is in fact a necessary feature of any theory with a fluctuating spin-2 
gauge field. 

It is thus tempting to wonder whether we can in fact derive the graviton in a 
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general spacetime as an emergent gauge boson in a strongly interacting quantum 
field theory. Some progress has been made on this front (c.f. [447]), but a detailed 
understanding remains elusive. 

(iv) Do we always need the complete holographic description? 

In many of the examples discussed above, and indeed almost generically in the 
holographic literature, most of the work is being done by a very small fraction of 
the total degrees of freedom in the holographic description. For example, in the 
computation of r\j s in A/=4 SYM, whose holographic dual, string theory in AdSsX 
S 5 , includes an infinite number of dynamical massive fields, the only fields that 
matter in the computation are the 5d metric and a single scalar mode encoding 
one particular fluctuation of this metric. In practice, keeping track of the full 
holographic description often seems to be overkill. s 

This raises two natural questions. First, when does it make sense to truncate 
the full holographic description to a small subset of modes? Second, might this 
truncated holographic model have a more straightforward QFT representation 
than the full holographic description? Said differently, might there be a "semi- 
holographic" formalism defined entirely within field theory which allows us to 
build conventional QFTs which reproduce the interesting phenomenology of fully 
holographic QFTs without using higher dimensions? 

Recent work on holographic Fermi liquids, holographic quantum critical points 
and holographic quantum liquids (see [309] and especially the work of Nickel and 
Son in this Focus Issue [311]) demonstrates that at least some of the interesting 
holographic effects can indeed be realized in a purely quantum field theoretic 
formalism. On the other hand, for questions where the non-linear dynamics of the 
geometry such as far-from equilibrium physics, as opposed to kinematics or linear 
response, are important, it seems unlikely than any semi-holographic description 
can usefully capture the salient physics. Just how far a semi-holographic approach 
can be pushed, then, is an extremely interesting question. 

(v) What predictions does holography make for generic strongly coupled systems? 
Once we let go of the idea that our models must be defined by traditional 
quasiparticle effective QFTs, we can start asking general questions about the 
strongly coupled QFTs illuminated by our holographic examples. What common 
properties do these theories share? What general lessons can we draw about their 

s Note that this is analogously true in many-body physics: to describe ground states and low energy 
dynamics, one generally only needs a tiny corner of Hilbert space. However, in the holographic 
description the physics of a single bulk field represents an enormously complicated and non-local set of 
interactions amongst the boundary degrees of freedom. So observing that the bulk description can be 
effectively simplified to a small subset of bulk fields does not immediately imply that the same is true of 
the boundary QFT. For example, in holographic models of non-Fermi liquids, the bulk physics can be 
truncated rather dramatically, but the dual QFT has no quasiparticle description and so has not proven 
amenable to any simple truncation; concretely, the spectral function does not have quasiparticle poles. 
With that said, it seems well worthwhile to search for relations between known effective truncations in 
quantum many-body and holographic models. 
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kinematics, dynamics, phase structures, critical behavior, equilibration, etc? For 
example, we already know to expect nearly inviscid hydrodynamic behavior at finite 
density and low temperature. What else should we expect in general? Perhaps most 
importantly, what lessons can we draw for systems we really care about? 

(vi) What is the holographic dual of a non-relativistic QFT? 

As discussed in Sec. 4.3.4, non-trivial progress has been made in studying non- 
relativistic QFTs and many body systems. This has led to some simple heuristic 
lessons, for example that we should expect the viscosity to entropy ratio to again 
be extremely low, and that the Bertsch parameter^ should be expected to be of 
order unity and remain relatively insensitive to the coupling constant near strong 
coupling. However, the state of current models is rather poor, with awkward 
thermodynamics and bizarre phase structure. What are we missing? Can we build 
a holographic dual in the same universality class as fermions at unitarity? 

(vii) How (and why) is quantum information encoded in the bulk geometry? 
Holographically, information about the thermodynamic state of the boundary 
quantum field theory is encoded in the detailed geometry and matter field 
backgrounds of the bulk theory. For example, the temperature and chemical 
potential are encoded in the mass and charge of the bulk black hole. Similarly, 
when the QFT is in a pure state, the bulk geometry should encode the quantum 
correlations of the boundary QFT at zero temperature. But how? 

One approach to this question is the Ryu-Takanayagi conjecture [449, 450], 
which states that the entanglement entropy associated to a region R in a holographic 
QFT is given by the area of the minimal-area surface in the bulk which ends on the 
boundary of the region R. Holographically, this conjecture matches the relation 
between the thermodynamic entropy of a black hole and the area of its horizon. 
Interestingly, it suggests an intimate connection between quantum and thermal 
entropy - and more generally, between quantum and thermal fluctuations - all of 
which are encoded by the physics of horizons in the holographic dual. 

While a general proof remains lacking, considerable evidence has been found 
for this conjecture, and conversely numerous predictions have been made for the 
entanglement entropy of strongly-coupled many-body systems by exploiting this 
conjecture (see e.g. [451, 452, 453, 454, 455, 456] and references therein). For 
example, in this focus issue, Albash and Johnson [457] study the evolution of 
the Ryu-Takanayagi entropy under a thermal quench, leading both to evidence 
supporting the conjecture and to predictions for the evolution of the entanglement 
entropy. 

(viii) How does spatial inhomogeneity change the story? 

As we have seen, holography translates quantum computations in the boundary 

" We defined the Bertsch parameter of the unitary Fermi gas in Sec. 2.2. Here we use the term to 
refer to ratio of the ground state energy at strong coupling to that at weak coupling in a general 
non-relativistic CFT. 
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QFT (i.e. Feynman diagrams and path integrals) into classical GR calculations 
in the bulk (i.e. solving nonlinear partial differential equation boundary value 
problems). As such, it is clear why most research to date has focused on 
homogeneous, isotropic problems. However, from the quantum many-body side 
we know that inhomogeneity can radically alter the dynamics, and that the precise 
way they do so provides a detailed probe of the physics driving those dynamics. The 
same is true for anisotropy. Moreover, much of what we know about inhomogeneity 
and anisotropy in QFT comes from some form of perturbation theory. It would 
thus be particularly interesting to study the dynamics of holographic systems in 
the presence of inhomogeneity and anisotropy. Will we find evidence of many-body 
localization? Do (2+l)d QFTs at strong coupling and in the hydrodynamic regime 
display a canonical turbulent cascade, or an inverse cascade as is usually though 
not universally expected in two spatial dimensions? 

Many of the issues discussed in Sees. 5.1-5.3 are addressed in the contributions 
to this Focus Issue, and we hope that this summary, as well as the insights contained 
in these contributions, will spur further advances in the intriguing subject of strongly 
correlated quantum fluids. 
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